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ABSTRACT 


The  report  promotes  an  understanding  of  V/STOL  analysis  for  stimulation 
purposes  and  develops  equations  of  motion  compatible  to  either  analogue 
or  real  tins  digital  solution.  I 

A  general  set  of  equations  of  motion  are  developed  In  vhtch  axis  systems 
and  aerodynamic  coefficients  a^e  minimized.  Equations  of  motion  are 
then  developed  for  five  different  V/S7CI.  aircraft  wherein  thq^need  for 
additional  axis  systems  and  lero-fernamic  coefficients  for  a  particular 
V/STQL  configuration  is  developed. 
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VOLUME  II 

V/STOL  ANALYSIS  REPORT 


Foreword 


With  increasing  interest  in  Vert ica 1 /Short  \lakc  Off  and 
Landing  (V/STOL)  aircraft  by  the  military  services  it  was 
determined  that  it  wiuld  be  to  the  advantage  of  the  U.S.  Naval 
Training  Device  Center  to  eathcr  together  in  onc.rejort  a  study  of 
axis  systems  used  for  rcr resontat ion  of  such  aircraft,  and  to 
develop  "eneral  aerodynamic  equations  of  motion  which  nay  be  used 
in  simulation  of  these  aircraft.  This  effort  hasbeen  accomplished 
in  this  report, for  five  experimental  categories  of  V/STOL\aircraft , 
each  of  which,  as  of  this  date,  appear  to  give  promise  of  becoming 
operat iona 1 . 


Gamine  C.  Castellano 
Aerospace  Staff  Conultant 
U.S.  NaVal  Training  Device  Center 
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SECTION  I 
INTRODUCTION 

The  purpose  of  this  report  is  to  promcte  an  understanding  of 
V/STOL  analysis  for  simulation  purposes  and  to  develop  equations  of 
motion  in  a  form  most  compatible  to  either  analog  or  real-time  digital 
solution.  Only  V/STOL  aircraft  configurations  are  considered  in  this 
analysis. 

In  Section  II  of  the  report  a  method  of  attack  for  the  general 
V/STOL  problem  is  developed  and  various  axis  systems  are  selected. 

The  method  developed  is  then  applied  to  a  number  of  different  types 
of  V/STOL  aircraft.  A  model  and  its  conceptual  representation  are 
developed  in  Section  III.  This  model  may  be  used  as  a  guideline  in 
developing  equations  of  motion  for  various  types  of  V/STOL  aircraft 
but  cannot  be  further  simplified  in  general.  The  various  terms  in  the 
final  equations  for  a  particular  aircraft  conf iguration  must  be  ana¬ 
lyzed  for  the  magnitude  of  their  contributions  before  final  sinplifica- 
tior.  is  possible. 

The  XC-LU2A  tilt— wing  and  the  VZ-iiDA  tilt-duct  aircraft  are  used 
to  illustrate  equation  development  in  Section  IV.  Equations  are  pre¬ 
sented  though  not  painstakingly  developed  for  the  X-19  tilt-prop,  the 
XV -5a  fan- in-wing  and  the  P.1127  rotating  thrust.  A  great  deal  of  the 
symbology  used  for  these  different  V/STOL  aircraft  is  used  Just  as  pre¬ 
sented  by  the  respective  manufacturer  in  order  that  verifying  data  might 
be  available  in  the  proper  form  at  the  earliest  possible  time.  At  pre¬ 
sent  It  appears  that  the  equations  developed  will  c  onprise  a  miniaaim 
requirement.  Terms  in  these  equations  which  prove  in  flight  test  to  be 
sufficiently  small,  will  be  discarded.  Analysis  of  flight  test  data  may 
enable  simulation  with  less,  rigorous  equations  but  no  such  sijqsllflca- 
tlon  can  now  safely  be  made. 

A  bibliography  and  Appendix  A  dealing  with  basic  concepts  of  axis 
systems  and  equations  of  motion  have  been  provided  for  the  convenience 
of  the  reader. 


(1)  ftefer  to  JflLVTOCIEVCEJf  126^-1  tor  an  analysis  of  single  and  tandem 
rotor  helicopter  configurations. 
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SECTION  II 
AXIS  SYSTEMS 


METHOr  OF  SELECTION 

We  will  initially  consider  that  the  aircraft,  is  operating  out  of 
the  V/STOL  region.  This  aircraft  whether  it  be  a  transport,  observa¬ 
tion  plane,  subsonic  interceptor  or  any  similar  type  of  vehicle  can  be 
considered  to  be  operating  in  a  conventional  flight  regime.  A  descrip¬ 
tion  of  both  the  equations  of  motion  and  necessary  axis  systems  for  this 
mode  of  operation  has  been  done,  notably  by  Connelly  and  others  in  Refer¬ 
ences  2  and  7.  Connelly,  in  Reference  2,  discusses  simulation  of  fixed- 
wing  aircraft  and  presents  optimum  representations  of  the  equations  of 
motion  which  may  be  used  for  construction  of  a  computer.  Following  the 
recommendations  of  Reference  2,  we  will  use  wind  axes  for  the  descrip¬ 
tion  of  forces  and  aircraft  body  axes  for  the  description  of  moments 
in  the  conventional  flight  regime. 

Appendix  A  is  provided  as  an  aid  to  the  reader.  In  Appendix  A  the 
definitions  of  inertial  and  body  axes  are  developed.  Transformations 
between  inertial  and  body  axes  as  well  as  diagrams  illustrating  their 
relationship  to  one  another  are  shown. 

The  minimum  set  of  terms  in  the  equations  of  motion  of  a  V/STOL  air¬ 
craft  that  enable  a  description  of  its  conventional  mode  of  flight  will 
be  designated  as  a  necessary  description  of  the  particular  V/STOL  air¬ 
craft.  When  the  aircraft  is  in  such  a  normal  flight  attitude,  fewer  axis 
systems  and  aerodynamic  coefficients  may  be  needed  for  description  than 
during  hover  or  transition  flight*  For  some  V/STOL  configurations  the 
axis  systems  available  for  fixed-wing  aircraft  will  be  sufficient;  how¬ 
ever  for  other  V/STOL  conf igurations  additional  axes  must  be  defined. 

The  capability  of  a  V/STOL  aircraft  to  change  its  physical  configuration 
during  vertical  and  transitional  flight  may  indicate  the  need  of  addi¬ 
tional  axis  systems.  In  addition,  the  consideration  of  available  aero¬ 
dynamic  data  for  a  particular  vehicle  may  indicate  the  inclusion  or 
exclusion  of  certain  aerodynamic  coefficients  while  traversing  the  V/STCL 
region.  The  additional  terms  generated  by  consideration  of  the  V/STOL 
region,  together  with  the  minimum  set  of  terms  developed  by  consideration 
of  the  conventional  flight  regime,  will  constitute  a  necessary  and  suffi¬ 
cient  dec°rirtion  of  the  equations  of  motion  for  the  V/STOL  aircraft. 

It  is  to  be  expected  that  any  physical  change  in  the  actual  config¬ 
uration  of  the  aircraft  will  probably  necessitate  the  definition  of  addi¬ 
tional  axes,  but  the  existence  of  a  change  in  aircraft  conf iguration  or 
the  use  of  additional  thrust  producers  in  the  V/STOL  region  of  flight 
does  not  necessarily  force  the  development  of  new  axis  systems  in  addi¬ 
tion  to  those  already  specified  for  the  fixed-wing  aircraft.  Each  type 
of  aircraft  must  be  treated  individually  and  new  axes  developed  only  when 
the  effects  of  the  special  features  of  the  aircraft  configuration  cannot 
be  accurately,  economically  and  adequately  handled  with  the  axis  systems 
used  to  describe  the  conventional  mode  of  flight. 
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The  most  simple  case  of  configuration  change  would  he  that  causing 
physical  rotation  of  the  thrust  vector.  This  can  be  accomplished  by  tilt¬ 
ing  one  or  more  wings  at  the  same  or  different  rates  or  by  tilting  two  or 
more  ducted  fans,  propellers  or  rotors.  Thrust  rotation  can  also  be  ac¬ 
complished  by  deflecting  the  flow  from  a  thrust  source  such  as  a  ducted  fan. 
Yet  another  method  of  achieving  V/STOL  capability  which  may  reqxsire  special 
axis  system  treatment  is  the  use  of  special  thrust  producers  in  the  V/STOI. 
region  of  fl'ght  that  are  completely  separate  from  the  thrust  producers 
utilized  in  the  conventional  mode  of  flight.  These  separate  sources  of 
thrust  take  the  form  of  jet  engines,  fans,  or  completely  discrete  paths  for 
jet  exhaust. 

Possible  axis  systems  for  aircraft  configurations  involving  such  de¬ 
vices  are  presented  below, 

PARTICULAR  V/STOL  CONFIGURATIONS 

We  will  now  consider  some  examples  of  the  use  of  additional  axes  to 
describe  particular  V/STOL  aircraft.  In  the  final  consideration  of  a 
given  V/STOL  axis  system  the  form  of  coefficient  and  performance  data  avail¬ 
able  for  the  given  aircraft  will  contribute,  to  a  certain  extent,  to  the 
choice  of  axis  systems  for  simulation  of  that  aircraft. 

TILT-WING.  A  current  example  of  tilt-wing  aircraft  is  the  Tri -Service 
XC-lli2A  Transport  under  development  by  Ling-Temco-Vought,  Hiller  Aircraft 
and  Ryan  Aeronautical  Comp  any.  Figure  1  is  a  line  drawing  of  this  air¬ 
craft.  The  XC-1U2A  is  powered  by  four  G.E.  T-6h  turbo-prop  engines. 

These  engines  are  mounted  on  a  wing  which  can  be  tilted  to  an  angle  of 
100  degrees  from  the  x-y  plane  of  the  aircraft  body  axes  such  that  hover 
in  a  tail  wind  is  possible.  The  tail  assembly  is  conposed  of  a  vertical 
fin  rudder,  a  movable  horizontal  stabilizer  (no  elevator)  and  a  tail 
rotor  for  pitch  stabilization. 

In  order  to  reduce  the  problem  of  keeping  track  of  the  relative 
orientation  of  the  wing  chord  and  the  wind  as  well  as  the  orientation  of 
resulting  force  vectors,  a  wing  axis  system  is  defined.  This  system  has 
its  origin  at  the  wing  pivot  point  and  has  its  x-a  plane  coincident  with 
the  body  axis  x-z  plane.  The  system  rotates  with  the  wing  such  that  tho 
wing  x  axis  is  always  parallel  to  the  wing  chord.  Further,  a  set  of  wing 
stability  axes  is  defined  analogous  to  stability  axes  of  fixed-wing  air¬ 
craft.  This  system  is  rotated  about  the  wing  y  axis  by  the  angle  of  attack 
of  the  wingj  i.e.,  the  angle  formed  by  the  projection  of  the  aircraft 
velocity  vector  on  the  x-z  plane  of  the  wing  axis  system  and  the  wing  x  axis. 

Also,  in  the  transition  and  hover  regions  the  forces  developed  at 
each  propeller  become  the  primary  components  of  the  stability  of  the 
aircraft.  Therefore,  an  axis  system  (the  propeller  axis  system)  is  speci¬ 
fied  for  each  main  propeller.  It  my  also  be  necessary  tc  develop  a 
propeller  axis  system  for  the  tail  rotor  of  some  similar  configurations, 
but  this  is  not  necessary  for  the  XC-lJb2A  as  shown  in  subsequent  portions 
ef  this  reports - - - 
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Figure  1. 


Three-View  of  XC-1),?A 
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These  axis  systems  enable  the  development  of  force  and  moment  ex¬ 
pressions  for  the  wing,  main  propellers  and  the  tail  rotor.  The  result¬ 
ing  forces  and  moments  are  transformed  to  aircraft  body  axes  and  there 
incorporated  into  the  total  force  and  moment  equations.  These  axes  in 
conjunction  with  inertial,  body  and  wind  axes  would  romprise  a  sufficient 
description  of  axes  for  simulation  of  a  tilt-wing  transport, 

TILT-DHCT.  An  example  of  the  tilt-duct  is  the  VZ-uDA  of  DOAK  Aircraft 
Company.  This  is  a  small  aircraft  with  tiltable  ducted  fans,  one  at 
each  wing  tip.  The  tail  assembly  is  composed  of  a  vertical  fin  and 
rudder,  a  horizontal  stabilizer  and  reaction  controls.  As  for  the  tilt¬ 
wing  it  is  desirable  to  compute  the  direction  and  magnitude  of  the  thrust 
vectors  from  each  ducted  fan  in  order  to  simulate  a  trim  aircraft  in  the 
V/STOL  region.  In  order  to  accomplish  this  each  rotatable  ducted  fan 
should  contain  an  axis  system.  The  resulting  force  and  moments  are  trans¬ 
formed  to  the  aircraft  body  axes  by  considering  the  angle  of  tilt  of  the 
ducts  with  respect  to  the  wing.  Figure  2  is  a  line  drawing  of  this  air¬ 
craft. 

Another  example  of  the  tilt-duct  is  the  oroposed  h  ducted  fan  VTOL 
transport,  the  Pell  D2127  of  Pell  Aerosystems Company.  Here  again  axes 
for  duct  velocities  and  duct  thrust  vectors  «*re  desirable.  Attached  to 
the  aft  side  of  each  duct  are  wing  areas  in  line  with  the  air  flow.  These 
wing  areas  rotate  with  the  ducts.  Forces  and  moments  created  by  these 
surfaces  can  be  computed  in  the  axis  system  defined  for  the  ducts,  and 
then  transformed  to  fuselage  axes  through  the  duct  tilt  angle,  thus 
enabling  summation  of  total  forces  and  moments.  In  general,  the  thrust 
output  of  a  ducted  fan  varies  with  its  direction  of  flight  with  respect 
to  axes  located  in  the  duct.  The  thrust  output  is  not  necessarily  col- 
linear  with  the  axis  of  the  fan.  The  components  of  force  perpendicular 
to  this  axis  become  significant  in  the  total  moment  equations  and  in 
general,  these  components  of  force  have  varying  lever  arms  on  which  to 
act  and  become  particularly  significant  terms  in  vertical  or  near  verti¬ 
cal  flight. 

TILT-PROP.  An  example  of  this  type  of  aircraft  is  the  Curtiss-Wright 
X-19,  Figure  3  is  a  line  drawing  of  this  aircraft  which  has  two  engines 
driving  four  propellers.  Four  tilting  propellers  are  mounted  in  nacelles 
at  the  tips  of  tandem  high  wings.  Stabilization  results  from  variance  in 
the  control  of  pitch  at  each  propeller.  Pitch  angle  control  comes  from 
simultaneous  change  in  the  pitch  of  the  forward  propellers  and  the  aft 
propellers  respectively.  Roll  angle  control  results  from  simultaneous 
change  in  the  pitch  of  starboard  propellers  and  the  port  propellers 
respectively.  Yaw  angle  control  is  developed  from  simultiuieous  change 
in  the  propeller  pitch  of  diagonally  opposite  propellers.  Since  the 
propeller  nacelles  rotate,  a  nacelle  incidence  angle  is  defined  with 
respect  to  the  aircraft  body  axes.  In  order  to  account  for  the  variance 
in  magnitude  and  direction  of  the  forces  developed  at  each  propeller 
individual  axis  systems  will  be  defined.  The  propeller  axis  systems  in 
addition  to  the  axis  systems  of  the  conventional  aircraft  enable  a  suf¬ 
ficient  description  of  the  equations  of  motion  of  the  X-19.  Propeller 
axes  are  required  because  of  the  importance  of  nonaxial  components  of 
thrust  in  each  nacelle. 
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Figure  3.  Three- View  of  X-19* 
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'•’V.'-  I '  -WIN",  An  example  of  thin  type  of  aircraft  in  the  Oenoral  Electric/ 
Evan  Aeronaut  leal  XV-rA  V/STOL  research  aircraft.  Figure  h  is  a  line 
drawing  of  ;nch  an  aircraft.  For  this  aircraft  exhaust  gases  are  emitted 
straight  out  the  tail  pipes  in  normal  flight.  "'or  transition  to,  and 
operation  in,  vertical  flirht  the  tailpipes  are  blocked  off  and  the  ex¬ 
haust  rases  are  diverted  to  three  free  wheeling  fans— one  in  the  nose  and 
a  pair  of  large  lift  fans  embedded  in  the  wings.  The  thrust  developed 
out  of  each  fan  is  controlled  by  fan  exit  louvers.  The  thrust  is  propor¬ 
tional  to  fan  rrn,  and  its  direction  is  controlled  by  the  setting  of  the 
louvers.  If  thrust  versus  fan  rpm  data  and  louver  positirn  versus  thrust 
data  were  available  such  that  the  direction  and  marnitude  of  the  thrust 
at  c?ach  fan  could  be  described,  then  forces  and  moments  could  be  developed 
in  existing  wind  and  body  axis  systems  so  that  the  aircraft  could  be 
simulated  by  methods  already  developed  for  subsonic  jet  aircraft.  How¬ 
ever,  if  necessary,  an  axis  system  describinp  louver  position  could  be 
developed  in  order  to  keep  track  of  fan  thrust  durinp  transition  and  hover 
conditions. 

ROTATING  THRi'ST  At"''  TTEFLFOTEf  THRHST.  An  example  of  rotating  thrust  is 
the  Hawker  P.1127  t’-unsonic  V/STOL  strike  fighter.  Fipure  5  is  a  line 
drawinr  of  such  an  a.rcraft.  A  deflected  thrust  aircraft  example  is 
the  VTPL  trainer  rt  aircraft,  X-liuA,  by  Hell  Aerosystems  Company. 

The  P.1127  is  powered  by  the  Rristcl  Siddeley  BS  53  engine  which  has  four 
simultaneous  rotating  nozzles.  The  thrust  being  ejected  from  each  nozzle 
is  capable  of  being  rotated  through  an  angle  of  100  degrees  from  a  posi¬ 
tive  X-aircraft  body  axis.  Nozzle  deflection  angles  over  90  degrees 
allow  for  braking  flight  during  a  STOL  maneuver.  For  this  system  we 
keep  track  of  the  thrust  vectors  from  the  exhaust  nozzles  of  the  power- 
plant.  This  can  be  done  by  tracking  the  rotation  of  the  main  exhaust 
nozzles  with  respect  to  body  axes.  No  extra  axis  systems  appear  to  be 
necessary. 

In  the  case  of  the  X-lUA,  thrust  diverter  position  and  power  output 
of  the  engines  may  give  enough  information  to  determine  thrust  magnitude 
and  orientation. 

For  both  the  P.1127  and  Xlli-A,  no  additional  axis  systems  are  neces¬ 
sary  for  the  development  of  the  equations  of  motion  for  the  purpose  of 
simulation.  Aerodynamic  data  should  enable  thrust  calculations  in  wind 
or  body  axes,  and  the  consequent  direct  conr: tation  of  aerodynamic  forces 
and  moments. 


-I 


Fipure  it. 


Th  ree-  V  lew 


f  XV-  ;A. 
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SECTION  III 
MATHEMATICAL  MODEL 

In  this  section  we  will  construct  a  mathematical  model  which  can  be 
applied  to  the  simulation  of  any  V/STOL  aircraft.  With  this  model, 
equations  of  motion  for  various  V/STOL  aircraft  configurations  are  evolved 
in  Section  IV.  The  approach  used  in  deriving  these  equations  of  motion 
will  enable  the  reader  to  develop  simulation  equations— as  needed— in  a 
systematic  and  simple  manner  for  V/STOL  conf igurations  that  are  evolved  in 
the  future. 

The  method  of  attack  is  to  start  with  the  equations  of  motion  repre¬ 
senting  the  conventional  mode  of  operation  of  the  V/STOL  aircraft;  to 
choose  appropriate  axis  systems  based  upon  the  particular  aircraft  which 
are  in  addition  to  the  inertial,  body,  stability  and  wind  axes;  to  deter¬ 
mine  the  V/STOL  aerodynamic  characteristics  not  simulated  by  the  aircraft 
in  its  conventional  flight  mode;  and  finally  to  write  the  equations  of 
motion  for  the  V/STOL  aircraft.  The  first  objective  is  to  develop  the 
idea  of  a  conventional  flight  mode. 

CONVENTIONAL  FLIGHT  MODE  OF  V/STOL  AIRCRAFT 

The  V/STOL  aircraft  (in  the  conventional  mode  of  flight)  is  de¬ 
scribed  by  a  set  of  necessary  equations.  There  is  no  attest  to  define 
whether  they  apply  to  single  or  multi-engine  propeller  or  Jet  aircraft. 

The  configuration  of  the  aircraft  is  not  defined  other  than  there  are 
control  inputs  for  rudder,  elevator,  flap  and  aileron.  There  is  no 
indication  as  to  low  or  high  performance  or  mission.  Modifications  of 
the  basic  simulation  equations  describing  this  model  will  depend  upon 
the  particular  V/STOL  aircraft.  Consequently,  specifying  the  particular 
aircraft  enables  the  addition,  if  warranted,  of  the  factors  which  allow 
the  necessary  equations  of  the  aircraft  operating  in  the  conventional 
flight  mode  to  become  a  necessary  and  sufficient  set  of  simulation  equa¬ 
tions.  In  conjunction  with  the  following,  the  development  in  Appendix  A, 
and  Reference  2,  we  are  able  to  write  the  equations  for  a  V/STOL  aircraft 
in  t.ie  conventional  flight  mode. 

MATHEMATICAL  MODEL  OF  CONVENTIONAL  FLIGHT  MODE  OF  V/STOL  AIRCRAFT.  The 
com 'ntional  flight  mode  is  represented  by  a  set  of  three  force  and 
three  torque  equations.  Linear  and  angular  rates  are  obtained  from 
these  equations. 

Linear  Motion.  The  following  are  the  force  expressions^^. 

FNi  -  F*11  (i  -  1,  2,  3) 

e  1  e  •  external 

I  •  internal 
N  -  inertial  axes 


HT  "See  Appendix  B  of  mm\mm  IV 05-1  and 

(2)  Connelly,  reference  2,  develops  similar  equations  of  motion. 
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r  Ni 

The  >F  are  the  total  external  forces  which  are  composed  of  aero- 

dynamic  and  gravity  forces  and  the  Jr  ^  are  the  total  internal  forces. 

We  can  view  the  aircraft  as  a  system  of  mass  points  such  that  the  dis¬ 
tance  between  pairs  of  mass  points  are  constant  throughout  the  aircraft’s 
motion.  This  is  a  rigid  body,  and  consequently  work  done  by  internal 
forces  between  pairs  of  mass  points  vanishes.  However,  since  we  have  a 
rigid  body  the  center  of  mass  of  the  mass  points  (m^)  is  concentrated  at 

a  single  point  and  the  radius  vectors  (r^)  from  the  inertial  frame  to 

each  arbitrary  mass  point  are  replaced  by  a  radius  vector  r  to  the  center 
of  mass.  We  then  have  the  following: 


Equation 


-  j>y?f  -  M  rN1 

M  -  Jmj  maos  of  the  aircraft 

(3.1) 

is  expanded  to  give: 

P? 

P? 

- 

(3.1) 


(3.2) 


To  facilitate  the  description  of  forces  aircraft  body  axes  are 
located  at  the  center  of  mass  of  the  aircraft.  Forces  described  in  the 
axes  are  transformed  to  the  X^  body  axes  as  shown  in  Appendix  A.  We 

have  then  the  following  equations: 


Jl^1  -  M(U  ♦  Wqx  -  V  r) 

£f£2  -  M(V  +  U  r  -  W  p)  (3.3) 

Jf£3  -  M(W  ♦  V  p  -  U  qx) 

In  equations  (3.3)  the  following  identifications  are  made: 

Jf**  ■  Aerodynamic  Forces  ♦  Gravity  Force  ♦  Engine  Thrust 
where  XI  ■  x,  X2  -  y,  X3  -  *,  T  -  engine  thrust,  and  o,p  is  the  angle 
of  the  thrust  line  with  respect  to  the  x-body  axis. 


I 


i 
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So  that  there  is 


X  -  M  g  sin  9  ♦  T  cos  aT 

at  1 

Y  ♦  M  g  cos  9  sin  $ 

a  (3.U) 

Z  +  H  c  cos  9  cos  4>  -  T  sin  a 
a  t 

and  (?.h)  the  force  equations  in  aircraft  body  axes  are? 

X  ■  M  (U  ♦  Wq  -  V  r)  ♦Mg  sin  0  -  T  cos  a_ 

ot  X  1 

Ya  -  M  (V  ♦  U  r  -  W  p)  -  M  g  cos  0  sin  *  (3.5) 

fll 

Za  -  M  (W  ♦  V  p  -  Uq^)  -Mg  cos  fe  cos  <f>  +  T  sin  a_ 

Angular  Motion.  Following  are  the  tcrque  express’  ons  for  the  aircraft 

in  conventional  flight.  It  is  practically  a  universal  practice  to 
develop  aircraft  angular  motion  in  body  axes  both  in  the  aircraft  manu¬ 
facturing  and  simulator  industry.  For  simulation  purposes  this  practice 
avoids  the  added  problem  of  transformations  of  the  products  and  moments 
of  inertia  which  would  be  required  if  other  axis  systems  were  used.  From 
Appendix  A  we  write  for  the  torque  components  the  usual  dynamic  terms  in¬ 
volving  products  and  moments  of  Inertia,  angular  velocities  and  accelera¬ 
tions,  and  the  coupling  of  angular  velocities. 

La  "  hi  P  *  I13  *  J>V  *  (I33  '  W  V 

Ma  •  I22  ^  -  I13  (r2  -  p2)  ♦  (In  -  Ijj)  pr  (3.6) 

Na  '  :33  ^  ‘  J13  (P  *  qlr)  +  (I?2  -  Ill)  pql 

In  the  equations  of  3.6,  the  x-z  plane  of  the  aircraft  body  axes  is  a 
plane  of  symmetry.  Consequently  products  of  inertia  (I^ ,  1^^,  I^2) 

are  zero.  L  ,  M  and  N  are  the  aerodynamic  coefficients  necessary  to 

ol  fil  A 

describe  ary  V/STOL  aircraft  in  conventional  flight. 


From  (3.3) 
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COEFFICIENT  FORCE  OR  70RCTTE  EXPRESSION 


X 

a 

Cx  (a,  Ma) 

i  o  V2SCx  (a,  Ma) 

1 

Cx(6) 

l  0  V2SC  (( 8) 

'xtF 

12 

1  0  V*SC  .  6F 

2  *6F 

Ya 

c 

y8 

f  °  vHb  8  v 

za 

(a,  Ma) 

\  0  V2SCa  (a,  Ma) 

1  P  V2SC  .  6F 

2  Z6F 

°Z6E 

1  0  V2SC  C  .  6E 

2  m6E 

C1 

1 8 

y  0  V2Sb  C.  *  8 

2  18 

C1 

\  o  V2Sb  Cn  .  6A 

2  i6A 

C1 

oR 

1  o  V2  Sbic.  .  6R 

2  \  26R 

C1 

16R 

1  0  V2  Sb  C.  .  6R 

2  oR 

C1 

P 

p  p  VSb2  C.  .  p 

U  P 

\ 

r  o  VSb2  C,  .  r 

a  s. 

Table  1. 

Hine ns  cn  less  Coeff 

icents  Used  for  Basic  Aircraft. 

Reproduced  From 
*  4-  Best  Available  Copy 
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COEFFICIENTS  FORCE  OR  TORQUE  EXPRESSION 


Ha 

C  (a,  Ma) 
m  9 

y  P  V2  SC  C  (a,  Ma) 

d  m 

C 

1  0  V2  Sc  C  .  6F 

2  "6F 

C 

1  p  V2  Sc  C  .  6E 

m6E 

2  "SE 

Cm 

i»VSo2C  .  qx 

qi 

ql 

Cm- 

a 

J  p  V  Sc2  C  .  « 

M  a 

N 

C 

1  2 

i  „  V£Sb  C  .  8 

a 

n6 

2  p  "s 

C 

\  P  V2Sb  C  .  6A 

n6A 

2  n6A 

c 

i  p  V2  Sb  C  .  6R 

n6R 

2  "6R 

% 

r  p  V  Sb2  C  .  p 

t ' 

C 

n 

C  »  V  **  Cn  *  r 

r 

-  r 

V  =  Aircraft  velocity 
p  =  Air  density 
S  =  Wing  area 
b  *  Wing  span 
c  =  Mean  aerodynamic  chord 

Table  1.  Dimensionless  Coefficients  Used  for  Basic  Aircraft  (Cont'd.) 


15 


NAVTRADEVCEN  1205-2 


EXAMT NATION  OF  COEFFICIENTS  TOR  V/STOL  AIRCRAFT  IN  CONVENTIONAL  FLIGHT 

MODE.  Following  the  work  presented  by  Connelly  in  Reference  2,  we  will  ( 

state  that  the  coefficients  listed  in  Table  1  constitute  a  necessary 

number  of  aerodynamic  coefficients  in  order  to  develop  the  aerodynamic 

forces  and  moments  for  a  V/STOL  aircraft  in  conventional  flight.  In 

many  cases  actual  manufacturers’  data  may  not  be  available  for  just  these 

specific  coefficients.  For  example,  data  are  usually  available  in  terms 

of  coefficient  of  lift  (C^)  or  drag  (C^)  which  in  turn  can  be  transformed 

to  give  Cx  and  C  .  Table  1  is  then  a  list  of  what  is  considered  to  be  a 

minimum  number  of  aerodynamic  coefficients  to  describe  the  aircraft  in 
conventional  flight.  Modification  of  this  list  is  possible  when  the 
flight  of  the  aircraft  is  described  in  the  V/STOL  region. 

For  an  actual  V/STOL  aircraft  the  equations  for  X,  Y,  Z,  L  ,  M 

A  fit 

and  Nft  may  contain  terms  due  to  the  specific  effect  of  various  aircraft 

components— the  wing,  propeller(s)  or  ducted  fan,  the  fuselage,  vertical 
tail  and/rudder,  horizontal  tail,  and  devices  to  insure  trim  flight  such 
as  a  small  variable  direction  thrust  propeller  or  ducted  fan.  Each  of 
these  forces  and  moments  due  to  the  various  aircraft  components  will 
have  associated  aero-Vnamic  coefficients— some  common  to  the  aircraft  in 
conventional  flight  and  other  aerodynamic  coefficients  applicable  only 
to  the  particular  V/STOL  aircraft. 

STATEMENT  OF  EQUATIONS  FOR  V/STOL  AIRCRAFT  IN  CONVENTIONAL  FLIGHT  MODE. 

From  the  force  equations  (3.5)  and  the  torque  equations  (3.6)  and  the 
data  of  Table  1  we  can  write  a  set  of  equations  for  the  aircraft  in 
conventional  flight.  It  must  be  remembered  that  these  equations  are 
only  guidelines.  They  give  us  an  insight  into  what  will  be  contained 
in  the  simulation  of  a  V/STOL  aircraft.  The  equations  which  follow 
(3.7  to  3.12)  are  written  in  aircraft  body  axes  and  do  not  include  ef¬ 
fects  due  to  landing  gear  or  other  variations  in  external  stores  or  the 
aircraft  configuration. 

X  Force  Equation 

\  P  V2S[C  (a,Ma)  +  C  ( £)  ♦  C  6F] 
c  6F 

m(U  +  Wq^  -  Vr)  +  mg  sin  0  -  T  cos 

2 

or  (U  ♦  Wqx  -  Vr)  -  [Cx(a,  ♦ 

T 

♦  —  cos  a_ 
in  a 

Y  Force  Equation 

o  V2S(yge  -  m(V  ♦  Ur  -  Wp)  -  mg  cos  0  sin  * 

•  %  v2S 

or  (v  ♦  Ur  -  Wp)  ■  p  C  .  6  ♦  g  cos  0  sin  $ 


(3.8) 


(3.7) 

aT 

C  C6 )  +  C  .  6F]  -  g  sin  0 
* 
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Z  Force  Equation 


or 


~  0  V2S  [C  (a,  Ma)  ♦  C  .  6F  ♦  C  .  6E] 

2  *  *6F  *6E 

m(W  ♦  Vp  -  Uq^)  -  mg  cos  0  cos  4>  *  ?  sin  a,p 

2 

(W  +  Vp  -  Uq.)  -  p  [C  (a,  Ma)  ♦  C  .  6F  ♦  C  .  6E) 
1  2m  z  *6F  Z6E 


+  g  cos  0  cos  ♦  -  ~  sin  a^, 


Roll  Equation 


+  bC,  •  r) 
r 


or 


(3.9) 


l  p  VSb(2VC,  .  0  +  2VC..  .  6A  ♦  2VC.  .  6R  ♦  bC-  .  p  (3.10) 

h  B  6A  6R  S 


“  TllP  *  Ii3  (*  *  P^)  ♦  (I33  -  I22)  qlr 

p  -  ♦  J1  (r  ♦  pq  )  -  —  ■  22  q^r 

*11  ill  * 

(2VC.  .  8  ♦  2VC-  ,  6A  ♦  2VC-  .  6R  ♦  b  C,  .  P  ♦  b  C,  .r) 

UI11  6A  16R  \  V 

Pitch  Equation 

1  p  VSc[2»Cm(a,  M.)  ♦  JVC^.  6F  *  ZTO^.  SE  *  c  C^.  (3.11) 

’  J22«l  ’  V*  '  pJ>  *  (I11  '  133)  >» 


or 


q  -  ♦  -il  (r2  -  p2) - - 

l  I22  I 


(•^n  "  ^33) 


22 


pr 


♦  o  rlS-  (27C  (a,Ma)  ♦  2VC  .  8F  ♦  2V0„  .  6E  ♦  cC  ,qn  -  cC.  .  a) 

hloo  n  1  m* 


22 


®8F 


"8E 


1? 
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Yaw  Equation 


r  0  VSb(2VC  .  8  ♦  2VC  .  6A  ♦  2VC  .  6R  ♦  bC  .  p 
U  P  n6A  n6R  np 

♦  bCn  •  r) 
r 

"  I33r  “  X13^  "  qlr^  *  ^22  “  111^  pql 


(3.12) 


or 


(X22  "  hl> 


i  ■  ♦  (p  -  V)  " 

*33  1  X33 


•f  P  rf—  (2VCL  •  8  ♦  2  VC  •  6A  +  2VC  *  6R  ♦  bC  .  p  ♦  bC  .  r) 
qI33  P  ^6A  ^6R  np  nr 

ADDITIONAL  AXIS  SYSTEM  SELECTION 

The  next  step  ir.  the  mathematical  model  is  to  choose  the  necessary 
additional  axis  systems  in  order  to  describe  the  particular  flight 
performance  of  a  given  V/STOL  aircraft.  In  Section  II  axis  systems  were 
discussed  in  conjunction  with  particular  V/STOL  configurations. 

V/STOL  AERODYNAMIC  COEFFICIENTS 

Once  a  system  of  axes  is  selected,  additional  coefficients  may  be 
specified.  These  could  relate  to  thrust  for  Jet,  propeller  of  ducted 
fan.  They  could  be  functions  for  speed  brakes,  slats,  external  stores, 
landing  gear,  drop  tanks,  bomb  bay  doors  or  louvers  for  control  of  air 
flow.  The  Inclusion  of  any  such  coefficients  is  best  illustrated  by 
the  tilt-wing  exanple  in  Section  IV. 

V/STOL  EQUATIONS  OF  MOTION 

With  the  equations  of  the  V/STOL  aircraft  in  conventional  flight, 
axis  systems,  and  a  study  of  the  particular  V/STOL  aerodynamics  dis¬ 
cussed  up  to  this  point,  we  are  now  reacty  to  write  the  equations  of 
motion  for  the  V/STOL  aircraft.  No  small  angle,  or  small  term  simpli¬ 
fications  have  been  made  in  the  force  or  torque  equations  of  the  air¬ 
craft  in  conventional  flight.  Siiqplifications  may  be  made  during  the 
development  of  the  equations  for  a  particular  aircraft. 
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SECTION  IV 

EQUATIONS  OF  MOTION  FOR  V/STOL  AIRCRAFT 

In  this  section  equations  of  motion  will  he  developed  for  five 
V/STOL  configurations.  These  configurations — tilt-wing,  tilt-duct, 
tilt-prop,  fan-in-wing  and  rotating  thrust— incorporate  individually 
and  collectively  the  designs  of  V/STOL  aircraft  currently  being  devel¬ 
oped.  A  particular  aircraft,  either  one  that  is  proposed  or  a  proto¬ 
type  will  be  used  as  an  example  of  each  of  these  configurations.  The 
first  type  to  be  considered  is  the  tilt-wing. 

TILT-WING 

The  Vought-Hiller-ftyan  XC-LU2A  Tri-Service  Transport  will  be  used 
in  the  example  in  the  development  of  tilt-wing  equations  of  motion. 

The  equations  from  Section  III  will  be  modified  so  that  they  are  appli¬ 
cable  to  an  aircraft  such  as  the  XC-lh2A. 

In  order  to  obtain  some  idea  of  the  XC-1U2A  let  us  consider 
Figure  1,  Table  2  and  Table  3.  In  Figure  1,  a  three  view  arrange¬ 
ment  of  the  XC-U12A  is  shown.  Table  2  contains  some  of  the  physical 
characteristics  of  the  aircraft  and  Table  3  is  a  list  of  definitions 
of  symbols  contained  in  the  XC-1U2A  equations. 

With  this  general  idea  of  the  XC-1U2A,  we  may  now  develop  the 
mathematical  model  for  the  XC-1L2A.  First,  define,  if  necessary, 
applicable  axis  systems.  Second,  develop  additional  aerodynamic  coef¬ 
ficients,  Finally,  write  the  equations  of  motion  for  the  XC-1U2A, 
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CHARACTER! STIC 
Normal  Gross  Weight 

Center  of  Gravity  (aft  of  the  leading  edge 
of  the  mean  aerodynamic  chord  (MAC)) 

Max  Forward 
Max  Aft 

Wing 

Total  Area 
Span 

Aspect  Ratio 
Dihedral  Angle 
Airfoil  Section 

Trail  Edge  Flaps  -  Double  slotted 
Maximum  Deflection 
Deflection  for  take  off  (STOL) 

Deflection  for  landing  (STOL) 

Leading  Edge  Flaps 
Deflection 

Ailerons  -  Plain 

Maximum  Deflection  (wing  up) 

Maximum  Deflection  (wing  down) 

Horizontal  Stabilizer  -  All  moving 
Area 
Span 

Aspect  Ratio 
Dihedral  Angle 
Airfoil  Section  (root) 

Airfoil  Section  (tip) 

Maximum  Deflection  (leading  edge  up) 
Maximum  Deflection  (leading  edge  down) 
Hinge  Line,  $  of  tail  mean  geometric  chord 

Vertical  Tail 

Area:  Fin  to  rudder  hinge 
Area:  Rudder  aft  of  hinge 
Aspect  Ratio 
Airfoil  Section  (root) 

Airfoil  Section  (tip) 

Fuselage 

Length 

Length  (including  tail  rotor) 

Outside  Height 
Outside  Width 

Maximum  cross-sectional  area 
Propellers 

Diameter 

Number  of  Blades 
Tall  Rotor 
Diameter 

Number  of  Blades 


VALUE 

37,L?L  pounds 


10$  MAC 
28$  MAC 


53L.37  ft2 
67.5  ft 
8.53 
-2.12° 

NACA  63-318  (Mod) 


37° 


♦  50° 

♦  15° 

163.5  ft2 
31.1L  ft 
5.08 
0° 

NACA  0015 
NACA  0012 
60° 

15° 

13$ 


95  ft2  2 
21.6  ft* 
1.87 

NACA  0018 
NACA  0012 


50  ft 
58.12  ft 
10.72  ft 
9.25,ft 
90  ft2 

15.5  ft 

L 

8  ft 

3 


Table  2.  Selected  Physical  Characteristics  1C-1L2A 
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Fll  Milage 


ENfir&tion 

1 


Wing 

i{ 

Main  propellers  n  ■  1#  2,3,  It  -  top  view 
left  to  right  \  ^ 

\ 

Vertical  Tail  '\ 

♦ 

Horisontal  Stabiliser 
Tail  Rotor 

V 

Angle  of  attack  -  fuselage 

\ 

Angle  of  attack  -  wing 
Angle  of  attack  -  horixontal  stabiliser 
Incidence  angle  -  horiaontal  stabiliser 
Downwaeh  angle 


Sideslip  angle  -  fuselage 
Sideslip  angle  -  wing 
Incidence  angle  -  wing 


«  "  Sr  * 


Total  relocity’  in  aircraft  body  axes 
Total  velocity  in  wind  stability  axas 


Mean  aerodyn^jnic  chord 
Hlng  span 

Total  area  of  wing 


Sp 

% 


Total  are a  of  rein  propeller  dick* 

Dmartc  pma«nra  due  to  p<v*er  on  wi  ng  effect# 

Pm  a- Ac  pr^eenre  -  frre  striae 

Eyne-ic  pr»*«rrre  at  horizontal  stabiliser 


Table  3. 


Definition  r>r  Herd  In  XC-3)i?A  fenatlone 


f-oTinicn  N'n.  1 
rH  1.0^ 


Reproduced  From 
Best  Available  Copy 


21 


KAVTHAr.Tc:  ■  UQ'y-Z 


Sc^si 

°T,S 


JTR 

*n 

NTR 

N 

o 

Bn 

Tn 

Nn 

Tn 

«n 

«n 

^TR 

B 

TR 

TTR 

te 

ITR 

n 

TR 

nm 

TR 


Definition 

Coefficient  of  thrust  at  ving  due  to  inflow  velocity 
Total  thrust  of  main  propellers 
Advance  ratio  main  propeller 
Advance  ratio  tail  rotor  \ 

RPM  -  main  propeller 

RTM  -  tail  rotor 

Nominal  PJ71  -  main  propeller 

Blade  pitch  angle  -  main  propeller 

Main  propeller  thrust 

Main  propeller  thrust  component  normal  to  TR 
Main  propeller  moment  (initially  turning) 

Main  propeller  moment  (initially  pitching) 

Main  propeller  torque 

Diameter  main  propeller 

Diameter  tail  rotor  , 

Blade  pitch  angle  -  tail  rotor  y 
Thrust  -  tail  rotor 

Inertia  of  rain  propeller  blades  a^id  shaft 
Angular  velocity  rain  propeller 
Angular  acceleration  main  propeller 
Inertia  of  tail  rotor  blades  and  shaft 
Angular  velocity  tail  rotor 
Angular  acceleration  tail  rotor. 


Table  }.  Definition  cf  Tyrbols  In  7T-1 Equatt one  (Ccr.'td.) 

•  ■  • 
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AXIS  SYSTEMS  FOR  X0-1U2A.  In  order  to  describe  the  dynamics  and  aero¬ 
dynamics  of  the  tilt-wing  aircraft,  one  axis  system — the  propeller  axis 
system— in  addition  to  the  inertial,  body,  stability  and  wind  axes  is 
required.  The  wind  axis  system  is  complicated  by  the  tilt  of  the  wing. 
Figure  6  shows  the  wing  tilt  angle  (iw)  and  the  variable  distances 

(x3  ^  and  ®  „  )  that  track  the  aerodynamic  center  (a.c.)  of  the  wing  as 

the  wing  is  tilted  through  the  angle  iw#  The  aerodynamic  center  of  the 

wing  as  located  in  the  x-z  body  axes  plane  is  the  origin  of  the  wing 
stability  axes.  The  wing  stability  axis  system  is  parallel  to  the  air¬ 
craft  stability  axis  system  and  has  its  origin  translated  by  x  and 

&•  c  • 

z_  „  from  the  aircraft  stability  axis  origin  (which  is  located  at  air- 
craft  c.g.  nominally)  as  in  Figure  6. 


The  additional  axis  system,  the  propeller  axis  system,  is  shown  in 
Figure  7.  The  axis  system  is  repeated  for  each  main  propeller,  so  that 
there  are  actually  four  main  propeller  axis  systems.  The  propeller  axis 
system  will  enable  the  development  of  the  propeller  thrust  (Tn)  propeller 

torque  (C^)  in  terms  of  propeller  power,  a  normal  force  (Nr)  perpendicu¬ 
lar  to  T  along  the  propeller  blade,  and  propeller  moments  (M  )  and  (T  ). 
n  n  n 

Further  development  of  these  propeller  forces  and  moments  will  be 
considered  in  the  discussion  of  aerodynamic  effects. 
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"n  la  In  tha  x-a  plana 
Tn  la  In  tha  x-y  plana 

(*)  f„  +  0,  lv  -  0 

1^  1 11  tad  from  x-a  plana  by  }  f)n, 

Tn  tlltad  fro*  x-y  plana  by  *(90-pn) 

Obtain  ooaponanta  of  and  KR 
by  Multiply lng  by  aln  fa  and  ooa  fa. 

<3)  ♦  0.  K  ♦  0 

In  addition  to  (2)  fa  rotatad  In  x-i  plana  by  *  1*. 

Mn  rotatad  fro*  x-y  plana  by  $  1*. 

Fifura  7.  Main  Fropallar  Axis  Syataai 
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At  this  point  in  the  discussion  we  will  recall  the  equations  of 
■otion  of  the  basic  aircraft  as  discussed  in  Section  III.  The  internal 
rrrvnts  (the  right  side  of  the  moment  equations)  now  include  gyroscopic 
effects  due  to  the  rotating  mass  of  an  engine.  The  gyroscopic  effects 
are  due  to  the  main  engines  and  the  tail  rotor.  For  the  main  engines 
we  have  the  following  terms  representing  gyroscopic  effects: 


For 


4  (W  cos  Sr  *  ql(IEnE)  sln  Sr 
♦  pdgflj.)  sin  iw  ♦rdgflg)  cos  iy 


-  <  W  Sin  Sr  -  VW 

the  tail  rotor  we  have 

*  <^lITRnTR 
0 


cos  iw 


fv 


for  L  term, 
a 

for  Mi  term, 
a 

for  N  term, 
a 

\>  ■ 

for  L  term, 
a 

for  M  term, 
a 

for  N  term. 

a 


In  these  terms  p,  q^  and  r  are  respectively  the  angular  velocities  around 
the  x,  y  and  *  bocjy  axes.  The  wing  tilt  angle  is  iw#  and  1,^  are 
respectively,  main  engine  and  tail  rotor  inertias  and  flg  and  are 
respectively,  engine  and  tail  rotor  rotational  speed. 


The  equations  of  motion  without  expansion  of  the  aerodynamic  terms 
and  incorporating  the  engine  gyroscopic  and  tail  rotor  terms  become: 

Xa  •  m(U  ♦  Wq^  -  Vr)  ♦  mg  sin  0 

Y  ■  m(V  ♦  Ur  -  Vp)  -  mg  cos  0  sin  ♦  , 

a  \' 

\ 

•  \ 

Z  ■  m(W  ♦  Vp  -  Uq. )  -  mg  cos  0  cos  0 

a  X 

La  -  I^p  -  ♦  pq1)  ♦  -  l22^1r  *  dgflg)  cos  Sr  “ 

sin  i„  ♦ 

Ma  "  T22^1  *  I13(l‘2  *  p2)  *  (I11  -  I33)pr  ♦  pdgflg)  sin  i^ 

♦  rdp.Rg)  cos  iw 

Ma  •  ^33^  -  *33(p  *  ♦  (Igg  -  ^ll^P^  ”  ®*n  Sr  “  ^l^E^E^ 

cos  i*  -  PItrSi  Reproduced  From 
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The  aerodynamic  force?  and  moments— X  ,  Y  ,  7.  ,  L  ,  M  and  N  will  be 

a  a  a  a  a  cl 

developed  for  the  tilt-wing  aircraft  so  that  a  complete  set  of  equations 
of  met. ion  is  presented. 


AEROrV’iAMTC  FORCES  At.'D  MOMENTS  -  XC-lli2A.  In  the  equations  of  motion 
developed  for  the  basic  aircraft  in  Section  III  a  set  of  aerodynamic 
force  and  moment  terms  were  developed  as  a  function  of  specific  dimension¬ 
less  and  aerodynamic  coefficients.  We  will  develop  aerodynamic  force  and 
moment  expressions  based  nron  the  particular  configurations  of  the  XC-1U2A 
tilt-wing  aircraffet—’^"'^ 

In  order  to  develop  expressions  for  Xa,  Y  ,  2  l  L  ,  M  and  N  ,  con- 

a  a  cl  cl  at  cl 

tributions  from  the  major  airframe  components  of  the  XC-lk2A  aircraft  will 
be  considered  separately.  The  major  components  to  be  considered  are  the 
wi  rg,  the  main  propellers,  the  vertical  tail  and  rudder,  the  horizontal 
stabilizer,  the  tail  rotor  and  the  fuselage.  After  the  aerodynamic  force 
and  moment  expressions  for  each  of  these  major  components  are  developed, 
they  will  be  summed  to  get  the  total  aerodynamic  force  and  moment  ex¬ 
pressions.  These  total  force  and  moment  expressions  can  then  be  visually 
checked  with  the  aerodynamic  force  and  moment  expressions  developed  for 
the  basic  aircraft  in  Section  III.  The  expressions  in  Section  III  will 
be  contained  in  the  aerodynamic  force  and  moment  expressions  developed 
for  the  tilt-wing  XC-11*2A. 


Before  continuing  we  will  develop  the  aircraft  body  axes  to  inertial 
axes  transformation.  The  matrices  are: 


N1  1 


N2 


N3 


COS  i|r  COS  0 

cos  >|r  sin  ©  sin  ♦ 

-  sin  ^  cos  4> 

cos  i  sin  0  cos  ♦ 
♦  sin  \|r  sin  6 

sin  cos  0 

sin  i  sin  9  sin  $ 

♦  cos  cos  4> 

sin  f  sin  0  cos  4> 
-  cos  t  sin  4» 

-  sin  9 

cos  0  sin  4> 

- 

/ 

cos  ^  cos  4* 

XI 


X2 


X3 


In  these  matrices  N1  is  north,  N2  is  east  and  N3  down  for  the  inertial 
axes,  and  XI  is  x,  X2  is  y  and  X3  is  z  for  the  aircraf|  body  axes. 

0  and  <t>  are  conventional  Euler  angles.  Their  definitions  in 
developing  the  above  matrices  are  shown  in  Appendix  A,  \ 


The  inertial  to  body  axes  rates  are: 


p 

«1 

r 


•  • 
sin  0  ♦  4 


i  cos  0  sin  i 


i|r  COS  0  COS  t 


♦  0  cos  4 
-  0  sin  ♦ 
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Avd  conversely:  . 

4*  •  p  ♦  q,  ten  0  sin  *  +  r  tan  0  cos  ♦ 

• 

0  •  q^  cos  ♦  •  r  sin  ♦ 

•  cos  ♦  sin  ♦ 

♦  -  r  - x  +Q'  - r 

COS  0  cos  0 


We  are  nor  ready  to  develop  the  aerodynamic  farces  and  nownta  for 
each  major  cosponent  of  the  aircraft.  The  first  to  be  considered  will 
be  the  affect  of  the  rain  propellers .  It  should  be  noted  that  the  ex¬ 
pressions  developed  are  for  hovering,  transition  and  conventional  flight. 


Main  Propeller*  There  are  four  mutually  similar  systems  of  axes 
used  to  describe  forces  and  norents  generated  ty  the  main  propellers 
daring  horer  and  low  aircraft  velocities.  An  analysis  which  was  developed 
by  Ling-Temoo-Vcmght  (Reference  6)  for  the  IC-1U2A  is  adopted  herein  for 
the  description  of  main  propeller  forces  and  moments.  The  subscript  n , 

(n  *  1,  2,  3,  U),  denotes  the  particular  propeller.  They  are  numbered 
left  to  right  looking  fro*  the  top— 1  and  2  are  port  propellers j  3  and  U 
are  starboard  propellers.  The  approach  used  to  develop  the  required 
propeller  equations  is  to  first  consider  the  propeller  geometry,  next 
state  the  aerodynamic  coefficients  and  finally  write  expressions  for 
the  force  and  moment  contributions  of  the  main  propellers. 

1.  Main  Propeller  Geometry.  Pro*  figure  7  the  wind  vector 

with  respect  to  each  propeller  is  formed  as  V‘  »  ♦  v*  ♦  w2. 

n  n  n  n 


For  each  propeller  wo  ruy  define  an  angle  Vhich  is  the  angle 
between  the  component  of  the  total  aircraft  velocity  (Vg)  in  the  x-s 
piano  of  the  body  ares  and  the  propeller  axis  line.  If  1^  *\0,  then 

■  v 


+o  *  (V*Li 2- 


Fro*  geometric.  con*ic!*rations  In  Figure  7  ye  pay  form  the  expression 

for  f  «hlch  is  the  angle  between  the  u  and  V  \ velocities, 
n  n  n 


cos  * 

Tn 

sin  ♦ 
Tn 


COS  flp  cos 
*  Tn>  1/2 


t 


o 


cos  B„ 

f 

CO* 


sin  i 

o 


n 


Ar  *  Pin"1 (5-) 


(h.l) 
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In  a  like  manner  expressions  are  formed  for  8^  which  is  defined 

as  the  angle  between  the  velocities  V  sin  ^  and  w  . 

n  n  n 


'in  8_  - 


V  sin  8 

n  _  _ F 

n  "  sin  ^n  "  s4n  * 


n 


W 


cos  S  * 


cos  gp,  sin  i|ro 


(U.2) 


n  _  ^ 

"  *  («“;  *  7?)  1/2  "  sln  *n 


n  n 


The  velocity  expressions  u  ,  v^  and  wn  will  now  be  developed.  The 

origins  of  each  propeller  axis  are  located  along  a  line  parallel  to  the 
y-body  axis  at  approximately  the  center  of  mass  of  each  nacelle.  In 
Figure  7  the  origins  of  the  propeller  axis  systems  are  rotated  from 
cos  8  by  angle  \|^o.  Each  propeller  axis  system  origin  is  located  by 

xn»  yn  anc?  zn  3x43  components  which  are  multiplied  by  the  appropriate 

body  axis  angular  velocity  in  order  to  give  tangential  velocity  components 

of  u  .  v  and  w  , 
v  ’  n  n 

un  "  VB  cos  6F  cos  *o  “  yn(p  sin  K  *  T  003  V  +  xn  ql  sin 

♦  znqD  cos  iw  (U.3) 

vn  "  VB  sin  ®F  (U.U) 


wn 


VR  cos  8p  sin  ♦  yn(p  cos  iy  -  r  sin  iw)  -  cos 

+  zn<l±  sin  iw  (U.3) 


In  equations  (U,3)  and  (U,5)  the  term  yn  is  a  constant.  In  fact  y^“  y^ 
and  y2  ■  y^  since  the  inboard  and  the  outboard  propellers  are  each  the 
same  distance  from  the  x-z  body  axis  plane.  The  component  xn  and  zn 
vary  as  a  function  of  wing  tilt. 


2,  Main  Propeller  Aerodynamic  Coefficients,  The  aerocfynamic 
coefficients  presented  for  the  main  propellers  follow  those  presented  by 
Reference  6,  First  let  us  define  the  advance  ratio  (Jfi)  for  each  main 

propeller  and  the  advance  ratio  normal  to  the  propeller  disk  (<T), 


60Vn 

J  ■  ..-tic  and  J'  ■  J  cos  \Jr 
n  n  D  n  n  n 


(U.6) 
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The  symbol  Nn  is  the  particular  propeller  RPM,  the  number  60  changes  RPM 

to  RPS,  D  is  the  diameter  of  the  propeller  and  Fn  the  blade  pitch  angle 

of  the  particular  propeller.  The  aeroctynamic  coefficients  are  then 
developed  in  terms  of  advance  ratio  and  blade  pitch. 


C. 


3Ct 

V  V  5F 


ac 


j»  ♦ 
n 


32S 
3J1  C 
2„ 


o 

<Jk> *  T 


n  ♦ 


3C 


o  cD  +  —E.  j»  ♦  — ( j»  )  +  — £■  b 

P  P„  3J*  n  2  n'  3B  n 

no  3  J* 


B  J» 


3P  n  3P3 J'  n  n 
2 


3  CP 

♦  - E-  B  J« 

3B3J*  n  n 


(J4. 7) 


Ch.8) 


a(C  cot  \|r) 

CN  ■  Jr  ( - TJI - 1  Vnsln  *n 

n 


a.9) 


3  (C  cot  >|r) 


3fC  cot) 


V  A  ^ — i  V„  sin  ♦„  *  i-ir— i  Vn  *„  «••») 


,  3Cm 

r  «  -i~  [ - Hi  j,  j‘ 

Mn  3J'  3V  n 


a. id 


C_  is  the  coefficient  of  thrust  (Tn),  C  is  the  coefficient  of  power 
Tn  "n 

used  to  express  torque  (C^),  CN  is  the  coefficient  of  normal  thrust 

n 

(N  ) — the  thrust  component  perpendicular  to  T  ,  C  and  are  the 
n  n  yn  n 

lateral  and  longitudinal  hub  moment  coefficients  that  appear  during  wing 
tilt. 


3.  Main  Propeller  Force  and  Moment  Expressions.  Before  ex¬ 
pressing  the  force  and  moment  contribution  in  body  axes  due  to  the 
propellers,  the  individual  propeller  forces  and  moments  developed  in 
propeller  axes  are  stated  in  terms  of  equations  (b.7),  (U.8),  (h,9), 
(h.lO)  and  (h.ll). 


T„- 

J* 

OJ 

N 

2CU 

(-£-) 

p 

CT 

CU.12) 

0 

0 

n 

Nn* 

J* 

K  2 

(r) 

CN 

a. 13) 

0 

0 

n 

Y  - 
n 

D5 

Nn  ? 

<r>2 

0 

❖ 

°*n 

u.iu 
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r  2 

M  =  rrs  (pH) 

n  u 

o 


(-2-)  0 
%  Mn 


( ii .  1 5 ) 


_  rJL) 

n  'N  ; 

o 


(—)  0 
S  j  P 

o  n 


(14.16) 


In  equations  (ii.12)  through  (U.16)  N  is  the  maximum  RPM  of  the  propellers 

and  >:  is  the  air  density  at  sea  level  on  a  standard  day.  These  equations 

then  enable  us  to  write  the  propeller  force  and  moment  contributions  in 
aircraft  body  axes.  Observe  that  each  equation  is  subscripted  by  n  so 
that  each  propeller  individually  influences  the  forces  and  moments. 


'l  a 

a'p 

V 

(T 

n*l 

h 

'  n 

a}p  " 

( -N 

n-1 

h 

n 

a]P  “ 

Z. 

(-Tn 

n-1 

w  n 


n 


w  n 


n 


w 


a. 1?) 

(u.18) 

(li.19) 


C4La) 


-  !«„>*  -  (“a'pj  h  -[(4ZaV  <4Za>p3>  *2 


-[(AY  )  ♦  ( AY  )  ]  a.  -[AY  )  ♦  (AY  )  ]  z? 

a  Pl  a  pu  i  a  p2  a  P?  2 

u  u 

-  r.  (Y  cos  8  )  sin  i  -  I,  (M  sin  8  )  sin  i 
n*L  n  n  w  n-1  n  n  \ 


(U.20) 


(AM  ) 
a'p 


M_  + 

™ PIVOT 

n-1 

T  (cos 
n 

V  ; 

2 PIVOT  ' 

f  L.  T 
n-1  n 

- 

cos 

®1 

sin 

i 

V 

cos 

sin  i  ) 

w 

-  (m2 

cos 

62 

sin 

i 

w 

*N3 

cos 

63 

sin  i  ) 

w 

>  (% 

cos 

81 

cos 

i 

w 

*NU 

cos 

eu 

cos  i  ) 

w 

♦  (n2 

cos 

fi2 

cos 

i 

V 

cos 

S3 

cos  i  ) 

w 

-  (Y  sin  6  )  ♦  E  (M  cos  6  ) 
n-1  n  n  n-I  n  n 


(14.21) 
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Whe re  \  •  1.625(T  «■  T,  )  ♦  1.092(T?  +  Tj 

PIVOT  u  c  3 


-  ((AXjn  -  (AXj^  ]  y,  -  [(AXa)„  -  (AXJ 


a  P 


p2  a  p3 


^(AYa)pi  -  (AYa)pJ  x,  ♦  [(AY^^  (AYa)p3J  x, 

It  U 

‘  nil  (Yn  cos  en}  cos  S,  *  n-l(Mn  sin  8n)  cos  ^ 


(U.22) 


In  equations  (U.17)  through  (U.22)  the  letter  p  denotes  the  effects  of 
the  main  propellers  and  p  subscripted  p^  where  n  «  1,  2,  3,  li  is  a 

particular  propeller.  For  example,  in  (AL  )p  the  term  (AZ  )pn  equals 

(-T-^  sin  i  -  N-  cos  6^  cos  i  ),  In  order  to  better  appreciate  these 

equations  let  us  consider  the  aircraft  in  normal  forward  flight  where 
the  propeller  wind  vector  is  parallel  to  the  x-z  plane  (Sn  *  0)  and 

there  is  no  tilt  of  the  wing  (i  «  0).  The  equations  (U.17)  through 

(U.22)  then  become: 

<"a>P  *  n*l  Tn  CU.23) 

(AYa)p  -  0  (U.2U) 

U 

(AZa)P  -  ^  -  Nn  (U.25) 


(ALa)p 

(AMa)p 


-  *  iy  y1  -  (-N2  ♦  N3)  y2 

^PIVOT  +  n-1  Tn  ^ PIVOT 


♦  (N1  +  ryx-L  +  (N2  ♦  N^)x2  + 


u 

I  M 
n»l  n 


(ANa)P  -  .(Tl  -  Tu)yi  -  (T,  -  T3)y2  -  ^  Yn 


(U.26) 

(U.27) 


(U.28) 


Equation  (U.23)  is  the  total  thrust  and  (U.25)  is  the  total  normal  force 
due  to  the  propellers.  Equation  (U.26)  is  the  rolling  moment  contribution 
which  will  be  zero  if  outboard  (n»l  and  U)  and  inboard  (n«2  and  3)  normal 
propeller  forces  are  balanced;  (U.27)  is  the  pitching  moment  contribution; 
and  (U.28)  is  the  turning  moment  contribution  which  will  be  negligible  if 
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the  outboard  (n«*l  and  U)  and  inboard  (n«2  and  3)  thrusts  are  balanced 
and  EYn  "0.  The  incremental  propeller  forces  and  moments  equations 

(h.17)  through  (U.22)  will  be  included  in  the  total  aerocfcynamic  forces 
and  moments* 

Wing.  The  calculation  of  wing  aerodynamics  forces  and  moments  is 
complicated  by  the  wing  tilt  during  vertical  and  transition  flight. 
These  forces  and  moments  are  developed  in  wing  stability  axes  by  first 
considering  the  wing  geometry  and  then  defining  wing  aerodynamic  coef¬ 
ficients  in  accord  with  Reference  6. 


1.  Wing  Geometry.  In  wing  stability  axes  there  occurs  an 
induced  velocity  (AV)  due  to  the  propeller  wash  across  the  wing.  This 
gives  the  effect  of  increased  lift.  In  order  to  describe  the  effect, 
a  coefficient  of  thrust  of  the  wing  (CT  s)  is  defined  as  a  function  to 
total  aircraft  velocity  (Vp).  * 


CU.29) 


where  S  is  the  total  disk  area  of  the  four  propellers  and  a  is 

r  ™ 

the  wing  dynamic  pressure. 

\  -  (q  ♦  jr-),  where  q  -  1/2  ov| 


At  low  forward  speed  during  transition  and  in  hover  the  effect  of 
CT  g  is  at  a  maximum.  The  particular  function  of  CT  g  is  dependent 

upon  wing  tilt  and  wing  flap  angle  which  should  be  determined  from 
manufacturer's  data. 


From  the  geometry  in  Figure  8  the  u^  component  of  wing  velocity 

is  u  "♦  AV .  The  velocity  u  is  the  average  of  the  u  velocities 
P  P  ^ 

developed  in  propeller  axes  and  is  the  predominant  velocity  effect  in 

wing  axes,  u  is  then  the  velocity  of  the  aircraft  (V  cos  B  )  and  as 
p  t?  i* 

such  defines  the  body  axes  velocities  J  and  W.  Any  V  in  body  axes  is 

equal  to  vw  since  we  are  concerned  with  wing  stability  axes  and  u^  is 

contained  in  the  x-z  bocty-  axes  and  x-z  wing  axes  planes.  The  total 
vcloc  ity  in  wing  axes  (V  )  is  as  follows: 

V  ■  (w2  +  lu  ♦  AV)2  +  v2)  ^2  ,  where  v  ■  V  (U.30) 

w  p  p  w  '  w  7 


From  Figure  8,  w  -  w  -  «  !J  sin  i  ♦  W  cos  i 

p  w  w  w 


and  u  •  U  cos  i  ♦  W  sin  i 
p  w  w 


a. 31) 
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From  the  foregoing  the  wing  angle  of  attack  (o ), 
angle  (8^)  can  be  stated  as: 


tan”1  ( 


w 


P  i 
TW>  * 


and  the  wing  sideslip 


0w  32) 


with  the  sign  of  aw  as  positive  down  from  plus  x^. 


8^  -  tan~x  — - -  .  — . —"■■■■  (U.33 ) 

'/(u  *  AV)2  +  Wp2 

The  expression  for  induced  velocity  (AV)  is  derived  from  momentum  theory.1 
AV  is  similar  to  the  mean  inflow  velocity  (VT  )  developed  for  helicopter 

inflow  analysis. ^  mean 


...  f  2  2T  x  1/2 

av  .  -  Up  +  (up  ♦  -jg-) 

Rearranging  we  have: 


(u  ♦  AV)  -  (u2  ♦  1/2 

P  P  PSp' 

but  Vw  •  up  ♦  AV 

and  u  •  V  ,  so  that 
p  B' 

substituting  for  (u  +  AV)  and  u  ,  and  multiplying  by  o  we  have: 

P  P  T 


(U.3U) 


CU.35) 


1/2  ovJ-1/2  oV\  .  |_ 

P 

2 

q  •  1/2  oV ,  where  q  is  the  dynamic  pressure 
but  q|f  ■  1/2  pV2 
so  that  we  have: 

^  ■  (Q  +  |-)  U»36) 


1.  For  example,  see  Airplane  Aerodynamics.  Reference 

2,  NAVTRACEVCEN  1205-1,  Section  3. 


NAVTRArEVCEN  1205-2 


Tn  this  manner  ^  develops  naturally  from  the  momentum  theory  definition 

of  AV.  The  wing  cfynamlc  pressure  is  the  pressure  used  in  equation  0j,29' 
for  C.p  Consequently,  we  have  now  related  the  coefficient  ^  with 

the  induced  velocity,  AV,  through  the  wing  dynamic  pressure,  q^.  We  are 

now  ready  to  develop  the  wing  aerodynamic  coefficients. 


2,  Wing  Aerodynamic  Coefficients,  In  order  to  develop  the 
wing  forces  and  moments  five  aerodynamic  coefficients  will  be  defined 
as  in  Reference  6.  Since  the  development  is  in  wing  axes,  rolling  (p) 
and  turning  (r)  rates  necessary  to  define  these  coefficients  are  trans¬ 
formed  from  body  axes.  The  pitching  rate  (q^  is  the  same  in  wing  axes 

since  q^  lies  in  the  x-z  plane  and  the  plane.  The  transformation 

of  the  angular  rates  (p  and  r)  is  accomplished  by  a  rotation,  t .  From 
Figure  fl  i  is  equal  to  (i  -  aw).  ^e  can  t,^en  write  for  the  wing  rolling 

rate  (pw)  and  the  wing  turning  rate  (rw)  the  following  equations. 


■  p  cos  £  -  r  sin  £ 
rw  *  p  sin  £  ♦  r  cos  £ 


a. 37) 
(I. 38) 


The  aerodynamic  coefficients  for  the  wing  are  C^,  Cj ,  (C  ),  (C  )  and 

(Cn  ).  They  are  defined  as  follows  in  accordance  with  their  development 
w 

in  Reference  6, 


2  dCD  *2°T)  2 

V  .  6„  ♦  — .  6  r 


cd  "  V  +  kcl  +  VST  *  6f  + 


(U.39) 


2 

The  strong  flap  dependence  (6  F)  in  the  expression  above,  is  due  to 
the  inportance  of  flap  during  transition, 

0\Jw 


CL  "  CL  *  CL  *  6F  +  Ct  •  aw  ♦  'TTJjT 
L  Lo  l6F  LaF  W  36F 


,  6F  ,  a 


(U.UO) 


Ccft)  - 


(Cn  ) 

\ 


3C 


rn 


36F 


m  .  6F  ♦  m£—  C 


KT  X  *  Ql 


(U.iil) 


(U.U2) 


(V 


C  ,  a.  ♦  c  ,  6A  ♦  T-—  C  .  p  ♦  -2-  c  ,  r 

\  K  "6*  "p  P»  \  w 


(hM) 
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3.  Wing  Force  and  Moment  Expressions.  Before  writing  the  force 
and  moment  expressions  for  the  wing,  equations  (1^.39)  through  (iI,U3)  will 
be  transformed  to  the  aircraft  body  axes  through  the  angle  &  . 

The  body  axes  wing  coefficients  are  then: 


(Cx)  -  -  CD  cos  £  -  CL  sin  £  (ii.Ui) 

CC2)  -  CD  sin  £  -  C  cos  £  CU.U5) 

w 

(C  )  -  (C  )  cos  £  +  (C  )  sin  £  (k,U6) 

^  w  w  w 

(cm>  -  K  )  -  *-f  (c*>  *  x-f  CU.U7) 

w  W  w  w 

(C„)  m  -  (r‘i  )  sin  £  ♦  (C  )  cos  £  CU .U8) 

w  w  w 


x  and  2  are  the  respective  distances  from  the  c.g.  of  the  aircraft 
SC  sc 

in  body  axes  to  the  aerodynamic  center  (ac)  of  the  wing  in  the  x-z  plane, 
c  is  the  mean  aerodynamic  chord. 

The  wing  force  and  moment  contributions  to  the  total  force  and 
moment  equations,  (U.UU)  through  (U.U8)  can  be  expressed  as  forces 

and  moments.  For  example,  if  we  have  a  coefficient  C  ,  a  force,  X  . 

x  'a* 

is  immediately  defined  as  Sq.  *  S  is  the  wing  area  of  the  aircraft 

and  q  the  d5rnamic  pressure.  l..e  coefficients  in  equations  (h.UO 
through  (U.U8)  are  immediately  expressible  in  terms  of  their  respective 
force  and  moment  contributions. 


■ 

(c  ) 

x 

s  tf(o  ) 

(U.U9) 

w 

w 

9 

(«,)  • 

“V 

S[f(CT(S)  %] 

(U.50) 

w 

w 

(AL*)  - 

bS  [f(C-  .)  o  ] 

1  * 

(U.51) 

V 

w 

C»Ma)  - 

cstf(0T(S)  q*] 

(U .52) 

w 

w 

(4N.)  " 

W 

u 

bS(CT(S) 

(U.53 

w  w 


T.  Use  of  coefficients  to  express  forces  and  moments  can  be  reviewed 
in  References  2  and  7. 
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Equations  (U.U9)  through  (U.53)  are  the  wing  force  and  moment  contri¬ 
butions.  [f(CT  s)  qw]  is  the  dynamic  pressure  term  incorporating  the 

effects  of  having  thrust  produced  by  the  propellers  and  having  wash  across 
the  wing.  Consider  what  happens  if  the  aircraft  is  flying  and  the  engines 
are  turned  off  (T  «  0).  Then  -  q  from  equation  (U.36)  and  f("T  g)  £5  1 

from  data  available  in  Reference  6.  Thus  [ f ( q  with  the  engines 

off  and  equations  (U.U9)  through  (U.53)  are  wing  force  and  moment  equa¬ 
tions  that  would  be  expected  to  occur  in  a  jet  aircraft. 


Vertical  Stabilizer  and  Rudder.  Forces  and  moments  for  vertical 
tail  (vt)  and  rudder  arise  from  the  relative  wind  pushing  against  the 
vertical  tail  surfaces  thereby  causing  a  turning  moment  due  to  control 
input  to  the  rudder.  This  gives  side  force,  as  well  as  rolling  and 
turning  moments. 

Wind  pushing  against  the  vertical  tail  and  rudder  yields  a  side 
force  (Ya)vt  and  rolling  moment  (ALa)vt  respectively  which  are  non- 

dimensionalized  in  terms  of  aerocjynamic  coefficients  as  C^.  and  C^.  The 

rolling  moment  is  coupled  in  rolling  velocity  p  and  turning  velocity  r. 
We  can  than  write  for  C  and  CL  the  expressions,  (U.5U)  and  (U.55). 


Br.  ♦  C  .  6R 
h  y6R 


(U.5U) 


c  is  the  change  in  side  force  coefficient  with  changing  sideslip  angle. 
y  0 

It  acts  as  a  damping  term.  C  is  the  change  in  side  force  coefficient 

y  5^ 

with  rudder  deflection  and  represents  the  controllable  term  in  the  ex¬ 
pression.  This  is  important  in  the  use  of  automatic  stabilization 
systems  for  aircraft. 


C, 


'0p 


^  ♦  C  .  6R  ♦  g— - 
f  *6R  *v° 


[C  .  p  ♦  C  .  r 1 
B  *p  r 


CU.55) 


C.  is  the  change  in  rolling  moment  with  variation  in  sideslip  angle. 
*6 

It  is  caused  fcy  the  position  of  the  vertical  tail  normally  above  the 
X  axis  and  other  factors  as  wing  sweepback  and  dihedral.  C,  is  the 

6R 

change  in  rolling  moment  coefficient  due  to  rudder  deflection.  is 

P 

the  roll  damping  derivative.  Ct  is  the  change  in  rolling  moment  co- 

r 

efficient  with  change  in  yawing  velocity. 


Rudder  deflection  will  give  a  turning  moment  (AN'a)vt  which  can  be 
expressed  by  the  aeroctynamic  coefficient  C^,  and  the  expression  for 
is: 


37 


n 


6F 


»F  *  Cn 


KR 
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.  ,R  rcn  .  P  +  Cn  .  r] 

p  p  r 


(a. 56) 

is  rudder 


0  is  the  weathercock  or  static  directional  derivative,  C 

r’li  nf,R 

effectiveness,  the  change  in  yawing  moment  coefficient  with  rudder 
deflection.  0  is  the  chanpe  in  yawing  moment  coefficient  with  varying 
np 


roliinc  velocity.  C  is  the  yaw  damping  derivative.  The  tail  is  the 

r 

main  contribution  to  Cn  . 

r 


The  forces  and  moments  for  the  vertical  tail  can  then  be  expressed 
in  the  following  equations: 


f*Ya!vt  1 

1  S  (— ) 
y  q  q 

a. 5?) 

(ALa)vt 

5i 

a. 58) 

(ANa)vt 

Cn  bSq(^t) 

(a. 59) 

Fere  q  is  the  free  stream  dynamic  pressure  and  is  the  vertical  tail 

cfynamic  pressure.  This  force  and  these  moments  will  included  in  the 
total  aerodynamic  forces  and  moments. 

Horizontal  Stabilizer.  The  equations  for  forces  and  moments  for 
the  horizontal  stabilizer  (hs)  will  be  developed  in  the  standard  manner. 

We  define*  ci^  *  dp  ♦  i^  -  c 

Here  is  the  angle  of  attack  of  the  tail,  Op  is  the  angle  of  attack  of 
the  fuselage,  i  is  the  angle  of  incidence  of  the  tail  and  e  is  the  down 
wash  anrle. 

In  accord  with  Reference  6  we  have 

c  -  [1*  -  .01*21  (\,  *  aF)l  [(CL)  f(CT  s)]  ♦  iw  -  aw  ♦  ap 

w  * 

for  the  yC-lU2A. 

The  lift  (CT  )  and  (Cr  )  coefficients  of  the  tail  can  then  be  ex- 
L't  tt 

pressed  in  the  following  relations 
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a 


t. 


+ 


K  >'•.  ) 

t  i._ 


The  horizontal  stabilizer  (hs)  can  contribute  forces  in  the  x  and  z 
directions,  and  a  pitching  moment.  The  equations  are  as  follows: 

fAXa)h^  *  -  [Op  ccs  (ifc  -  at)  +  CL  sin  (it  -  )  ]  Sq( - )  (U.6G) 


'hs, 


(4Vh=  *  -  ;-cr.  sin  (it  *  at)  *  CL.  cos  (it  -  V1  (U,61) 


(AM  ).  =  -  (AX’).  .  h.  +  (AZ„ ).  .  1, 

a'hs  a'hs  hs  a'hs  hs 


c2S  0 


U~  ‘Cm  •  Vs  *  C„.  •  W’ 

r\  n 


a. 62) 


lhs  is  distance  from  aircraft  c.g.  to  the  aerodynamic  center  (a.c.) 

of  the  horizontal  stabilizer  and  h^g  is  the  height  of  a.c.  above  the  c.g. 

Both  lhg  and  h^g  are  measured  in  the  x-8  plane  of  the  aircraft  bocty  axes. 

S  is  the  wing  area,  p  is  the  air  density,  c  is  the  mean  aerodynamic  chord 
and  the  angle  (i^  -  a^)  is  used  to  transform  and  to  body  axes. 

t  t 


C  is  the  pitch  ds 
m 

«1 


ing  derivative —the  change  in  pitching 


moment  coefficient  as  pitch  velocity  is  changed.  is  due  to  the 

a 

fact  that  when  the  wing  undergoes  a  change  in  angle  of  attack,  there  is 
a  time  lag  before  the  change  in  downwash  is  fslt  at  the  tail  surfaces, 

(AX  ),  ,  (AZ  ).  and  (AM  )  will  be  included  in  the  total  aero- 
dynamic  forces  and  moments. 


Tail  Rotor.  The  z  force  (AZa)TR  and  the  pitching  moment  (AMa)TR 

developed  at  the  tail  will  be  obtained  by  finding  a  tail  rotor  (TK) 
advance  raoio  (JTp).  From  J^R  the  tail  rot  r  thrust  coefficient  (CT 


) 


and  tail  rotor  power  coefficient 


(CpJ 

fTR 


will  be  found. 


In  turn,  the 


tail  rotor  thrust  (T^)  and  torque  is  obtained  and  consequently 

( A"a^TR  and  ^AMa^TR*  TTR  Posltive  the  “z  direction. 
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Wp  define  the  total  tail  rotor  velocity  (VTR)  as  VTp  -  [(upp)2  +  (wTp)2]]/ 
Here  upp  and  w^p  are  defined  as: 

UTR  "  VR  Sin  Uo}TR 

WTR  Vn  cos  <VTW  “  1TRql 

I K 

lmR  is  the  distance  from  the  center  of  the  tail  rotor  hub  to  the  aircraft 

o .  p ,  and  (\|^  )  locates  the  tail  rotor  with  respect  to  the  aircraft  booty 
TR 

axes. 

c*0)  ■  90°  ♦  (aF  -  e) 

TR 

Here  e  is  again  the  downwash  angle  and  is  defined  as  for  the  horizontal 
stabilizer. 


The  advance  ratio  for  the  tail  rotor  is 


60  VTR 

NTR  ^TR  ^ 


6°C-Wtr) 


where 


is 


the  RPM  and  D_„  the  diameter  of  the  tail  rotor. 

rK 


We  now  define  (C  '  )  and  (C  )  as  is  done  in  Reference  6. 

TR  PTR 


^rp  “  (  B_  ) 

TR  lTR 


^T 

-h*  (j.  ) 

3.%.  WTR; 


32C  2 

C  -  (blr)z 

PTR  3P^ 

Here  BpR  is  the  collective  pitch  of  the  tail  rotor  blades  and  C^,  (Rpp) 

TR 

is  developed  from  a  curve  of  C^,  plotted  versus 

TR 


The  thrust  (T^)  and  torque  (Q^p)  of  the  tail  rotor  is  then  written 


T  ■  1 

TR 

j.  N„  2 

nijL  (— £— )  (  nr*| — )  C_ 

Po  o  TR  TTR 

(U.63) 

^TR  "  J 

*. 1  ’  - 

(U.6M 

:.0 
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Consequently  the  force  and  moment  terms  can  be  written  directly. 

<*VTR  ■  *TH 


a.  65) 


"  TTRlTR 
TR 


(U.66) 


(ANJ  - 
TR 


(U.67) 


Fuselage.  In  a  very  direct  manner  we  can  write  the  effects  of  the 
fuselage  (F)  on  the  total  aerodynamic  forces  and  moments. 

We  have  for  the  forces 


0VtSCD 
F  o 


CU.68) 


is  the  equilibrium  drag  coefficient, 
o 

<«.>r  ■  *  f  pVB  S  °y  •  SF 


(U.69) 


C  is  the  change  in  side  force  with  respect  to  a  changing  sideslip  angle. 

y<v 


a*2.’  ■ -?  °'4sol  .  «f 

r  On 


CU.70) 


CT  is  the  change  in  lift  coefficient  with  varying  angle  of  attack. 
aF 

This  is  also  known  as  the  lift  curve  slope. 

We  have  for  the  moments: 


(AL  )  -  0 

F 

(4Ma)  •  ?  0  S  c  0B  ♦  \  0  v|  Sc  Ca  .  «F 

F  o  Op 

C  is  the  aerodynamic  pitching  moment  coefficient  in  equilibrium  flight 
o 

and  Cm  is  the  longitudinal  static  stability  derivative. 

°F 


(4Na)  -  \  0*1  S  b  C  .  SF 


U.71) 


hi 
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is  the  static  directional  or  "weathercock"  derivative* 


Tn  the  above  expressions  q  ■  j  o  VR  which  is  the  cfynamic  pressure* 

b  is  the  win£  span,  c  is  the  mean  aerodynamic  chord  and  S  is  the  wing 
area. 


For  hovering  and  in  transition  regions  Bp  is  assumed  small  so  that 

COS  g  *  1. 


EQUATIONS  OF  MOTION  -  XC-1U2A.  The  aerocfomamic  force  and  moment  terms 
developed  for  each  of  the  major  aircraft  components  will  be  combined 
and  finally  expressed  in  the  equations  of  motion  for  the  total  aircraft. 


To1 al  Aerodynamic  Forces  and  Moments,  The  total  aerodynamic  forces 
and  moments  are  as  follows: 


X,  =  (AX.)  *  (AXj„  *  UXa)hg  ♦  (AXJ 


aw 


a  F 


Ya  "  (AYaV  +  (Vvt  *  <AYa>F 

Za  “  (AZa*p  *  (AZa}w  *  (AVha  *  (ftZa)TR  ♦  (AZft)p 

La  -  (ALa)p  ♦  (ALa)w  ♦  (ALa)yt 


Ma  -  (AMJ^  ♦  (AMJ„  -  (AMa)hs  ♦  (AM.)to  +  (AM.) 


a  P 


a'w 


a  TR 


a  F 


N.  -  fANj^  ♦  (ANJ„  ♦  (ANa)yt  ♦  (ANJtb  *  (ANJ 


a  p 


aw 


a  TR 


a'F 


Equations  of  Motion  Expanded.  The  forces  and  moments  are  pre¬ 
sented  in  body  axes.  These  equations  are  subject  to  simplification  in 
representation  of  aerodynamic  propeller  forces  and  moments,  when  flight 
test  data  are  available  from  which  to  etc  ■’mine  typical  magnitudes  of 
terms. 


1*  X  Force  Equation 


rxB,(CT,S)S.  ’  *(Tncos  \  '  Nn  oos  6n  sin 
-ffljj  cos  (it  -  at)  ♦  CL  sin  (it  -  a^)  ]  Sq(-ip)  -  Ea  C_, 


m(H  ♦  Wq1  -  "r  mg  sin  0 


NAVTRA.IEVCEN  1205-2 


?,  Y  Force  Equation. 

l  (-N  sin  B  )  +  C£  bSq  (— -1  ♦  Sq  C  .  ep 

n*l  bf 

■  m(V  +  Ur  -  Wp)  -  mg  cos  0  sin  $ 

3,  Z  Force  Equation, 

h 

Cz)5f^T  s)qw  +  l  (-T  sin  iw  -  Nr  cos  Bn  cos  iw> 
w  *  n*l 

+[CD  sin(i.  -  a.)  *  C.  cos  (it  -  at)  3  Sq(— —■)  ♦  TTR  -  SqU.  .  ap 
t  t  q  ap 

-  m(W  ♦  Vp  -  Uq^)  -  mg  cos  0  cos  ♦ 
k.  Roll  Equation. 

(Cl>bSf'<CT,S>''w 

♦[(♦T^  sin  iw  +  N1  cos  ^  cos  i^)  -  (+T^  sin  cos 

cos  iw)]  y1 

+[(+T2  sin  \  +  N2  cos  ®2  008  “  ^+^*3  s*n  *w  *  N3  cos  83 

cos  VJ  *2 


+  (+1^  sin  8^  ♦  sin  8^)5-,  ♦  (+Ng  sin  sin  B^)  z ^ 

h  U  %rt 

-  I  (Yn  cos  Bn)  sin  iw  -  I  (Mn  sin  Bn)  sin  i  ♦  C£  b  Sfcj(~~~) 

n*l  n-i  q 

■  InP  ♦  I13(r  ♦  pqx)  ♦  (I33  -  I?2)  q3r  +  (Ijf^.)  cos 

“  ^l^E^E^  8in  Sr  *  ql  ITrf,TR 
5.  Pitch  Equation. 


(Cm)cSf(CT  s)q|# 
w  * 

\ 


♦  n-  ♦ 

PIVOT 


h 

z 

l 


Tn(co»  iw)  «plT0t  ♦ 


U 

I  Tn  (sin  ij  xpiT0t 


U3 
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-  (N^  cos  8^  sin  ±v  *  cos  8^  sin  iw)  a^ 

-  (Mg  cos  8g  3in  iw  cos  8^  sin  i^)  a>2 

+  (N-^  cos  61  cos  cos  8^  cos  i^) 

h 

♦  (N2  cos  &2  000  **  ♦  ^3  cos  63  cos  x2  "  I  (Yn  sin 


4  I,  (Mn 

COS 

6n} 

n*l 

‘  tcDt 

cos  (i^  - 

at5 

♦  CL 
Lt 

sin 

^t  " 

*  tc»t 

sin  (i^  - 

a0 

1 

0 

etT 

cos 

(1t  - 

c2Sp 

*  m r 

-  [Cm  •  «1VB  4 

Gm*  ‘ 

n 

r  W] 

4  ttr' 

,\i 


*1 

“  ^22^  *  Ii3^r  ~  P  )  +  “  I33^pr 

♦  p(lEV  sin  *V  *  r(IEnE)  C0S  *V 

6.  Taw  Equation. 

<Pn)bSf(CT  ,sK 

w  • 

-  [ (T^coa  iw  -  Nj  cos  sin  i^)  -  (T^  cos  i^  -  cos  B^sin  Srlyl 

-  [(T2  cos  i^-  H2  cos  Bg  sin  5^)  -  (T^  cos  cos  S^sin  iw]y2 

-  (+N^  sin  8^  ♦  sin  B^)x^  “  sin  8g  ♦  sin  ^3)^ 

u  u 

-  r  (I  cos  8  )  cos  i  -  I  CM  sin  8„)  cos  i 

n  n  w  n  ii  w 

n»l  n-1 

4  Cn  bSq(^  4  q  ^  Cn.  •  SF  4  ®TR 


*  hf  4  X13(P  ‘  V5  4  (I22  -  I13.)pql 
-  (h°l}  ,ln  Sr  -  VW  co*  Sr  *  PITR‘,TR 


No  simplifications  hare  bean  nada  on  the  right  side  of  these  six 
equations,  fccteraal  stores,  rough  air,  or  landing  gear  conditions  haae 
not  been  developed  in  these  aquations. 
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TILT-DTCT 

The  VZ-hDA.  VTOL  aircraft  of  DOAK  Aircraft  will  be  used  in  the  ex¬ 
ample  of  simulation  equations  of  motion  for  a  tilt-duct  aircraft.  Much 
of  the  development  which  follows  has  been  based  upon  information  gleaned 
from  References  1  and  8.  Soule  and  Baumgarten  in  Reference  8  present  a 
method  for  calculating  the  aerodynamic  forces  and  moments  for  the  VZ-UEA 
and  Rreul  in  Reference  1  presents  a  study  of  handling  qualities  for  a 
tilt-duct  aircraft.  From  this  information  and  the  presentation  set 
forth  for  the  tilt-wing  example,  a  well-defined  set  of  engineering  simula¬ 
tion  equations  for  the  tilt-duct  will  be  presented  realizing  that  an 
actual  simulation  will  be  dependent  upon  the  type  and  usefulness  of  the 
aerodynamic  data  available  from  the  manufacturer  of  the  aircraft. 

In  order  that  we  may  have  some  idea  of  the  specific  details  of  the 
VA-UIA,  Figure  2  is  a  three-view  drawing  of  the  VZ-hEA,  Table  U  is  a 

selected  list  of  physical  characteristics  of  the  aircraft  and  Table  5 
is  a  list  of  symbols  used  ±n  developing  the  equations  of  motion.  As 
much  as  possible,  repetition  of  symbols  will  be  used  in  the  various 
examples  of  V/STOL  aircraft  in  this  section  of  the  report. 

From  this  general  idea  of  the  VZ-liEA,  we  may  now  develop  the  mathe¬ 
matical  model  for  the  VZ-UDA.  First,  define,  if  necessary,  applicable 
axis  systems.  Second,  develop  the  aerodynamic  coefficients.  Finally, 
write  the  equations  of  motion  for  the  VZ-iiDA.  As  needed  we  will  borrow 
from  the  results  obtained  in  the  tilt-wing  equations  of  motion  and  from 
the  results  obtained  in  References  1  and  8. 

AXIS  SYSTEMS  FOR  THE  VZ-UDA.  In  order  to  describe  the  dynamics  and  aero¬ 
dynamics  of  the  tilt-duct  aircraft,  one  axis  system,  the  duct-axis  system, 
in  addition  to  the  inertial,  body,  stability  and  wind  axes  is  required. 

The  wind  axis  system  is  conqplicated  by  the  tilt  of  the  ducts.  Both  of 
the  ducts  tilt  through  the  same  angle  simultaneously.  Figure  9  shows  the 
wing  tilt  angle  (!„)  and  the  duct  incidence  angle  (ip).  Duct  incidence 

is  described  in  the  aircraft  stability  axes.  The  force  and  moment  coef¬ 
ficients  for  the  ducted  fans  are  described  in  the  duct  axis  system.  There 
are  two  duct-axis  systems— one  located  at  each  duct.  This  is  similar  to 
the  repeated  propeller  axis  systems  developed  for  the  XC-LU2A.  We  define 
the  duct  axis  system  by  the  following  three  mutually  orthogonal  velocity 
vectors.  First,  there  is  a  side  duct  velocity  (Vg)  parallel  to  the  y 

body  axis  and  in  a  minus  y  direction.  Second,  there  is  an  axial  velocity 

(V  )  into  the  center  of  the  duct.  Finally,  there  is  a  normal  velocity 
A 

(V^)  parallel  to  the  x-z  plane  of  the  aircraft  body  axis  system  and 

perpendicular  to  Vg  and  V^.  If  there  is  no  duct  incidence,  the  duct-axis 

systems  as  shown  in  Figure  10  are  similar  to  the  x,  y,  z  body  axes.  The 
duct-axes  (xm#  ym,  and  Xjg,  yK,  *tt)  are  right-handed  systems  with 

y  parallel  to  y^  and  y^g  and  x  in  the  same  direction  as  and  x ^  and 

z  in  tiie  same  direction  as  and  z^g  when  1^  is  zero. 
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Airplane: 

—  ..A  —  — — 


Over-all  length 

31.2  ft. 

Over-all  span  (including  ducts) 

25.5  ft. 

Over-all  height 

Propulsive  Device: 

10.3  ft. 

Fan 

“TVpe 

fixed  pitch 

Number  of  blades 

8 

Diameter 

a  ft. 

Blade  root  chord 

0.312  ft. 

Blade  tin  chord 

0.18U  ft. 

Solidity 

0.117 

Blade  angle  ( . 7  5R ) 

0.306  radians 

Hub  to  tip  diameter  ratio 

0.333 

Fan  location,  position  of  hub  from  duct  lip 

0.806  ft. 

Duct 

Centerbody 

Diameter  f maximum)  1.333  ft. 

Over-all  length  5.79  ft. 


Trim  Vane 
Area 

Span 

Airfoil  section 

Flap  chord,  25  percent 

Movement 


5.66  ft.2 

U.525  ft. 

NACA  0009 
0.313  ft 

0  to  ♦  0.3U9  radians 


Inlet  Guide  Vanes 
Number  of  vanes 
Area  (per  vane) 

Span 

Chord  (constant) 

Airfoil  section 

Vane  location,  position  of  C/U  from  duct  lip 
Movement 


1 h 

0.333  ft. 

1.333  ft. 

0.25  ft. 

NACA  652  -  015 
O.li  ft. 

♦  .17h5  radians 


Wing  Group: 


Wing 

Area  (excluding  ducts) 
Span  (excluding  ducts) 
Mean  aerodynamic  chord 
Aspect  ratio 
Airfoil  section 
Dihedral 


96.0  ft.2 
16.0  ft. 

6. 08  ft. 

2.67 

NACA  2U18  (mod.) 
0  radians 


A  r craft  Characteristics 


.!/■ 


Table  h 


Mr 
V  i 
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Aileron 

Area  (per  aileron) 
Span 

Movement 


6.0  ft.2 
3.0  ft. 

♦  a2b2  radians 


Empenage  Group: 


Horizontal  Tail 
Area  (projected) 

Span  (projected) 

Mean  aerodynamic  chord 
Aspect  ratio 

Airfoil  section 

Incidence 

Dihedral 

Elevator  movement 

(stabilizer  incidence,  0  radians) 

(stabilizer  incidence,  +1.7U5  radians) 


28.5  ft.2 
11.8  ft. 

2.5  ft. 

U.9 

NACA  0012  (Mod.) 
variable,  0  to  +  0,17 U5 
0.17ii5  radians 

♦0.u72  radians  (T.E.  down) 
— 0.2i7 2  radians  (T.E.  up) 
♦0.297  radians  (T.E.  down) 
-0.576  radians  (T.E.  up) 


Vertical  Tail 


Area 

13.9  ft.2 

Span 

5.2  ft. 

Mean  aerodynamic  chord 

2.8  ft. 

Aspect  ratio 

1.95  (geometric); 

Airfoil  section 

3.02  (effective) 
MAC A  0012  (Mod.) 

Fuselage: 

Length  (excluding  reaction 

control  vanes) 

29.3  ft. 

Maximum  height  (approx.) 

U.12  ft. 

Maximum  width 

3.0  ft. 

Reaction  Control  Vanes 

Area  (within  15  inch  tail 

pipe  diameter) 

1.012  ft.2 

Span  (within  15  inch  tail 

pipe  diameter) 

1.25  ft. 

Aspect  ratio  (effective) 

1.32 

Airfoil  section 

NACA  0009 

Movement 

♦  0.612  radians 

Table  li. 


Selected  Aircraft  Characteristics  -  VZ-UDA  (Cont’d.) 
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SYMBOL 


DEFINITION 


A 

c„ 


C, 


r* 

'"S 


yM 


m 


c 

D 

g 

iD 

*W 

i-HT 
1 VT 

<h 

m 

N 

r 

p.. 


Aspect  ratio 

Thrust  coefficient  (for  velocity  *  v^) 

Thrust  coefficient  (for  velocity  ■  V^) 

Normal  force  coefficient  ■  normal  force/on^E^ 

Side  force  coefficient  ■  side  force/on^D^ 

Pitching  moment  coefficient  =  Pitching  Moment/ pn^P'’ 

Total  rolling  moment  coefficient  for  ducted-fan 

Total  pitching  moment  coefficient  for  ducted-fan 

Total  yawing  moment  coefficient  for  ducted-fan 

Lift  Coefficient 

Drag  Coefficient 

Duct  chord 

Wing  chord 

Fan  diameter 

Gravitational  acceleration 
Dieted  fan  incidence  angle 
Wing  incidence  angle 
Horizontal  tail  incidence  angle 
Vertical  tail  in^uif  e  angle 

Axial  advance  ratio  (V^/E) 

Mass  flow  rate  (slups/sec.) 

Fan  speed  (R.P.S.) 

Roll  rate  about  X-Axis  (Rad. /Sec.) 

Roll  rate  about  v  (Rad, /Sec.) 
w 

Diet,  man 


ar  . 


T r> p  I  •  -tion  of  Symbols  -  V/.-.'.I) A 
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5YMP0L 

h 


*1 

r 

S 


4-  C 

P 


% 

0 

D 


DEFINITION 

Wing  spam 

Pitch  rate  about  Y-Axis  (Rad. /Sec.) 

Yaw  rate  about  Z-Axis  (Rad. /Sec.) 

Area  -  Wing 

Residual  Thrust 

Axial  velocity  into  duct 

Side  velocity  at  duct  entrance  parallel  to  X-Z  plane 
Side  velocity  at  duct  entrance  parallel  to  Y-Axis 
Jet  exhaust  velocity  (see  Appendix  G) 

Moment  arm  from  center  of  gravity  along  X-Axis 
Moment  arm  from  center  of  gravity  along  Y-Axis 
Moment  arm  from  center  of  gravity  along  Z-Axis 
Angle  of  attacK 
Arc  tan  (vr_Aa_) 

Sideslip  angle 
Arc  tan  (VS_/VA_) 

Angular  displacement 
Downwash  (♦)  or  upwash  (-) 

Mass  density  of  air 
Mass  density  of  jet  efflux 
Sidewash  angle 
Duct  (Subscript) 

Table  5.  Definition  of  Symbols  (Cont’d.) 
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SYMHOL 


E 

F 

HT 

I 

IV 

J 


M.G.C. 


N 

OL 

R 

RC 

T 

TV 

VH 

VT 

W 

W 


EEFIN1TI0N 

Elevator 

Fan,  fuel  or  fuselage 
Horizontal  tail 
Interference 
Inlet  vanes 

renotes  port  (j  ■  1)  or  starboard  (j  ■  2) 

Mean  geometric  chord 

Centerbody  nose,  nose  or  normal 

Zero  lift 

Radial  or  rudder 

Reaction  control 

Test  or  thrust 

Trim  vane 

Horizontal  vane 

Vertical  tail 

Vertical  vane 

Wing 


Table  5.  Definition  of  Symbols  (Cont'd.) 
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S>  "  *»  *  6D 

i  “  wing  incidence  angle 
w 

6^  *  change  in  duct  angle 

i  -s  duct  incidence  angle 
D  ' 

Figu-e  9.  Duct  Tilt  Axis 
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VA  Axial  velocity  into  duct 

Vu  parallel  to  x-z  plane 
rt 

V(,  parallel  to  y-axis 


\ 

V, 


f  -1,\ 
aD  “  tan  (r) 

A 

-1  VS 
^  -  tan  (~) 

A 

x,  y,  z  are  body  axes. 


Y 

z 


Figure  10.  Duct  Fan  Axis 
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Before  developing  the  aerodynamic  forces  and  moments  it  should 
be  noted  that  the  VA-UDA  has  gyroscopic  effects  due  to  the  rotating 
masses  In  the  duct  fan.  In  the  following  expressions  i  is  the  fixed 

wing  incidence  angle  and  6^  is  the  variable  duct  angle.  For  the  ducted 

fans  we  have: 


*  (IEV  cos<1v  *  V  -  ‘>l(IEnE) 

*  Pfljllj.)  sin(iK  *  8d)  -  r(IEnE) 

“  sln(iw  *  ^T?  ”  P1^E^E^ 


sin  (iM  ♦  f»D) 
nos  (iw  ♦  »„) 
cos  (iy  ♦  6jj) 


f or  L  term, 
a 

for  M  term, 
a 

for  N  term, 
a 


These  expressions  will  be  added  to  the  internal  moments  in  the  final 
equations  for  the  VZ-UDA,  Next  the  forces  and  moments—  X  ,  Y  ,  Z 

L  ,  M  ,  and  N  will  be  developed  for  the  tilt-duct  so  that  a  complete 

cl  a  ol 

set  of  equations  of  motion  can  be  written. 


AERODYNAMIC  FORCES  AND  MOMENTS  VZ-liDA.  In  order  to  develop  expressions 
for  the  aerodynamic  external  forces  and  moments,  contributions  from 
the  major  airframe  components  will  be  considered  separately.  The  major 
components  to  be  considered  are  the  wing,  the  ducted  fans,  the  vertical 
stabilizer  and  rudder,  the  horizontal  stabilizer  and  elevator,  the 
fuselage  and  the  reaction  controls.  As  with  the  XC-lii2A  after  the  aero¬ 
dynamic  force  and  moment  expressions  are  developed  for  each  of  these 
major  components,  thoy  will  be  summed  to  get  the  total  aerodynamic 
force  anti  moment  expressions. 


Wing.  The  calculation  of  wing  aerodynamic  forces  and  moments  is 
relatively  straightforward.  In  Reference  8  a  definition  of  the  wing 
aerodynamic  coefficients  is  developed.  These  coefficients  contain 
interference  and  damping  effects.  The  aerodynamic  coefficients  for 
the  wing  in  accordance  with  Reference  8  are  as  follows: 


(Cd\i 

(CL>« 

<Vv 

(cl\i 


p 

^CD  )  ■"  ^7  CL  *  aw 

o  w  3cl  « 


(CL  )  ♦  (CL  )  .  aw 

0  w  a 


"Vw  •  8.  ♦  5C  •  Pw  *  (Cy_\, 


p  w  w 


r  w 


(C‘.,-^‘^(CVVP’'^(°*rVVV 


~ra  '  57“ 

o  w  w 


(c-  )  ♦  ^  (V  )  ’  «l  ♦  S?  •  6F 


(U.72) 

CU.73) 
(U.7U) 
6A  (U.75) 

(U.76) 
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"Vw  *  (Cn  >  •  6v  *  5T  (Cn>  •  P.  *  ST  (Cn  >  •  rw  *  °n  ‘  6A  (h'77) 

6  w  w  p  w  w  r  w  6A 

In  Figure  8  we  see  how  a  is  defined.  For  purposes  of  clarity  the 

wing  incidence  angle,  i  ,  has  been  exaggerated,  AV  is  the  induced 

velocity  previously  defined  in  equation  (I4.3U)  which  is  derived  from  wash 
effects  across  the  wing.  In  Figure  8,  we  see  that: 

2  (  2  2 
vw  •  (uw  ♦  AV)  *  vw  .  WK  . 


So  that: 


Sf  *  ^  ^uw “T  Avw'  ■  sin  vw  cos  ^ 


a. 78) 


and: 


Bw  «  tan-1  [ 


V(u  *  Av)Z  ♦  wf 


(U.79) 


In  equation'  (a. 78)  and  (a. 79)  vy  *  V  and  at  low  velocities  where  AV?*0, 


then  V 


v 


In  order  to  transform  wing  data  to  body  axis  we  rotate  through  an 
angle  where  £  «  (i*  -  \)  • 

Before  writing  force  and  moment  expressions  for  the  wing,  equations 
(U.72)  through  (a. 77)  are  transformed  to  body  axes  by  rotating  through 
the  angle  £  as  in  the  similar  analysis  for  the  XC-lli2A  (Section  IV). 


x'w 


z  n 


m  w 


sin  C  - 

(r  ) 

K  r'w 

cos  £ 

-  (CL) 

» cos  c 

-  CcD)w  sin  C 

(c  ) 

l  w 

COS  £  - 

(C  ) 

n  w 

sin  £  - 

(c  ) 

y  w 

(C)w 

m  w 

,  Za.c. 

(C  ) 

x  w 

• 

O 

• 

(C  ) 

z  w 

c 

c 

(c  ) 

x  n'w 

cos  £  ♦ 

<«»>* 

sin  £  ♦ 

(c  ) 

y  w 

a.c, 


■) 


CU.80) 

(a. si) 

CU.82) 

CU.83) 

a. 8a) 


x  and  z  are  the  respective  distances  from  the  c.g.  of  the  air- 

craft  in  aircraft  body  axes  to  the  aerodyn&nric  center  (a.c.)  of  the  wing 
in  the  x-z  plane,  c  is  the  mean  aerodynamic  chord. 

As  in  the  tilt-wing  example,  (Section  IV)  Equations  (a.80)  through 
(a. 8a)  are  readily  expressed  as  force  and  moment  contributions  of  the  wing. 
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<*Vw  '  <V.  S  % 

(U.85) 

(Jr,),,  *  Ccy)w  s  ^ 

a. 86) 

*  (<Vw  S  S, 

a. 8?) 

(4Lal-  <C  >w  b  S  ‘’w 

(U.88) 

'  KV  S  S, 

(U.89) 

<A,,a’w  *  (Vv  b  S  \ 

(U.90) 

In  equations  (U.85)  through  (U.90)  S  is  the  wing  area,  c  is  the 
mean  aerodynamic  chord^  b  is  the  wing  span  and  qw  is  the  wing  dynamic 

).  These  equations  will  be  incorporated  in  the 

total  force  and  moment  equations. 


pressure  (qw  *  1/2  p  V 


Ducted  Fan.  In  the  calculation  of  ducted  fan  aerodynamic 
forces  and  moments  wake,  interference  and  damping  are  included  in  tfce 
definition  of  the  aerodynamic  coefficients.  Six  coefficients  C_  ,  CM 

3  3 

C-  ,  C  ,  C.  ,  and  C  are  defined  in  Reference  8  to  describe  the 

S3  W)  h  *3 

behavior  of  the  ducted  fan.  These  coefficients  are  as  follows: 


cT  ■  c_  +  c  +  act 

1  *  JI  XIV 

3  3  3 


a.9i) 


CT  is  the  total  thrust  coefficient.  C_  is  the  change  in  thrust  due 

1  i  XT 

J  _  3 

to  varying  thrust,  CT  is  the  change  in  thrust  due  to  varying  inter- 

h 

ference  effects  and  ACT  is  the  change  in  thrust  due  to  the  inlet 

ivj 

vanes  to  the  duct.  j«*l  for  the  port  duct  and  j«*2  for  the  starboard 
duct. 


Cm  "  ♦  Cm  ♦  CM  +  Cw  ♦  C„  ♦  CM 

Ni  NF  nd  ni  nn  niv  ntv 

3  3  3  3  3  3  3 


(U.92) 


C„  is  the  total  normal  force  coefficient.  C„  results  from  fan 

N1  NF 

3  3 
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effects,  CM  results  from  duct  effects,  CM  comes  from  interference 

u  NI 

_  3  3 

effects,  C.j  is  the  change  in  normal  force  due  to  varying  normal  force, 

Nj 

CM  is  the  change  in  normal  force  due  to  the  inlet  vanes,  and  finally 

_Ivj 

CM  is  the  change  in  normal  force  due  to  the  trim  vanes. 


C«  m  C-  +  C«  +  Cq  C— 

>  FJ  D1  \  s 


(U.93) 


Cc  is  the  side  force  coefficient,,  C-  results  from  fan  effects, 

. J 

Ce  results  from  duct  effects,  C„  is  the  change  in  side  force  due  to 

s  .  \ 

varying  normal  force  and  G«  is  the  change  in  side  force  due  to  inlet 

ivj 


vanes, 


cl  "  cM  +  Cw  ♦  CM  ♦  ^  +  <L  +c„ 

^  %  >,  %  \  Avj  'S 

Cm  is  the  total  pitching  moment  coefficient  for  the  ducted  fan. 


a.9h) 


C„  ,  Cu  ,  C„  -  result  respectively  from  fan,  duct  and  interference 

”f  ^  ^ 

J  3  3 

effects. 

5^  is  the  change  in  pitching  moment  due  to  vaiying  normal  force, 

and  C}^  SJ"®  respectively  the  changes  in  pitching  moment  due 

3  3, 

to  inlet  and  trim  vanes. 

h.  ’  -  (-1,J  'it5  [5t  8in  1r  *  5n.  008  V  -(f  8inV)(Sn  8ln  V 

J  J  J  j  Dj 


X  XX 

-  [(-£)  cq  ♦  (-£)  C  +  (-—)  Ce  ]  Sin  in+  K..  sin  i  U.95) 
D  5f  D  Sd  D  5iv  d  1  a 

J  J  J 

Cr  is  the  total  rolling  moment  coefficient  for  the  ducted  fan. 

3 
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is  a  function  of  axial  advance  ratio  of  the  ducted  fan. 

%  ■  -(-D3(-|)[C  cos  lD  "  V  sin  iD1  -(7  sin  cos  V 

J  J  J  J  V  4 


Xp  _  Xy»  —  ^TU 

"  I  (-5)%  +  (T)CSn  *  (TT)CStw  ]cOS  iD+  K1  cos  iD  (U*96) 
r  *  D  .  I V  . 


j 


j 


j 


Cjj  is  the  total  yawing  moment  coefficient  for  the  ducted  fan. 

Before  writing  the  force  and  moment  equations  these  coefficients 
are  re-expressed  an  coefficients  representing  both  ducted  fans. 


"  K2  E  ^CT.  C0S  iD  “  CN.  sin  V 

D  j-l  J  J 


a. 97) 


(C.\ 


(V, 


"  K2  E  Sin  iD  +  C0S 

j-l  J  J 


K2  E  % 
j-l  J 


(.4-) 

TT  s 


a.98) 


a.99) 


(cv, 


-  K2  D/b  1  Cf 

J-l  **J 


CU.100) 


(0"V  ** D/c  &  \ 


a.  101) 


^n}D  -  K2  D/b  £  C1 


J-l 


a.  102) 


The  equations  for  duct  forces  and  momsnts  are  then  as  follows: 


(AX)  -  K,(C,  ) 

*  P  J  X  P 

(AYa)  -  K  (C  ) 
a  D  *  y  D 


(U.103) 

(U.10U) 
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(AZ  )  - 

a  D 

K3  K3  '  t  (’n  D'>  S 

(U.105) 

(AL  )  - 

a  D 

K,b(Ct) 

3  D 

(Li .  106) 

(4M»>  - 

K  c(Gm) 

(li. 307) 

D 

D 

(AN  )  - 

a  n 

K  t(C  ) 
i  71  n 

a.  108) 

Equations  (U.103)  to  (U.108)  will  be  included  in  the  total  aerodynamic 
force  and  moment  expressions. 


Vertical  Stabilizer  and  Rudder.  Forces  and  moments  for  the  vertical 
tail  and  rudder  will  be  defined  in  the  same  manner  as  for  the  XC-lli2A  and 
Reference  8.  The  associated  coefficients  are  as  follows: 


.  8  ♦  C 


.  6R 


5R 


C,  -  C  •  B  +  C  • 

1  lB  l6R  2VB 


♦  E*L.  c  ♦  c 


rb 


C  -  C  •  8  ♦  C  .  6R  +  Sr y—  C  +  Ary1  o 
n  nfl  n6R  np  V 


0i. 109) 

a. no) 

(mill) 


The  forces  and  moments  for  the  vertical  tail  can  then  be  expressed 
in  the  following  equations. 


(4Y  )  -  C  Sq  <3Ej 

VT  y  H 


4vt, 


(AL  )  -  C  bSq  (—•) 

VT  H 


^VT 

(AN  )  -  Cn  bSq  (-—) 

a  VT  q 


(1.112) 

(U.113) 


a.llli) 


This  force  and  these  moments  will  be  included  in  the  total  aerodynamic 
forces  and  moments. 

Horizontal  Stabilizer  and  Elevator.  The  aerodynamic  forces  and 
moments  for  the  horizontal  stabilizer  will  be  defined  in  a  similar 
manner  as  for  the  IC-1L2A.  Elevator  movement  is  considered  in  the  z-force, 
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and  the  pitching  moment.  The  ancle  of  attack  of  the  tail  will  he  a 
function  of  interference  damping  terms,  and  considerations  of  downwash. 
The  lift  and  drag  coefficients  of  the  tail  can  he  expressed  in  the 
following  relations. 


r 


.  a 


t 


(lull?) 


cD  ■  ♦  k  a. 116) 

% 


The  horizontal  stabilizer  (HS)  can  contribute  forces  in  the  x  and 
z  directions  and  a  pitching  moment.  The  equations  are  as  follows: 


'  [~\  co"  (it~at)  -  \  sin  (it“at)]  Sq(-r} 
HS  t*  ^ 


a.117) 


(AZJ 


HS 


[+CDt  Sin  (it’at}  “  CLt  cos  (Vat)] 


♦  qScC„  .  6E 


<“a> 


(U.118) 


HS 


*  hHS  *  (AZa)  *  lHS  +  TT"^  *  qlV  +  Cm.  •  ^ 

Ho  Ho  q^  a 


♦  qScC  .  6E  (U.119) 

6E 

^S  is  the  distance  from  the  aircraft  c.g.  to  the  aerodynamic  center 

of  the  horizontal  stabilizer  and  h„_  is  the  height  of  the  a.c.  above 

Ho 

the  c.g.  it  is  the  incidence  angle  of  the  tail.  Equations  (U .117 )  to 

(U.119)  will  be  included  in  the  total  aerodynamic  force  and  moment 
expressions. 

Fuselage.  In  a  very  direct  way  we  can  write  the  effects  of  the 


fuselage  (F)  on  the  total  aerocfynamic  forces  and  moments, 
follows: 

They  are  as 

(«Xj  ■  -q  S  CD 

F  o 

(U.120) 

(AY  )  •  M  S  C  .6 

a  f  ye 

(U.121) 

(AZ  )  -  -  S  CT  .a 

a. 122) 

F  a 


IMVTRADEV'EN  I/O1:-? 


(W,)  "  <5  Sc<Cm  ♦  cm  .  >) 


(an  )  =  q  Sh  c  .  e 

F  g 


(h.121) 

(lulPh) 


Equal  i  ons  (h.120)  to  (l;#l?h)  will  he  included  in  t,he  total  aerodynamic 
forces  and  moments. 

Reaction  Controls.  The  control  vanes  (horizontal  (VH)  and  vertical 
(W)  are  located  in  the  exhaust  flow  of  the  ,1  e t  engine,  The  reaction 
control  (RC)  vanes  provide  longitudinal  and  lateral  control  in  the  hover 
mode.  The  following  coefficients  are  defined  as  in  Reference  R. 


LVH  'VH 


CLW 


w 


Cdvh  "  K  Clvh 
C°w  "  K  C*w 


Tr  »  m  (VE  -  U)  - 


-  residual  thrust  (T^) 


vH 


JVH 


"RC 


\h 


a. 125) 


vv 


vv 


Hrc  v  H'V  v  v  Hyv 


Oi.126; 


XVH*  XVV*  ant^  •V’w  are  momen^  arms  to  control  vanes. 


(n  )  m  _  r  r  +  r 

v  /  l  J  n  o 


(O 


RC 


“VH  "W 

C. 


(C  ) 


-  c 


RC 


LW 

Wh 


(14.12?) 

(U.12P) 

a. i29) 


The  force  and  moment  equations  then  follow  from  equations  (b.125)  to 
(U.129)  as: 
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<«.) 


RC 


(AYa) 


(A"a> 


RC 


rc 


%  V  V  4 

RC 

qE  v' 

•  PC 


qF  3V^zc 


P  p.  =  .  OOP' 7 7  0  slugs  /f  t 


7p  ■  3*;h.‘  EPS 


7» 


i  ? 

qF  =  2  P F  VE 


(Ji.no) 
(ii.i  n) 

( J  i  •  1 3  2 ) 


(\’ftn  *  "K  Sv  C*  ^  ”m)-(C 

(ANa^  ■  %  WV^ 


(h .133) 
(Ji.131) 


Equations  0i,130)  to  (ii.l3h)  will  he  included  in  the  total  aerodynamic 
forces  and  moments. 


EQUATIONS  OF  MOTION  VZ-iiDA.  The  aerodynamic  force  and  moment  terms 
developed  for  the  VZ-UDA  for  each  of  the  major  aircraft  components 
will  be  combined  and  set  equal  to  their  corresponding  internal  force 
and  moment  expressions  developed  in  axis  system  section  for  the  VZ-UDA. 


Total  Aero  (dynamic  Forces  and  Moments.  These  expressions  are  as 
follows: 


*  (AX J  4 

a 

(AX  )  4 

a 

<“a> 

+ 

(AXa)  - 

(AXa) 

w 

D 

HS 

F 

RC 

-  (AY  )  4 

(AY  )  4 

B. 

(4Ya 

)  - 

(AYa)  4 

(AYa) 

w 

D 

VT 

F 

RC 

■  * 

(AC  )  4 

cl  __ 

)  + 

(AZa)  4 

(AC  ) 

,r\ 

w 

r 

HS 

F 

RC 

-  (AL  )  4 

(AI  )  4 

(•'! 

> 

/ 

a 

w 

a  r 

VT 

-  (AM)  4 

( AM  )  « 

(am 

)  4 

(AM  )  4 

(AM) 

a 

w 

a  r 

a 

HS 

1  F 

3  RC 

N  -  (AN  )  4  (AN  )  4  (AN)  ♦  (."')  4  (AN  ) 

w  D  a  VT  F  3  P7 


Equations  of  Motion,-  These  equations  are  subject  to  simplification 
in  the  interference  rerion  and  in  the  rerresentati on  of  duct  fan  forces 
and  moments  when  further  flight  test  data  is  available  for  the  determina¬ 
tion  of  typical  magnitudes  of  terms. 
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1.  X  Force  Equation 

(-C  cos  (it  -  at)  -  sin  (it  (it  -  at)  Sq(— ) 

U  t* 

♦  rx)  S  *  -CSS/,E  VCxL  ♦  Tr 

w  ™ 


D 


RC 


n(TI  +  Wq-j  -  Vr)  ♦  nyr  sin  6 


2.  Y  Force  Equation 


(c)  S  +  K.j(C  )  .Uk.i  v  ~ ~ 

y  *  s,  3  y  D  y  q  q 


s  (-XI)  +  q  s  C  . fl  ♦  q  S„(c  ) 

Q  q  yp  P  E  Vv  y'RC 


HS, 


■  m(V  +  Ur  -  Wp)  -  mg  cos  0  sin  4> 

3.  Z  Force  Equation 

[Cp  sin  (it  -  at)  -  CL  cos  (i  -  at)  ]  Sq(-~) 
t  t 

*  (C.)  *  k3(52)d 

*  q  Sc  C  .  6E  -  qS  CT  .  a  ♦  q_  S..(C  ) 

"SE  La  E  ^  Z  RC 

-  m(W  ♦  Vp  -  Uqx)  -  mg  cos  0  cos  4> 
hm  Roll  Equation 

-  ^VT 

(Ct)  bSq^  ♦  K3  b(c£)  ♦  Ct  b  Sq(-jji) 

w  D 

-  TnP  -  I13  (r  +  pqx)  ♦  (I33  -  I22) 

*  (Ie4)  cos  (iw  *  6„)  -  qx  (IEnE)  sin  (i„  ♦  6„) 

5.  Pitch  Equation 

(cm)  c  Sq^  4  K3c(C  )  ♦  q  Sc(Cm  ♦  C  .  a)  ♦  q  Sy  cv(Cm> 
w  D  o  a 

-(AV  *  hHS+  (AZa)„_#  lHS+  IT’  [Cm.  *  qlV  *  Cm.‘  *] 


RC 


HS 


HS 


♦  qSc  C_  •  6E 
oE 


1 22^1  -  1«<r*  ♦  *2>  *  (In  -  Jj3>  pr 
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♦  pCgfig)  sin  (iw  ♦  RD)  ♦  rdjPg)  cos  (iw  ♦  fip) 

6.  Yaw  Equation  » 

-  9\m* 

(C  )  bSqw  +  K  b(C  )  ♦  Cn  bSq(— ~— *)  ♦  q  Sb  C  .  g 

w  D  q  b 

*  %  SV  MBn> 

RC  , 

'  X33f  *  X13  *  qlr)  *  X22  -  hi'1  V**l 

-  (Ijr^g)  s^n  +  bp)  ”  9^  (T£*^;)  CCS  Ciy  ♦  6^) 

| 

No  simplifications  have  been  made  on  the  rights  side  of  these  six 
equations. 

\. 

TILT  PROPELLER 

The  Curtiss-Wripht  X-19  V/STOL  aircraft  will  be  used  as  the  exanqple 
in  the  development  of  tilt  propeller  (tilt  prop)  V/STOL  simulation 

equations  of  notloru _ The  development  which  follows  is  based  upon  the 

type  of  arguments  presented  for  the  tilt-wing  and  information  gleaned 
from  Reference  h. 

\ 

In  order  to  gain  some  idea  of  the  X-19  let  us  consider  Figure  3 
and  Table  6.  In  Figure  3,  a  three-view  arrangement  of  the  X-19  is 
shown.  Table  6  is  a  list  of  definition  of  symbols  contained  in  the 
X-19  equations. 

With  this  general  idea  of  the  X-19,  we  may  now  develop  the  mathe¬ 
matical  model  for  the  X-19.  First  we  define,  as  necessary  applicable 
axis  systems.  Second,  develor  aeroc$mamic  coefficients.  Finally,  write 
the  equations  of  motion  for  the  X-19. 


AXIS  SYSTEMS  FOR  X-19..  It  arpears  that  the  forward  and  aft  wing  chord 
arfe- parallel  to  the  x-bocfy  axis  of  the  aircraft.  In  this  aircraft  the 
aft  wing  functions  as  the  horizontal  stabilizer  with  elevator.  The 
total  velocity  of  the  wing,  V  ,  can  be  expressed ^easily  in  aircraft 

body  axes.  In  addition  separate  axis  systems  to  describe  propeller 
velocities  are  not  necessary  since  the  propeller  velocities  can  be 
expressed  in  terms  of  aircraft  boc$y  axis  velocities  (U,  V,  W).  Conse¬ 
quently,  no  additional  axis  systems  will  be  necessary  other  than  the 
conventional  inertial,  body,  stability  or  wind  axes. 

In  addition,  there  will  be  no  gyroscopic  effects  from  the  rotating 
masses  of  the  propellers  and  their  drive  shafts.  The  four  propellers 
are  -echarically  "irked  to  rotate  at  the  same  speed  and  the  starboard 
rrcr  lie.-  tur  1  r  .re  or:  os’*.-  d  recti  on  cf  the  sort  rropellers.  If 
*re  ~r  :  tj  tel  >.  fail*  •  -  -r  n.lr  -e  gyroscr r‘. -  effects. 


Reproduced  From 
Best  Available  Copy 
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Symbol 

F 

w 

n 

vt 

aF 

°Vf 

Sfa 

6f 
8wf 
6  wa 


V 

v 

c 

b 

S 


<W 

Sra 

q 

<*vt 

T 

Jn 

JA 

Nn  •  H 


ffci  ini  tier: 

Fuselage 
W  i  np 

Main  propellers  n  ■  1,  2,  Front.  n  ■  3,  ii  Aft. 

Tcp  view,  left  to  ripht 

vertical  Tail 

Angle  of  attack  -  fusel ape 
Angle  of  attack  -wing  forward 
Anrle  of  attack  -  wing  aft 
Sideslip  anple  -  fuselage 
Sideslip  angle  -  wing  forward 
Sideslip  angle  -  wing  aft 
Incidence  angle  -  propeller  nacelle 
Total  velocity  in  aircraft  body  axes 
Total  velocity  in  wind  stability  axes 
Mean  aerodynamic  chord 
Wing  span 
Total  area  of  wing 

Eynamic  pressure  at  forward  wing  due  to  power-on  effects 

Eynamic  pressure  at  aft  wing  due  to  power-on  effects 

Eynamic  pressure  -  free  stream 

Eynamic  pressure  at  vertical  tail 

Total  thrust  of  main  propellers 

Advance  ratio  main  propeller 

Advance  ratio  normal  to  propeller  disk 

RPM  -  main  propeller 

Table  6,  Definition  of  Symbols  Used  in  X-19  Equations 
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Symbol 

No 

EL 


N. 


n 


Def ini tion 

Nominal  RPM  -  min  propeller 
Blade  pitch  an^le  ••  main  propeller 
Main  propeller  thrust 

Main  propeller  thrust  component  normal  to  T  , 
Main  propeller  moment  (initially  pitching) 
Main  propeller  torque  coefficient 
Diameter  main  propeller 


Table  6.  Definition  of  Symbols  Used  in  X-19  Equations  (Cont*d.) 
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AERODYNAMIC  FORCES  AND  MOMENTS  -  X-19.  In  order  to  develop  expressions 
for  Xa(  Ya,  Za,  La,  M  ,  and  N  ,  contributions  from  the  major  airframe 

components  of  the  X-19  will  lie  considered  separately  as  was  done  for 
the  tilt-wing  aircraft.  The  major  components  to  lie  considered  are 
the  propellers,  the  forward  winp  and  aft  wings,  the  vertical  tail  and 
rudder,  and  the  fuselage. 


Propeller.  The  inclination  of  the  propeller  nacelles  from  the 
x-body  axis  is  defined  by  propeller  shafts,  angle  of  tilt,  i  .  In 

Figure  3  we  see  the  angle  i  and  the  peometry  defininp  the  total  pro¬ 
peller  velocity,  V^*  For  the  velocity,  V  ,  n  denotes  the  particular 

propeller  (n  «  1,  2,  3,  U).  The  front  left  propeller  is  number  1,  the 
front  right  propeller  is  number  2,  the  rear  left  propeller  is  number  3 
and  the  rear  right  propeller  is  number  U.  Each  propeller  nacelle  is 
located  by  x^,  yn  and  zn  body  axis  components.  Due  to  the  symmetric 

placement  of  the  propellers;  ■  x2»“x3  "  "XU*“^1  "  "  ^U>“zl 

and  ■  -2,  .  These  body  axes  moment  arms  (s  ,  y  and  z  )  are 
*  3  4  nF  n  n 

multiplied  by  the  appropriate  body  axis  angular  velocity  (p,  q^  and  r) 

in  order  to  give  tangential  velocity  components  of  the  propeller 
velocity  consonants. 


u_  ■  U  -  y (p  sin  i  ♦  r  cos  i  )  -  q, (x  sin  i  -a  cos  i  )  CU.135) 

n  p  pinpnp 

vn  -  v  +  xnr  (U.136) 

wn  ■  W  ♦  y^(p  cos  ip  -  r  sin  ip)  -  q1(xn  cos  ip+zn  sin  i  )  (U.137) 

We  will  now  define  the  propeller  aerodynamic  coeffients  C,j.  , 

n 

Cp  ,  and  .  These  coefficients  are  the  same  as  those  for  the 
n  n  n 

main  propellers  of  the  XC-142A,  however  there  is  no  CY  for  the  X-19, 

n 

since  the  port  prooellers  turn  in  an  opposite  direction  to  the  star¬ 
board  propellers.  First  let  us  define  the  advance  ratio  (Jn)  for 

each  propeller  and  the  advance  ratio  normal  to  the  propeller  disk 

CJA) 

60  V 


J' 

n 


jn  00.  *n 


(u.138) 
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The  symbol  N'n  (the  particular  propeller  RPM  and  since  all  propellers  ar< 

mechanically  linked  to  turn  at  the  same  rate)  The  number  6 0 

n 

chanpes  RT?!  to  RPS,  D  is  the  diameter  of  the  propeller  disk  and  n  is 

n 

defined  as  the  blade  pitch  ample  of  the  particular  propeller.  The  aero¬ 
dynamic  coefficients  can  then  be  developed  in  terms  of  advance  ratio 
and  blade  pitch.  '’Tie  anple  \|r  is  the  anple  between  the  propeller  thrust- 

vector  and  v  as  shown  in  Fipure  7.  ^  is  defined  the  same  as  for  the 

n  n 

tilt  winp  as  is  Sn — equations  (U.l)  and  (h.2). 


n 


acT  ac  s2c 

*  — L  j>  +  — L  n  +  — —— 

Tq  3J»  n  3  P  n  3B3J» 


P  J‘ 
n  n 


a. 139) 


Cr 


3C  3C 

'  J»  ♦  —L  F 


32C 


F 


3J» 


n 


30  n  3  B3J’ 


B  J* 
n  n 


(h  .lliO) 


'N 


n 


3  (C  cot  \|r) 

JL  [ _ 2 _ 

30  3  J* 


1  Bn  Jn  sln  *n 


(U.liil) 


'M 


n 


3  J 


3  C, 


Mi 


J* 

3  n  n 


(U.1U2) 


CT  ,  Cp  ,  CM  are  defined  similar  to  equations  (U.7),  (1.8),  (h.9) 
n  '  n  n 

and  (U.ll).  Since  blade  pitch  is  the  prime  control  parameter  in  the 
X-19,  the  second  partial  derivative  32  is  retained  in  equations 

3  2  3  B3J’ 

(ii.139)  and  Oi.lliO)  whereas  ye ,  g-  is  neglected.  is  the  coefficient 

*  n 

of  thrust  (Tn),  Cp  *s  the  coefficient  of  power  used  to  express  torque 

n 

(0>  C„  is  the  coefficient  of  normal  thrust  (N  )  and  is  the 
11  i.  n  m 

n  n 

longitudinal  propeller  hub  moment  coefficient,  (M^). 


Before  expressing  the  force  and  moment  contributions  in  bocjy  axes 
due  to  the  propellers,  we  will  define  Tn,  N^,  and  (Reference  6). 


T 


n 


&  (A  (-£-) 


N 


(U.1U3) 


N 


n 


*  <A  (T>  = 

‘  o  % 


N 

n 


a.iiii*) 
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*n 


-  D5  (A  C— )  C, 


FP  '77  ~m 

0o  r 


«n  •  r  (r»  <f>  CP 

O  T) 


n 


(U.U45) 

(li.lU6) 


NQ  is  the  maximum  RPM  of  the  propellers  and  pq  is  the  air  density  at 
sea  level  on  a  standard  day. 

From  equations  (1*.U,3)  through  (U.l|j6)  we  can  then  write  the  pro¬ 
peller  force  and  moment  expressions  in  a  manner  similar  to  equations 
(U.17^  through  (U.22)  of  the  XC-lii2A. 


(AX  )  «  £  (T  cos  i  -  N  cos  6  sin  i  ) 


a  P 


n-1 


n  p  n  n 


(AY  )  ■  r  (-N  sin  8  ) 

v  a'p  “  1  n  n' 

K  n-1 


(«,jp  *  t  (-Tn  8ln  ip  -  Nn  008  Bn  008  ip> 
K  n-1  r 


(AVp 


-  I (AZ  )  -  (AZ  )  ]  -  [(AZ  )  -  (AZ  )  L 

a  p!  a  p2  yl  a  p3  PU  y3 


♦[(AT)  ♦  (AT  )  J  +[(AY)  ♦  (AY )  1 

a  p1  a  p2  z1  a  p^  a  p^  z^ 


-  r  M_  sin  B„  sin  i 
n-1  n  n  p 


(U.iii7) 

a.us) 

(U.U*9) 


a.  150) 


(“.Ip  •  *i(Ti  *  V  !ln  ip  -  x3(T3  *  V  sin  ip 

-  21(t1  ♦  Tg)  COS  i  -  *3(T3  ♦  )  cos  ip 

♦  z^(N^  cos  6^  ♦  Ng  cos  Bg)  sin  i 

cos  By  *  cos  8^)  sin  ip  (lul5l) 

♦  x^(N^  cos  6^  ♦  Ng  cos  Bg)  cos  ip 
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-  x^(N^  cos  cos  8^)  cos  i 


nil 


M  cos  8„ 
n  n 


<AVp 


f(X  )  -  (XJ 


a  Px  a  P2  yl  a  P 


♦  t(x  )  -  (XJ„  ] 


ph  y3 


♦  [0„  +  (Y.)n  1  v  -  [(*.)«  *  <YJn  ]  * 

a  Pi  a  p2  *1  a  p3  a  ph  x3 


a. 152) 


-  r  M  sin  8  cos  i 
S  n  n  j 

n-1  r 


Note  that  xj»  y^,  y 3,  z^  and  z^  are  positive  values.  Th  propeller 

force  and  moment  equations  (li.lU7)  through  (U.152)  will  be  included  in 
the  total  aerodynamic  forces  and  moments. 


Forward  Wing.  The  forward  wing  forces  and  moments  are  developed 
as  for  the  XC-1U2A  in  Section  IV  .  It  is  assumed  that  the  wing 

incidence  angle,  is  zero  and  that  the  induced  velocity,  AV,  is 

defined  as  equation  (U.3U)  when  i  is  zero.  AV  will  change  as  a 

P 


function  cos 


and  flap  position  when  the  aircraft  i3  in  hover  and 


transition  regions.  Consequently,  as  a  preliminary  estimate,  the  wing 
velocity  u^  is  as  follows: 


uwf  "  up  *  (AV)TP  CU .153 ) 

where  (AVj^p  «  K  AV  cos  i 


For  forward  flight  where  i  equals  zero  then  u  equals  V_  and  K  »  1, 

p  p  B 

K  is  a  function  of  flaps  and  interference  effects  as  the  propeller 
wash  hits  the  wing  at  high  i  angles.  K  and  cos  will  tend  to 

decrease  AV  as  i^  increases.  As  in  equation  (b.36)  q^  is  defined 

similarly  as: 


v 


,<T,  ♦  T,) 

(q  »  1  ,  ;  ) 

2TTir 


(Ji.lSU) 


The  forward  wing  aerodynamic  coefficients  are  Cp,  (C^)^, 
(Cm)wf  ^  ^Cn^wf  CE  ^  CL  define  ^Cx^wf  and  ^Cz^wf  resP*ctively» 
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Since  the  deveii  jpment  of  these  coefficients  is  in  wing  axes,  rolling 
and  turning  rates  necessary  to  define  these  coefficients  are  transformed 
from  bocty  axes.  The  only  difference  is  in  the  angle  of  attack  of  the 
wing,  which  is  from  equation  (h.3<?)  as  follows: 


1  W 

a  -  -  tan"  [ - P r  ...i — 

wf  Up  ♦  (AV)Tp 


(U.155) 


where  w  _  -  w  ;  w  -  u  sin  i 
wf  p  p  p 

and  u^  -  u  ♦  CAV )Tp ;  U  -  u  cos  i 


In  a  similar  manner  6^  is  defined  as: 


-1  V 

6“r  '  tan  vr 

Wf 


CU .156) 


where  Vp  *  ■  V 

and  Vwf  -  [W2  ♦  (U  ♦  AV)Tp)2  ♦  V2]  1/2 


We  write  for  the  wing  rolling  rate  (p^ )  and  the  wing  turning  rate 
(r^)  the  following  equations: 


wf 


*  p  cos  +  r  sin  a 


wf 


(U.157) 


rwf  ■  -p  sin  q^  ♦  r  cos  a, 


’wf 


(U.15B) 


These  are  just  equations  (U.37)  and  (U,3B)  where  £  -  since 

iw  •  0  in  ^  -  Ciw  -  o^). 

Following  the  tilt  wing  the  aeroctynamic  coefficients  for  the 
wing  are  as  follows: 


CD  "  CD  +  CD  *  6F 
o  SF 


Cl  ’  \  *  %f  ^  ^  *  \  ‘  ^ 


(U.159) 

(U.160) 


Cp  and  do  not  reflect  as  strong  a  flap  dependence  as  in  the  XC-lli2A. 
This  is  a  preliminary  estimate  of  Gp  and  Cj  which  may  be  modified  by 
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further  manufacturer  data.  (C  )wf  ^m^wf  ani3  ^n^wf  defined 

as  fnr  the  XC-llt2A. 


^  t>wf 


•  «wf  *  •  6A  *  2—  c,_-  }wf  ‘  5TT  c.  •  rwf  a-I61) 


•g  '6.  ct/wf  *p  ^wf  *r 


(c-)-  ‘  -  *  c-  *  6F  +  2C  V  •  qi 


'm  wf  'm  mK„ 

o  or1 


a. 162) 


K, JU 

wf  q 


(r  ) 


n'wf  "  Cn  *  8wf  +  Cn  *  6A  *  2FT  Cn  *  Pwf  *  2V“V  Cn  *  rwf  (1*.l63) 
gp  bA  wf  p  wf  r 


Equations  (U.159)  through  (U.163)  will  now  be  written  in  body  axes 
by  rotating  through  an  angle  (-a  )• 


^CzW 

^C£  W 

^m^wf 


-  CD  cos  ♦  CL  sin 


-  Crj  sin  cos 


(Vwf  cos  awf  ’  (Cn}wf  Sin  awf 


(CmKf  “  c  (CxW 


ir  <Vwf 


a.i6L) 

a. 165) 

a.i66) 

a .167) 


(Cn}wf  ’  (Vvf  sin  awf  +  (Cn}wf  cos  awf  (h-l66) 

X  and  Z  are  the  respective  distance  from  the  c.g.  of  the  aircraft 

SC  Cl  c 

in  body  axes  to  the  aerodynamic  center  (a.c.)  of  the  forward  wing  in 
the  x-z  plane,  c  is  the  mean  aerodynamic  chord  of  the  forward  wing. 

The  wing  force  and  moment  contributions  can  then  be  written  from 
the  coefficients  expressed  in  equations  (U.16ii)  through  (U.166).  S  is 
the  forward  wing  area  and  the  dynamic  pressure  at  the  wing. 


(Vvf 

s 

V 

a.  i69) 

<4Za>vf  ■ 

"Vwf 

s 

(U.170) 

(ALa)wf 

<V«f 

bS 

a. i7i) 

■ 

(Cm>wf 

cS 

a.  i72) 

(iV*r  ■ 

"nV 

bS 

-Stf _ 

a. 173) 
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Vertical  Stabilizer  and  Rudder.  The  discussion  is  the  same  as  for 
the  XC-II42A  since  there  is  obvious  similarity  in  configuration  of  the 
vertical  stabilizer  for  the  propeller  wash  effects  on  the  vertical 
stabilizer  for  the  XC-lii2A  and  the  X-19.  Forces  and  moments  for  verti¬ 
cal  tail  (vt)  and  rudder  arise  from  the  relative  wind  pushing  against 
the  vertical  tail  surfaces  thereby  causing  a  turning  moment  due  to 
control  input  to  the  rudder.  This  gives  side  force,  as  well  as  rolling 
and  turning  moments. 


Wind  pushing  against  the  vertical  stabilizer  and  rudder  yields  a 
side  force  and  rolling  moment  (ALa)v^.  respectively  which  are 

non-dimens ionalized  in  terms  of  aerocfynamic  coefficients  as  G  and  Cj . 

The  rolling  moment  is  coupled  in  rolling  velocity  p  and  turning  velocity 
r.  We  can  then  write  for  C  and  G  the  following  expressions* 

J  ' 


C 


0  ♦  C 


6R 


fiR 


(U.17U) 


C  is  the  change  in  side  force  coefficient  with  changing  sideslip 
y  6 

angle.  It  acts  as  a  damping  term.  C  is  the  change  in  side  force 

ygp 

coefficient  with  rudder  deflection  and  represents  the  controllable  terra 
in  the  expression.  This  is  important  in  the  use  of  Automatic  Stabiliza¬ 
tion  Equipment  CASE). 


Cj  ®  8  ♦  C 


6R  ♦ 


16R 


pb  n  rb  - 

2v~  W~  # 
*VB  p  *VB  *r 


ca.175) 


C,  is  the  change  in  rolling  moment  with  variation  in  sideslip  angle. 

e 

The  vertical  tail  is  normally  above  the  X  axis,  hence  a  side  force 
on  the  tail  gives  a  rolling  moment,  C  is  the  change  in  rolling 

l6R 

moment  coefficient  due  to  rudder  deflection.  C  is  the  roll  dairying 

P 

derivative.  C  is  the  change  in  rolling  moment  coefficient  with 
tr 

change  in  yawing  velocity. 

Rudder  deflection  will  give  a  turning  moment  (ANa)v_^  which  can  be 


expressed  by  the  aerodynamic  coefficient  C^. 


C  -  C 
n  n„ 


.  &*  C  .  6R  ♦  C  ♦ 


6R 


2V 


B 


n 


2V 


B 


(U.176) 


C  ia  the  weathercock  or  static  directional  derivative.  C  is 

B  n6R 

rudder  effectiveness,  the  change  in  yawing  moment  coefficient  with 
**dder  deflection. — C^.  la  the  change  In  rawing  moment  coefficient 

P 
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I 


I 


with  varying  rolling  velocity.  Cn  is  the  yaw  danping  derivative. 

r 

The  tail  is  the  main  contribution  to  C„  . 

n 


The  forces  and  moments  for  the  vertical  tail  can  then  be  ex¬ 
pressed  in  the  following  equations: 


(AYJ 


a'vt 


C  S  (— ) 
y  q  q 


(AVyt '  't  iS 


(AN  )  .  C  bSq(— ) 

a'vt  n  q 


U.177) 
(U.178) 
a. 179) 


Here  q  is  the  free  stream  dynamic  pressure  and  is  the  vertical 

tail  dynamic  pressure.  These  force  and  moments  will  be  included  in 
the  total  aerodynamic  forces  and  moments. 

Aft  Wing.  The  aft  wing  (wa)  has  two  propellers  placed  in  a 
similar  configuration  as  the  forward  wing  and  acts  as  a  combination 
wing  and  horizontal  stabilizer.  When  the  four  propellers  are 
simultaneously  rotated  by  the  angle,  i^,  the  aft  wing  is  either  in 

or  out  of  the  front  propeller  wash.  If  the  aft  wing  is  out  of  the 
front  propeller  wash,  the  aft  wing  will  be  practically  a  duplicate 
of  the  forward  wing  in  aerodynamics.  If  the  aft  wing  is  in  the 
front  propeller  wash,  the  aft  wing  has  induced  velocity  effects  from 
its  own  and  the  front  propellers.  The  description  of  downwash,  as 
for  the  XC-1U2A  Horizontal  Stabilizer  is  not  applicable  to  the  X-19. 
However,  since  the  aft  wing  has  elevator  control  and  no  aileron 
control,  the  forces  and  moments  developed  will  lie  in  the  x-z  body 
axes  plane  as  for  a  horizontal  stabilizer. 


As  for  the  forward  wing,  the  aft  wing  incidence  angle  is 
assumed  to  be  zero.  Induced  velocity  effects  for  the  aft  wing  are 
as  follows: 

(U.180) 


(14.181) 


'S/a 

S/a 


Up  ♦  UV)TP 

T,  ♦  T«  ♦  T,  ♦  T, 

(q  ♦  - - 


wD 


Equation  (U.180)  is  the  same  as  (h.153),  but  is  for  the  aft  wing  out 
of  front  propeller  wash.  Vw  is  the  same  as  V^.  out  of  wash,  and 

Vwa  differs  from  Vvf  bX  the  e^ect  of  induced  velocity  in  wash  effects. 

For  the  aft  wing  in  front  propeller  wash,  equation  (U.180)  becomes: 


r 
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uwl  •  up  *  3/?(4V)TP 


(it. 182) 


where  (&V)^p 
Notice  that  if  i 


K0  AV  cos  i 

2  p 


3K? 

o,  (av)tp  -  -r 


AV  - 


3AV 

~r 


(K  ~  1  at  i 
2  P 


0) 


Kg  accounts  for  interference  effects.  From  momentum  theory  the  induced 

velocity,  AV,  does  not  double  on  the  aft  wing  from  the  effect  of  two 

propellers  blowing  on  the  wing  in  parallel  but  AV  increases  by  one-half. 

The  coefficients  (Cj,  .  (C  )  and  (C  )  are  defined  similar  to 

u  wa  l  wa  m  wa 

equations  (U.159),  (U.160)  and  (U.162),  but  using  elevator  controls 
instead  of  flap  controls. 

'  CD  *  CDst.  •  6E  (1*-183) 

o  ob 


(CLJwa 


♦  C, 


6E  ♦  C. 


J6E 


°F 


wa 


(L.18U) 


m 


c 

5v 

wa 


rn 


Qn  ♦  C 


m- 


W 


(U.185) 


The  coefficients  (Cjj)wa» 
axes  by  the  aft  wing  angle  of 
axes  are  similar  to  equations 


(C, )  and  (C  )  „  are  rotated  to  body 
l  wa  m  wa 

attack  (a  ).  The  coefficients  in  body 
wa' 

(U.16L),  (U. 165)  and  (1.167). 


(cr) 


x'wa 


-  (CD>m  003  awa  * 


sin  a 


wa 


(U.186) 


-  (Vi* 


sin  a 


wa 


“  (fywa 


cos  a. 


wa 


CU .187 ) 


From  the  coefficients  defined  in  equations  (U.185),  (U.186)  and 
(U.187)  the  force  and  moment  contributions  of  the  aft  wing  are  as 
follows : 


■  (Cx>w.  S  S* 

^AZa\ra  "  ^Cz^wa  ®  Sfa 

Z  X 

(AM  )  -  (C  ) - —  (AX  )  ♦  —  (AZ  ) 

v  a'va  ra'wa  c  a'va  c  a'va 


(U  *188) 
(U .189) 

(U.190) 


In  equation  (U.190),  X  and  Z  are  the  respective  distances  from  the 

ac  ac 

aircraft  c.g.  to  the  aerodynamic  center  (ac)  of  the  aft  wing  in  the 
x-z  plane,  c  is  the  mean  aerodynamic  chord  of  the  aft  wing. 


NAVTRAEEVCEN  1205-2 


Fuselage.  In  a  very  direct  manner  we  can  write  the  effects  of  the 
fuselage  (F)  on  the  total  aerodynamic  forces  and  moments  as  for  the 
XC-U2A. 

We  have  for  the  forces 


(“.)  •  -  \o  v?  s  cD 

F  o 


a. 191) 


Cp  is  the  equilibrium  drag  coefficient. 


(4Y  )  •  ♦  f  pV2b  S  Cy  .  6F 

F  ^ 


CU.192) 


C  is  the  change  in  side  force  with  respect  to  a  changing  sideslip 
angle. 


(«.>  -  *  f  »VI  S  CL  •  “F 

F  aF 


(U.193) 


C  is  the  change  in  lift  coefficient  with  varying  angle  of  attack. 

This  is  also  known  as  the  lift  curve  slope. 

We  have  for  the  moments: 


(ALa)  -  0 


'  I  0  VB  S  c  Cm  *  ?  0  VB  S  c  °m 
F  o  a 


‘  aF 


a.i9U) 


C  is  the  aerodynamic  pitching  moment  coefficient  in  equilibrium 
mo 

flight  and  Cm  is  the  longitudinal  static  stability  derivative. 
aF 

(ANft)  -  \  o  v|  S  b  Cn  .  8f  Ch. 195) 

F  Bp 

xs  the  static  directional  or  "weathercock''  derivative. 

1  2 

In  the  above  expressions  q  »  0  Vg  which  is  the  dynamic  pressure, 

the  forward  wing  span,  c  is  the  mean  aerodynamic  chord  and  S  is 
the  forward  wing  area. 


n 


h 


b  is 
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For  hovering  and  in  transition  regions  8p  is  assumed  small  so  that 
cos  &p  "  1. 

EQUATIONS  OF  MOTION  -  X-19.  The  aerodynamic  force  and  moment  terms 
developed  in  this  section  for  each  of  the  major  aircraft  conporents  of 
the  X-19  will  be  combined  and  finally  expressed  in  the  equations  of 
motion. 

Total  Aerodynamic  Forces  and  Moments,  The  total  forces  and  moments 
are  as  follows: 


X.  -  (4Xa)p  *  <«,)„,  ♦  (Ma)wa  ♦  («a)F 

Ta  -  («a)p  *  <4Vvt  *  (4Vf 

Za  -  <«a)p  ♦  («,).*  ♦  (42a)„a  *  (AZa)p 

la  -  (4La)p  .  (41a)rf  ♦  (iLa)rt 

Ma  ■  (“.Ip  *  <“a>wf  *  *  (4Ma>F 

N.  -  («Vp  *  (4Vvf  *  (4Na>Tt  *  f4Na>F 

Equations  of  Motion  Expanded.  The  forces  and  moments  are  expressed 
in  body  axes.  These  equations  are  subject  to  revision  when  further  data 
are  available  concerning  the  X-19. 

1.  X  Force  Equation. 

U 

K  (Tn  co5  V*n  cos  Bn  eln  V  *  (Cx>wf  ^ 

n-i 

*  (CA.  SS«  -  CD  Sq 


-  m  (U  ♦  W„  -  Vr)  ♦  mg  sin  0 
ql 


2.  7  Force  Equation. 

U 

E  (-N  sin  8  )  ♦  C  Sq.  ♦  Sq  C 
n-i  n  n  7  ~Vt  y^ 

-  m(V  *  Vr  -  WJ?)  -  mg  cos  0  sin  ♦ 


8 


F 
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3.  Z  Force  Equation. 

J  (-Tn  sin  lp  -  Nn  cos  en  cos  1  )  ♦  Sq^ 


*  «.*»  Sqwf  '  S<!  CL  •  aF 

aF 

C 

■  m(W  ♦  Vp  -  Oq-^)  -  mg  cos  0  cos  ♦ 
h.  Roll  Equation. 

-[(Tx  -  Tg)  sin  lp  ♦  0^  cos  ^  -  N?  cos  Bg)  cos  i  J  y1 

♦[  (T^  -  T^)  sin  ^  ♦  (N^  cos  Bj  -  cos  8^)  cos  ip] 

-[Nj  sin  81  ♦  N2  sin  Bgj  *3  ”  N3  3in  8  3  *  ^  sinB^J  *3 

h 

.  ^  Mn  sin  6n  sin  ip  MO  bSq^  ♦  C  bSq^ 

+  J13  ■*  ^33  *  ^22^  ^1** 

-  IUP 


5.  Pitch  Equation. 

[xi(Ti  ♦  T2)  -  x^T^  ♦  T^)}  sin  ip-  [sjCTj+Tj)*  z^T^T^O  cos  i^ 

♦  [z1(N],  cos  N2  cos  Bg)  ♦  *3(^3  cos  6^  ♦  cos  8^)]  sin  i 

♦[x^N^  cos  Ng  cos  8j)  -  x^fN^  cos  83  ♦  cos  0^)]  cos  ip 

*  I  Mn  008  6n  +  oSV  *  (Cm>,rf  cS^cSfC.  ♦  CB  .  y) 


o  a. 


*  Ii3  (r2  -  p2)  -  (Iu  -  I33)  pr 


X22ql 
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Yaw  Equation. 

U\  - 

Tp)  cos  i 

♦  (Np  cos  Pp 

-  N1 

cos  p.^) 

sin 

V,yi 

((t3  - 

T,  )  cos  i 
h  P 

*  <Nl4  C°S  h 

-  N3 

COS  P3) 

sin 

V)y3 

*  (n3 

sin  B..  *  N. 

3  h 

sin  p^)  - 

<N1 

sin  p^  ♦ 

n2 

sin  Pp)  x. 

-  ^ 
n=l 

M  sin  B 
n  rn 

cos  5  *  ‘Vwf  bSqwf  *  Cn 

bSq 

♦  qBSC 
vt  H  1 

+  *13^  -  qxr)  "  (-122  ~  Xll^ 


External  stores,  rough  air  or  landing  gear  conditions  are  not  included 
in  these  equations. 
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FAN-IN-WINO 


The  Oeneral  Electric/ltyan  Aeronautical  XV-5A  V/STOL  airplane  will 
be  used  as  the  example  in  the  development  of  fan-in-wing  V/STOL  simula¬ 
tion  equations  of  motion.  The  equations  of  Section  TIT  will  be  directly 
applicable  to  the  simulation  equations  for  this  aircraft. 


In  order  to  obtain  some  idea  of  the  XV-5A  consider  Figure  la  and 
Table  7.  In  Figure  U,  a  three-view  arrangement  of  the  XV-UA  is  shown. 
Table  7  contains  some  of  the  physical  characteristics  of  the  aircraft. 

With  this  general  idea  of  the  XV-5A,  we  may  now  develop  the  mathe¬ 
matical  model  for  the  XV-F A.  First,  we  will  define,  if  necessary, 
applicable  axis  systems.  Second,  we  develop  additional  aerorjynamic 
coefficients.  Finally,  we  write  the  equations  of  motion  for  XV-5A.. 


AXIS  SYSTEMS  FOR  THE  XV-5A.  The  XV-5A  has  no  change  in  its  physical 
corf iguration  during  the  hover,  transition  and  normal  flight.  The 
only  gross  change  from  hover  to  normal  flight  is  from  thrust  gener¬ 
ated  from  the  nose  and  wing  fans  to  thrust  from  the  tail  pipes.  The 
thrust  developed  by  the  fans  during  hover  and  transition  should  be 
available  as  a  function  of  fan  revolutions  per  second  (rps)  and  fan 
louver  position  which  directs  the  thrust  of  the  fans.  The  thrust  of 

the  aircraft  is  described  sufficiently  in  the  axis  systems  used  for 
ordinary j et  aircraft  simulation  (Reference  2).  Consequently,  no  axis 
systems  in  addition  to  the  inertial,  bod&r,  stability  or  wind  axes  will 
be  necessary,  to  develop  equations  of  motion  for  the  XV-5A* 

In  addition,  there  will  be  no  gyroscopic  effects  from  the  wing  fans 
since  they  remain  in  the  same  geometric  plane,  are  free  wheeling  in 
opposite  directions  and  are  driven  from  a  common  gas  source  at  approxi¬ 
mately  the  same  rps.  The  fan  in  the  nose,  also  free  wheeling,  does  not 
give  any  appreciable  lift,  but  is  used  as  a  pitch  control.  The  doors 
in  the  fuselage  ventral  to  the  nose  fan  are  linked  to  the  control 
stick  in  order  to  control  the  flow  of  air  thrcugh  the  nose  fan  thus 
giving  pitch  control. 


tnf  ^F  qi 


,  L  term. 

a 


I 


NF 


M  term, 
a 


I 


NF  fiNF 


N  term. 

a 


It  may  be  determined  after  consideration  of  data  from  the  manufacturer 
that  these  terms  will  be  negligible,  but  for  the  present  they  will  be 
included  in  the  equations  of  motion. 
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Weight 

Empty- 

Design  Gross 

Overload  Ferry  (Conventional  Take  off) 
Dimensions 

Wing  Span 

Fuselage  Length  (over-all) 

Height 

Powerplant 

Lift  Fan  System  -  G.E.  X353-5B 

Fan  Diameter,  Installed 
Fan  Thickness 

Turbojet  -  G.E.  J85-5 

Over-all  Length  with  Diverter  Valve 
Maximum  Diameter 


9,200  lbs. 
12,000  lbs. 
13,600  lbs. 


29.8  ft. 
UU.5  ft. 
lli. 7  ft. 


76  inches 
lli .5  inches 


67.2  inches 
17.7  inches 


Table  7.  Selected  Physical  Characteristics  -  TV -5k 

AERODYNAMIC  FORCES  AND  MOMENTS  XV-5A.  It  is  proposed  that  since  the  TV -5k 
haa  no  special  problems  in  development  of  the  equations  of  motions  that 
equations  such  as  (3.7)  through  (3.12)  or  those  by  Connelly,  in  Reference  2, 
be  used  for  simulation  purposes.  In  equations  3.7  and  3.9  the  engine 
thrust  will  be  defined  to  conform  with  the  fan-in-wing. 

The  total  engine  thrust,  T,  is  composed  of  either  fan  thrust,  Tp, 
or  Jet  thrust  Tj. 

Tp  ■  T„  ♦  T-.  1  is  port  fan  (U-196) 

1  r2  2  is  starboard  fan 

Due  to  cockpit  control  of  the  diverter  valve,  T  is  equal  to  Tp  or  Tj. 

For  example,  when  the  aircraft  has  enough  forward  speed  to  be  supported 
on  wing  lift  the  diverter  valve  is  moved  so  that  all  the  incoming  air 
passes  through  the  Jet  engine  for  forward  thrust  (i.e.  Tj)  and  the  "an 

inlet  and  exits  are  closed.  There  is  a  three  to  four  second  delay 
before  the  fan  develops  100*  lift  when  the  diverter  valve  is  writched 
from  Tj  to  Tp.  This  delay  is  easily  simulated  and  rids  the  thrust  of  a 

discontinuity  between  T^  and  Tp. 
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Hover  Flight  XV- 5 A.  The  equations  for  hover  flight  are  equations 
(3.7)  through  (3.12)  vrith  the  aerodynamic  coefficients  considered 
negligible .  We  can  then  write  the  following  expressions: 

n  (TJ  +  Wq^  -  Vr)  •  cos  a,p  ♦  Tp  sin  0^ 


m  (V  +  ,Tr  -  Wp)  »  g  cos  0  sin  f 


m  (W  ♦  Vp  -  Hq^)  --Tj  sin  aT  -  Tp  cos  0^  ♦  f (T^p) 

♦  ? —  (Tp  cos  9t  "  cos  9t  ) 

xn  1  Tn  Fi  h  f2 


I13  •  IMP  %F 
V  “  T^~  r  “  • 
jll 


*  ^  CL  .  6A 
X11  6A 


t  „  m 

I 


NF  ^NF 


*NF  TNF  C0S  eNL  t  qSc  c 


22 


L22 


*22  “bE 


6E 


r  .  hi  p  +  . ^  (Tp  sin  0T  -  TP  sin  0T  ) 


I 


33 


33 


*33  '  F1 


♦ 


T  C  .  6R 
T33  6R 


C.'ur- ) 


a.  iq?) 


a. 200) 


(ii.201) 


In  these  equations  0L  is  the  angle  the  fan  thrust  is  rotated  (0^  is 

a  function  of  louver  angle)  and  y  is  the  moment  arm  from  the  x-axis 
(parallel  to  the  y-axis)  to  the  center  of  either  wing  fan.  0.  and 

0T  refer  to  louver  angles  at  the  port  and  starboard  fans.  Cross- 

coupling  of  angular  rates  is  assumed  negligible  for  the  hover  condition. 
"Hie  gyroscopic  effects  of  the  nose  fan  (NF;  are  included  in  equations 
(h.109),  (h.200)  and  (It. 201).  These  gyroscopic  terms  will  be  zero  in 
conventional  flight.  Reaction  moments  are  included  in  equations  (It. 19°), 
(It. 200)  and  (U.201)  to  account  for  aileron  (6A),  elevator  (6E)  and 
rudder  (f>R)  control  motion  during  hover.  The  reaction  moments  are  con¬ 
tained  in  C  ,  C  and  C  so  that  they  phase  naturally  into  the 
1  bA  m6E  n6R 

transition  and  conventional  flight  equations  of  motion.  Pitch  control 
is  represented  by  nose  fan  thrust  (TNF).  The  pitching  moment  created 
bv  the  nose  fan  is  fd^).  tT  where  is  the  moment  arm  from  the 

center  of  the  nose  fan  to  the  c.g.  of  the  aircraft. 
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t 


Transitional  Flight.  As  in  the  equations  for  transitional  flight 
!'or  the  T.1127  (Section  1^ )  ,  the  transitional  equations  for  the  XV- r  A 

are  composed  of  the  equations  (3.7)  through  (3»12).  Consequently,  w»* 
can  write  the  equations  of  motion  for  the  fan-in-winp  directly. 

F  'NATIONS  OF  MOTION  XV-5A,  The  equations  of  motion  are  expressed  in 
a.rcraft  bc<ty  axes. 


X  Force  Equation 


m  (n  ♦  Wq-^  -  Vr)  * 


0VBS 


(C  (a,  Ma)  +  C  (B)  +  C  .  6F] 
x  x  X6F 


cos  dp  +  T.,.  sin  0^  -  mp  sin  0 


Y  Force  Equation 


2 

vis 


•  '  TV"' 

n  (V  *  Ur  -  Wp)  ■  [C  .  0  ♦  C  .  6R]  +  mg  cos  0  sin 


2  y 


6R 


Z  Force  Equation 


V?s 


m  (W  ♦  Vp  -  Uq-,)  -  £-|~  [C  (a,  Ma)  ♦  C  .  6F  ♦  C  .  6E) 


'6F 


'6E 


-  Tj  sin  aR  -  Tp  cos  0^  ♦  f(Tjjp)  ♦  mg  cos  0  cos  <t> 


Roll  Equation 
I 

”n 


•  13  /♦  .  ,  zI33  "  I22n  _  INF 

v  -  7 —  (r  +  pq-, )  -  (-^p - )  qnr  ♦  — - q. 


fcn 


Ln 


Sb 

♦  £*£= —  [c  .  e  ♦  c  .  6A  ♦  C,  .  5R  +  (c 
£  ill  £0  1 6A  oR 


^  ^ 


p  ♦  .  r)] 


(T„  cos  e.  -  Tp  cos  0T  ) 


T11  '  F1 


4 
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Pitch  Equation 

T  T  -  T  T  A 

i13  /_2  _2,  /ll  x33>  __  NF  NF  _ 

q  -  - —  (r  *  p  )  -  ( - - - )  pr  -  — = - p 

i  I22  I22  I22 

vise 


♦  D5T—  (a,  Ha)  ♦  Cw  .  6F  4  C  .  6E  ♦  gL  (C  .q,  ♦  0  .  .  a)] 
“22  m  m6F  m5E  2VB  V  1  ma 


22 

f^TNp)  *? 
t22 


Yaw  Equation 

•  T13  /•*  _  >  /I22  “  I11 
r  -  =-=■  (p  -  Qnr)  -  ( - T - 

J33  X33 

v|sb 


)  pqx 


i 

1NF  NF 
X33 


+  [C  .  B  +  c  .  6A  +  C  .  6R  ♦  5§~  (C 

2I33  n6  6A  n6R  **B  "p 


P  +  C  •  r)l 

r 


♦  JL.  (t  sin  e  -  T„  sin  0T  ) 


*33  F1 


h 


External  stores,  rough  air  or  landing  gear  conditions  are  not  expressed 
in  the  above  equations. 
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ROTATING  THRUST 

The  Hawker  P.1127  will  be  used  as  an  exanple  In  the  development 
of  rotating  thrust  V/STOL  simulation  equations  of  motion.  The  equations 
of  Section  III  as  well  as  information  from  Hawker  Aircraft  Ltd,  (Refer¬ 
ence  9)  will  be  directly  applicable  to  the  development  of  the  simulation 
equations. 

In  order  to  obtain  some  idea  of  the  P.1127  let  us  again  consider 
Figure  5  and  Table  8.  In  Figure  5,  a  three-view  arrangement  of  the  P.1127 
is  shown,  whereas  Table  8  contains  some  of  the  physical  characteristics 
of  the  aircraft. 

With  this  general  idea  of  the  P.1127,  we  may  now  develop  the  mathe¬ 
matical  model  for  the  P.1127.  First,  we  define,  if  necessary,  applicable 
axis  systems.  Second,  we  develop  any  additional  aerodynamic  coefficients. 
Finally,  we  write  the  equations  of  motion  for  the  P.1127. 

AXIS  SYSTEMS  FOR  THE  P.1127.  The  argument  presented  here  is  similar  to 
that  of  the  fan-in-wing.  The  only  physical  change  in  the  configuration 
of  the  P.1127  from  hover,  through  transition  to  normal  flight  is  the 
rotation  of  the  thrust  nozzles.  This  rotation  can  be  adequately  describ¬ 
ed  in  the  aircraft  body  axis  system.  In  accord  with  Reference  9,  body 
axes  can  be  used  during  hover  flight  and  stability  axes  can  be  used 
during  transition  and  conventional  flight  in  order  to  specify  stability 
and  control  derivatives.  Consequently,  no  axis  systems  in  addition  to 
the  inertial,  body,  stability  or  wind  axes  will  be  necessary  to  develop 
equations  of  motion  for  the  P.1127. 

AERODYNAMIC  FORCES  AND  MOMENTS  -  P.1127.  In  Reference  9,  Hawker  Aircraft 
lists  equations  of  motion  for  hover  and  transitional  flight.  During 
hover,  control  is  maintained  by  the  reaction  controls  which  obtain  con¬ 
trol  power  from  the  cockpit  stick  and  rudder  pedals.  There  are  two  roll 
reaction  controls  (one  near  each  wing  tip)  which  operate  in  conjunction 
with  the  ailerons  from  command  signals  generated  by  the  cockpit  stick. 

Next,  there  are  two  pitch  reaction  controls  (one  under  the  nose  and  one 
under  the  tail  of  the  aircraft)  which  operate  in  conjunction  with  the 
all  movable  tail  plane  from  command  signals  generated  t>y  the  cockpit 
stick.  Finally,  there  are  two  yaw  reaction  controls  (one  on  each  side 
of  the  tail  of  the  aircraft)  which  operate  in  conjunction  with  the  rudder 
from  command  signals  generated  by  the  rudder  pedals. 
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I 


’.'.'eirht 


I 


f 


Dasi  c 

Design  Cross 

V-1  yir, 

10,280  lbs 
15,719  lbs 

‘■'V  1  11^ 

Span 

Incidence  to  fuselage  datum 

Dihedral  Angle 

Aerocftmamic  Mean  Chord 

186. hO  ft.2 
22.82  ft. 
1.75° 

-12° 

8.667  ft. 

Horizontal  Tail  (all  movable) 

Area 

Span 

Dihedral  Angle 

Uh.20  ft.2 
12,00  ft. 
.if.  SO' 

Fuselage 

Overall  Length 

Overall  Height 

1*2.00  ft. 
10.75  ft. 

Vertical  Tail 

Area 

Span 

Height  from  fuselage  datum 

26.10  ft.2 
5.625 

8.167  ft. 

Underfin 

Area 

Distance  of  tip  from  fuselage  datum 

6.0  ft.2 
1.67  ft. 

Flaps 

Area 

Movement 

13.25  ft.2 

0  to  50° 

Reaction  Controls 

Distance  of  roll  reaction  nozzles  from 
fuselage  centerline 

11.20  ft. 

Distance  of  forward  pitch  reaction  nozzle 
from  moment  reference  center 

17.55 

Distance  of  rear  pitch  reaction  nozzle 

20.55  ft. 

from  moment  reference  center 

Table  8,  Selected  Physical  Characteristics  P.1127 
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Reaction  Controls  (Cont'd.) 

Pi  stance  of  yaw  reaction  nozzles  from 
moment  reference  center. 

Aerodynamic  Controls 

Aileron  Area  (Port  and  Starboard) 

Movement 

Rudder  Area 

Movement 

Tailplane  Net  Area 

Movement 

Trim  Range 

Intake 


Area 

Powerplant 

No.  and  Model 

Type 

Manufacturer 
By-Pass  Ratio 
Length 
Diameter 

Maximum  Engine  Rating 
Maximim  Fan  RPM  (%  design) 


20.15  ft. 


8.75  ft 
♦  126 
“5.25  ft 


2 

2 


♦  15°  j, 

38.08  ft* 

+  120 f  .  10o 

♦  5  1/2°,  -  2  1/2° 


9.3  ft.2 


(1)  Pegasus  B.S.  Pg.  5 
Ducted  fan  lift/thrust 
engine 

Bristol  Siddeley 

1.U 

8.25  ft. 

U.00  ft. 

18,000  lb. 

101.0 


Table  8.  Selected  Physical  Characteristics  P.1127  (Cont'd.) 
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From  Reference  9  further  information  concerning  thrust  is  obtained. 
There  are  four  jet  nozzles  operating  from  the  engine.  These  nozzles  are 
rotated  simultaneously.  The  front  port  and  starboard  nozzles  are  cold 
.lets  while  the  aft  port  and  starboard  nozzles  are  hot  jets. 

The  cold  jets  are  oriented  5°  out  from  the  x-z  body  axis  plane. 

The  hot  jets  are  a  function  of  Pj  *nd  are  given  by  the  following 

o  xt- ression:  ^ 

%  -  7.?°  ♦  0.0?  0^  (|u202) 

The  total  thrust  vector  (four  nozzles)  is: 

T  "  TGC  cos  ^  +  T0H  cos  *H 


T  «  0,996  Tgc  ♦  Tqh  cos  <f»K 


CU.203) 


where  is  the  gross  thimst  from  tre  two  cold  jet  nozzles  and  Tqh  is 

the  gross  thrust  from  the  two  hot  jet  nozzles.  We  will  turn  to  aero¬ 
dynamic  forces  and  moments  duri eg  hover  and  transitional  flight. 


Hover  Flight.  The  equations  used  by  Hawker  in  Reference  9  for 


hover  are  next  stated.  Nomenclature 
to  agree  with  usage  in  this  report. 


H  ■  S.  T  cos  (6j  +  a^)  -  g  sin  0 


has  been  changed  in  these  equations 

CU.20I*) 

(U.205) 

(U.206) 
(h. 207) 


V  -  g  cos  0  sin  4> 

W  « ■  S.  T  sin  (0j  +  dj.)  +  g  cos  0  sin  <t> 


a. 208) 

(U.209) 


In  the 


above  equations 


P 


and  r 


«E 


» 


and  Np  are  aerodynamic  moments  due  to  reaction  controls  where  A  is  the 

aileron,  E  is  the  tailplane  and  R  the  rudder.  6.  is  the  angle  of 

J 
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deflection  of  the  engines  nozzles  from  the  centerline  of  the  engine, 

aT  is  the  angle  between  and  the  x  bocfy-  axis  and  T  is  engine  thrust. 

Equations  (U.20U)  through  (U • 209 )  will  be  accepted  as  hover  equations 
and  incorporated  in  the  final  equations  of  motion. 

Transitional  Flight.  Again  we  turn  to  Reference  9  and  quote  here 
the  transitional  flight  equations.  These  equations  will  be  stated  in 
the  nomenclature  used  in  this  report.  Appropriate  comments  will  be  made 
for  each  equation  so  as  to  compare  it  with  the  math  model  of  Section  III, 
There  are  no  explicit  thrust  terms  in  the  X  and  Z  forces.  The  thrust 
terms  will  be  in  the  final  equations  —  Ch.216)  through  (U.<21). 

X  -  Force. 

mU  -  qS  [C  ♦  C  +  C  .  6F]  -  mg  (6)  (U.210) 

*U  a  6F 

Comparing  equation  (U.210)  to  (3.7),  there  is  no  Cx(p)  coefficient 

in  (U.210).  In  (U.210),  C  ♦  C  corresponds  to  C  (a,  Ma)  of  (3.7) 

*U  xa  x 

and  mg  sin  0  ^  mg  (0)  where  0  is  a  small  angle, 

Y  -  Force. 

mV  »  qS  [C  .  6  ♦  C  .  5A  ♦  C  .  6R]  ♦  mg  ($)  (U.211) 

y8  y6A  y6R 

Comparing  equation  (U.211)  to  (3.8),  there  is  no  C  of  C 

y6R 

in  (3.8).  In  jet  aircraft  C  and  C  are  usually  small.  They  are 

y6R  y6A 

included  in  (U.211)  because  of  their  importance  in  the  transition 
region.  The  term  mg  cos  0  sin  ♦  »  mg  ($)  where  0  and  $  are  considered 
small  angles. 

Z  -  Force. 

niW  »  qS  [C  +  C  +  C  •  6F]  +  mUq,  (U.212) 

ZU  Za  Z5F  1 

Conparing  (U.212)  to  (3.9),  there  is  no  C  .  6E  term  in 

Z6E 

(U.218).  This  term  should  be  in  the  final  equations  since  6E  does 
affect  the  Z  -  Force.  The  term,  C  +  C  corresponds  to  C  (a,  Ma). 

*TT  Z 


In  the  moment  expressions  Hawker  has  assumed  that  cros|-coupling  of 
angular  rates  is  negligible  so  the  terms  like  pq^^  0  and  r  0. 
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P  - 


21 

[11 


VBSb 

P  2— [C£  •6+C£c.  -6A  +  C. 
6 


6R 


•5A 


6R 


+  w:cl  •  p  *  &:ci  •  r] 

“  p  D  J7 


(U.213) 


Equation  (U.213)  is  the  same  as  equation  (3.10)  if  cross-coupling  of 
angular  rates  is  neglected. 


<Pitch. 


q-i *  P  57—  [C  +  C  ♦  C  .  5E  ♦  xg-  C  .  q, 
1  2I22  ma  2VB  mq,  1 


*  \  • i] 


(U.2Ui) 


Yaw. 

l33 


VBSb 

P57—  •  8  ♦  Cn 

2I33  n0  n6A 


6A  +  C  .  6R 
n6R 


+  cn  •  P  ♦  2I-  cn  .  r] 


2V 


B  p 


'B  “r 


(U.215) 


Equation  (U.215)  is  the  same  as  equation  (3.12)  if  cross-coupling 
of  angular  rates  is  neglected. 

From  the  hover  equations  and  transitional  equations  final  force  and 
moment  expressions  will  be  written. 

EQUATIONS  OF  MOTION  -  P.1127.  The  equations  of  motion 

herein  presented  are  written  by  considering  the  equations  U  .7  )  through 
(3.12)  and  the  Hawker  hovering  and  translational  equations.  In  equations 
(3.7),  (3.8)  and  (3.9),  the  velocity,  U,  is  predominant.  The  equations 
are  written  in  bocfcr  axes. 


X  Force  Equation. 


m  (U  +  Wq.  -  Vr)  -  o  [C  ♦  C  ♦  C, 
1  *  *U  *a 


.  6P  ♦  C  (0)] 
X5F  * 


♦  T  cos  (9j  *  a^)  -  mg  sin  0 


(U.216) 
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Y  Force  Equation 


* 


m  (V  +  Ur  -  Wp)  -  p-£-  [c  .  6  ♦  .  6A  +  Cv  .  6R] 

‘  y i 


6 


6A 


6R 


♦  mg  cos  0  sin  <f> 


(It  .217 ) 


Z  Force  Equation, 


m 


(w  ♦  v  -  uqx)  “ 


2 

0  -£ —  [G  +  C  +0  .  6F  ♦  C 

2  ZU  za  Zf>F  6E 

-  T  sin  (0j  *  aT)  ♦  mg  cos  0  cos  $ 


6E] 


(U.  218) 


Roll  Equation, 
I 


“  *99  lA  VjR> 

p  -  =-=■  (r  ♦  pq, )  -  (-^s - =0q,r  ♦  —  [C£  .  0  ♦  Ct  .  6A 

*11  1  *11  1  hi  11  6  6A 


♦  C,  ,  6R  ♦  *k—  •  p  ♦  jw—  C  ,  r]  (U, 219) 

«  £Wp  i 


B  p 


B  1  r 


Pitch  Equation. 


q  -  hi  (r2  -  p2)  -  (hLlii)  pr  ♦  •  E+PI§^-  [C  ♦  C  ♦  C  •  6E 

ql  t22  P  1 22  *22  **22  "V  ma  "SE 


6F 

Yaw  Equation 
'13 


♦  C_  ,  6F  ♦  =£-  C.  .  q,  ♦  55-  C  .  a] 
2Vo  m„  1  2Vg  mQ 


'B  q 


*22  ’  *11 


*R 


(U.220) 


r  •  =-^  (p  -  qr,)-  ( - - - )  P  qx  ♦  = —  •  R  ♦  r= —  [C  • 

*33  1  *33  1  *33  2I33  n8 


♦  C  .  6A  ♦  C  .  6R  ♦  C 
n6A  n6R  2VB  "p 


P  *  J*T"  •  H 
nr 


(U.221) 


For  W£?  0,  V  0  and  a  snail  value  of  U,  equations  (it, 216)  through 
(U.221)  are  equivalent  to  the  hover  equations  (U.20U)  through  (U.209) 
since  all  the  aerodynamic  coefficients  in  brackets  are  practically  sero. 
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As  H  becomes  greater  (6j  approaches  zero)  L^,  Mg,  and  NR  phase  out  and 

the  aerodynamic  terms  are  predominant  as  inequations  (U.210)  through 
CU.215).  The  transitional  equations  then  carry  directly  into  conventional 
flight*  External  stores,  rough  air,  or  landing  gear  conditions  are  not 
developed  in  equations  (£.216)  through  (li.221). 
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DERIVATIJr  OF  BASIC  EQUATIONS  OF  MOTION 


The  derivation  of  the  equations  of  notion  for  un  aircraft  flying  in 
a  conventional  manner  ha 3  been  done  by  means  of  a  notation  that  is 
economical  of  space.  Naturally,  there  is  a  sophisticated  degree  of  ab¬ 
stractness  in  manipulating  equations  developed  in  thi3  notation;  however, 
no  more  than  the  simplest  best  known  expressions  in  mechanics  and  the 
representation  of  direction  cosines  are  actually  being  used.  The  idea 
being  developed  in  an  equation  can  be  rerceived  as  a  whole  U3ing  this 
notation  instead  of  seeing  bits  and  pieces  of  derivation.  The  main  terms 
of  an  equation  are  immediately  apparent,  and  initial  assumptions  as  to 
importance  of  various  terms  in  the  equation  can  be  applied  in  order  to 
give  an  orderly  routine  of  simplification  of  the  general  equation  for  a 
%  part icular-se^of'assumpt ions. 

NOTATION 

The  notation  employed  uses  the  device  of  subscripting  and  super¬ 
scripting  a  particular  symbol  with  an  algebraic  variable  which  will  vary 
over  a  t rescribed  range.  It  is  actually  the  range  and  summation  tech¬ 
niques  developed  in  Tensor  Analysis.  More  precisely,  what  is  used  and 
manipulated  are  Cartesian  Tensors.  These  are  particularly  simple  tensors 
in  three  dimensions  using  Cartesian  coordinates.  In  fact,  for  the  purposes 
of  these  derivations  the  only  tensor  quality  is  the  adaption  of  the  range 
and  summation  conventions  plus  a  few  manipulative  methods,  and  the  advan-i 
tage  gained  from  Tensor  Analysis  is  a  concise  notation. 

e  -  SYMBOL 

or  eat)C  is  defined  in  the  following  manner* 

■  *1  if  ijk  is  an  even  permutation  of  12} 

■  -1  if  ijk  is  an  odd  permutation  of  123 


eijk 

eyk 
6  ilk 


°ijk  • 0  u  1  •  J 
For  example: 

*  e123  *  e23l 
-1  -  e132  -  e321 

0  *  *112  "  e122 


or  j  *  k  or  i  ■  k  o»  i  •  j 


"  *312 
“  *213  . 

’  *121  **111 


k 


93 


Reproduced  From 
Best  Available  Copy 


NATTHADaVCEN  120^-2 


KH  NECKER  DELTAS  (ft) 

F.r.it  Order  Kronecker  Deltas 

ft*  =  -1  if  i  =  j 
j 

ft  j  -  0  otherwise  (i  4  .1 ) 


Second  Order  Kronecker  Deltas 


6'lb  - 
i.i 

>1 

if  a 

4  b, 

i  4 

j.  a  -  i,  b  =  j 

,ab 

O  ,  .  ~ 

ij 

-1 

if  a 

4  b> 

i  4 

j  ,  a  =  ,i ,  b  =  i 

ab 

ft .  .  * 
VI 

0 

otherwise 

This  scheme  can  be  generalized  to  nb  order  Kronecker  Deltas.  How¬ 
ever,  since  the  coordinate  systems  being  used  in  this  report  are  three 
dimensional,  the  highest  order  Kronecker  Delta  encountered  will  be  three. 

Third  Order  Kronecker  Deltas 

ftab£  ■  *1  if  a  f  b  ^  c,  a  f  c  and  ijk  is  an  even 
J  permutation  of  abc. 

ft...  *  -1  ifa4h4c»a4c  an^  ijk  is  an  odd 
permutation  of  abc. 

^abc  ^  q  otherwise 
Ljk  * 


Product  of  Two  e  Symbols 


eabceijk 


abc 


In  particular  a  special 

e  e  -  ftaJc  -  ftac 
ajc  ijk  ijk  ik 


case  is 


Example* 

no«  ex23e3?1 


6 


13  , 
31  ' 


(♦1)  (-1)  «  -1 

-1  *  -1  ^.E.D. 
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Expansion  of  Second  Order  Kronucker  Delta 


6a6b  -  .«,a0 
1  J  J  1 


Example: 

12  12  12 
Does  6^  *  6^1  "  6i^2  ? 

(-1)  -  (0)  (0)  -  (1)  (1) 

(-1)  =  (-1)  Q.E.D. 


Example: 

,11 


*1,1 

*2*3 


,1,1 

6362 


23 

0  -  (0)  (0)  -  (0)  (0) 
0  =  0 


e  -  Symbol  Time?  Kronecker  Delta 

ca 

6,e.  ,  ■  e, 
d  jmd  jma 

Example: 

6lejml  *  ejml 
But  6^  =  »1 

6jml  "  ejml 
Example: 

62ejm2  =  ejml 
Possibilities 

62®312  "  ®3H 
0-0 


6*  “  0 


e. 


31 


or 

6?®1 <2 
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RANGE  AND  SUMMATION  CONVENTIONS.  The  definitions  are  as  follows* 

Range  Convention.  When  an  uncapitalized  letter  (superscript  or  sub¬ 
script)  occurs  unrepeated  in  a  term,  it  is  understood  to  take  all  the 
values  1,  2  and  3,  where  3  is  the  number  of  dimensions  of  the  space;  i.e„, 
three  dimensional  space. 

Summation  Convention.  When  an  uncapitalized  letter  (superscript  or 
subscript)  occurs  repeated  in  a  term,  then  summation  with  respect  to  this 
letter  is  understood,  where  the  extent  of  summation  is  1...3. 

For  example,  consider  the  following  expression* 

-  M(VX‘)  ♦  VxkcoXCe,  .  ) 

J  CK 

In  F3^ ,  the  j  indicates  a  ranges  j  *  1,  2  or  3 


So  there  is: 

F*1  -  M(VX1  ♦  VxkojXCelck) 

F*2  -  M(VX2  ♦  VXVCe2ck) 

F*3  -  M(V*2  ♦  Vxko)XCe3ck) 

Notice  now  VJtkuxc.  In  f*1,  rx?  and  F»  there  are  two  repeated 
indices— c  and  k.  This^Sould  indicate  a  summation  process  according  to 
the  summation  convention* 

F*1  «  M  (V*1  ♦  VXko)XCelck) 
c,k«l 


F*2  -  £  M  (V*J  * 
c,k»l 

F»  ■  £  M  *  A«e3ck) 
c,  k-1 

Expanding  F*1  ,  it  is  written* 

-(v*1  ♦  ♦  Mu2  ♦ 

A  .XL  X2  „X2  X2  .  „X3  X2 

*  *121  +  ^  w  *122  *  ^  w  *123 
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Collect  non-vanishing  terms  so  that 

F*1  .(?'  .  VXi,«  12J  •  /2„X,e1J2)M 

F*1  -(v*1  .  PKP  .  M 


F*2  and  F*3 


are  treated  in  a  similar  manner. 


In  the  above  expansion  of  F^  ,  the  summation  and  range  conventions 
were  used.  It  should  be  noted  that  there  is  an  easier  method  to  arrive 
at  the  expanded  FX1. 


F*1  -  M(VX]  ♦  vXka)XCelck) 


Note*  excj{  to  be  non  zero  can  only  have  c  and  k  equal  to  2  or  3.  That  is* 
e]xk  ‘  e123  ,nd  e132 

Then  for  F11  we  write  the  ncr-vanishing  terms  immediately* 

F*X  -  M(VX1  •  •;XJuX2ei23  ♦  VJt2uX,e132) 

F*1  -  (V*1  .  VX2»X3)M 


Thus,  by  noting  the  non  zero  properties  of  the  e  -  symbol,  the  sum¬ 
mation  process  can  be  shortened. 


CONTENTS  OF  NOMENCLATURE.  This  is  a  description  of  the  more  important 
individual  designations  that  will  appear  in  the  terms  of  the  equations. 


Position  -  Particle  to  Particle.  (1)  This  is  designated  by  a 
vector  r.  In  the  indexed  notation  ?  becomes  r^i,  which  indicates  that 
t  originates  at  the  origin  of  the  t-axes.  There  may  be  other  vectors 
such  as  r3^-  or  rn^. 

(2)  i  »  1,  2,  3  gives  the  particular  component  of  r*1*. 

(3)  To  differentiate  between  various  r*'*,  a  subscript  is  added. 

There  results  r^  ,  r^,  .  .  .  .  .  or  for  a  particle*  r^i  . 

Velocity.  This  is* 


Vti  - 
1  rl 


and  for  a  particle  V^ 

P 
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All  dots  over  symbols  refer  to  derivatives  with  respect,  to  time. 

Remember  that  ti  can  be  replaced  by  yi ,  zi,  si  or  some  other  set  of 
coordinates. 

Acceleration.  This  is* 

si  =  uai  t-si 
a  “V,  *  r 

o  1  o 


a_Si  -  Vsi 
*1  1 


i31 


and  for  a  particle  aS^  *  Vs*  **rsi 

p  p  p 

Mass.  A  mass  point  is  designated  as  m  .  When  deriving  equations  of 
motion  a  summation  of  mass  points  is  made.  ^ 

Example* 

For  fuselage  particles  /m  ■  Mp(mass  of  fuselage) 

Force.  In  general*  F  *  m  a  (vector  form) 

Now  in  indexed  notation  for  example*  f*’1  *  m  a^1 

P  P  P 

But  the  equations  are  developed  by  knowing  what  particles  are  being 
manipulated.  If  the  expression  for  the  fuselage  particle  is  being 
sought  the  final  form  will  be* 

f*’*  *  m  a*'*  ,  where  a*1*  is  acceleration  and  is  a  function  of 
P  P  P  P 

direction  cosines,  position  vectors,  and  rates  of  change  of  both. 

Integrating*  /f^  *  Fp  Subscript  F  fuselage 

/■V  "  «F  *F=  Mass  of  fuselage  particle 

Torque.  Vectorial*  -  "r  x  "? 

New  in  indexed  notation*  tSm  ■  m  as^rsde.  ,  ■  fs^rsd 

p  p  p  jmd  p  p 

Integrating*  /tSm  *  Tpm  if  is  the  acceleration  associated 
with  the  fuselage  particle.  p 

Inertia,  Vectorial*  The  angular  mementum  J  -  ^  ^miri  x  (“  x  ri) ] 


or 


7  -  T[“x  £  mt(y^  .  ,*)  -  ^ 

*  -  “x  t-1  Wt  *  “y  £  *  zl>  - 

*  '  “x  A  Vt't  -  “y  "iyi2i  *  w  ”i<xi  *  y?)) 
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Momenta  of  inertia  are  terms  1  ike: 


<r  ,  2  2, 

lVyi  +  zi>  ■  I 


XX 


Products  of  inertia  are  terms  like: 

Zm.x.y,  =  I 
i  ri  xy 

Now  in  indexed  notation:  I 


pq 


r  xa  xd_ 

/mprp  rp  epakeqtlk’  P,  a  fixed 


r  X8  xd  .  T 
or  ym  r  r  *  -I 
P  P  P 


ad  4  ^Xll  4  X22  4  X33)6 


I  defines  both  moments  and  products  of  inertia.  To  make  tnis  clear 

pq 

consider  the  following  examples: 

T  r  xa  xd  r  xa  xd.qd 

I  *  /m  r  r  e  .  e  ■  ym  r  r  6^ 
pq  J  p  p  p  pak  qdk  P  p  p  pa 

Example:  Moment  of  inertia 

-  -  U  *a  xdld 

*  J m  r  r  6, 

11  p  p  p  la 

Now  6Xd,  6X?t  and  6X^  are  non  zero 
la*  12  *  13 


Then: 


*11 

f  X2  X2  12 
*  y  m  r  r  6.  _ 
PP  p  12 

♦  A 

.12 

612 

-  >1,  6 1 j  - 

+1 

*11 

„  /  ,,  X 2*2 
-/mp[(rp  ) 

> 

■  + 

(r« 

P 

Example: 

Product  of 

inertia 

/•  xa  xd. 

.2d 

X12 

■  y  m  r  r  f 
P  P  P 

’la 

X3  X3.13 


2d  21 

Now  6la  -  6XJ  -  -1 


Then: 


T  r  X2  XI 

*12  '  --Vp  rp 


E:  '"y  j  f-i 

1-  v  t . 


a  ;;d 

r 


T. 


<a  *o 
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Let  a  =  1  and  d  =  1 

/•  XI  XI  .  1,_  T  T 

ymprp  rp  ■  -hi  *  ?(I11  *  hi  *  hi'6! 


6^  -  ♦! 


'll 

,  n  *2>2 
‘  >pt  <  rp  > 

*  (■/);>] 

X22 

r  r  /  X3  v  2 

•■/y(rp  > 

*  (rf)2] 

*33 

XI  2 

»/mp[(rp  ) 

*  <rpVi 

So  that* 

r  „X1  XI 

>Prp  rP 


.  fu  ♦  !§*  *!gi 


.X3  X3 


r  XI  XI  ,  X2  X2  r  AJ  / 

/m  r  r  *  -  /in  r  r  -  /m  r  r 
PPP  PPP  ppp 


X3.X3 


*  /nvYrp  *  /mprp 


/mrfrf 


/■  n  xi 

♦  /  m  r  r 

ppp 


*  ^ 


X2  X2 
r  r 


AXES  SYSTEMS 


Two  different  axis  systems  are  described — the  inertial  system  and 
the  fuselage  system  in  terms  of  direction  cosines.  Direction  cosines 
are  the  angles  which  define  a  vector  in  a  particular  axes  system.  The 
matrices  developed  in  this  appendix  relate  the  measurement  of  a  unit 
length  in  one  direction  as  measured  in  a  different  direction.  This  is 
accomplished  by  the  reorientation  of  one  axis  system  with  respect  to 
another  axes  system.  The  method  used  is  to  consider  that  the  axes 
systems  have  conmon  origins  and  then  by  rotating  the  axes  of  one  system 
with  respect  to  the  other  achieve  coincidence.  The  direction  cosines 
measure  this  rotation. 


In  this  report  objects  in  three  dimensional  space  are  described  by 
the  use  of  Cartesian  coordinates  which  are  orthogonal  by  definition. 
Orthogonality  allows  the  transpose  of  the  direction  cosine  matrix  to 
equal  its  inverse.  It  is  sufficient  to  describe  the  rotation  of  one 
axes  system  with  respect  to  another  by  the  use  of  three  angles  (Euler 
angles). 
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INERTIAL  AXES  SYSTEM.  The  system  is  identified  by  Na  axes  (a  ■  1,  2,  3). 

N1  is  in  the  horizontal  plane  with  positive  direction  North. 

N2  is  in  the  horizontal  plane  with  positive  direction  East. 

N3  is  normal  to  the  N1N2  plane  with  positive  direction  down. 

THE  FUSELAGE  AXES  SYSTEM.  The  system  is  identified  by  xa  axes  (a  «  1,2,3). 
This  axes  system  corresponds  to  the  body  axes  system  used  in  normal  air¬ 
craft  parlance.  The  origin  is  fixed  in  the  fuselage  (usually  at  the 
nominal  center  of  mass).  The  origin  of  the  x-axes  is  separated  from  the 
origin  of  the  N-axes  by  a  vector  rQ — see  Figure  Al. 

xl  is  in  the  horizontal  plane  with  positive  direction  noseward. 

xl  is  not  necessarily  parallel  to  the  water  line  axis  defined 
by  the  aircraft  manufacturer. 

x2  is  in  the  horizontal  plane  with  positive  direction  right. 

x3  is  normal  to  the  xlx2  plane  with  positive  direction  down. 

INERTIAL  TO  FUSELAGE  AXES  DIRECTION  COSINE  MATRICES.  Three  angles,  the 
Euler  angles  are  sufficient  to  locate  the  fuselage  in  inertial  space. 

The  inertial  system  is  translated  such  that  it3  origin  coincides  with 
the  origin  of  the  fuselage  (i.e.  the  nominal  c.m.  of  the  aircraft). 

Specific  rotations  of  the  N-axes  are  then  made  so  that  a  unit  of 
measure  in  the  N-axes  can  be  expressed  in  the  x-axes.  In  each  case  the 
particular  rotation  is  in  a  plane  so  that  there  is  a  rotation  about  the 
third  axis— that  is  the  rotation  written  as  a  matrix  is  simj  'j  sines 
and  cosines  of  the  rotating  angle  with  the  orthogonal  propel  c y  of  the 
axes  systems  stating  that  sin2  ♦  cos2  »  1.  This  will  become  clear. 


The  three  rotations  are  ijr,  9  and  $. 

Yaw  by  Angle  t  Around  N3  Axis.  Here  an  intermediate  system  N1  is 
defined  by  the  rotation  of  t  around  W3. 

Define  C  as  direction  cosine  matrix. 
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N3 

Down 


Figure  Al.  Inertial  and  Fuselage  Axes 


» 
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Pitch  by  Angle  0  Around  n'2  Axis.  Here  an  intermediate  system  N" 
is  defined  by  the  rotation  of  0  around  N  2. 


N  "l 

N  "2 

N”3 

N'l 

COS0 

0 

sin  8 

-  N'2 

0 

1 

0 

N*3 

-sin0 

0 

COS0 

Roll  by  Angle  *  Around  N  nl  Axis.  This 
the  x-axes. 


Combination  of  Three  Rotations  (t,  0,  #).  To  get  the  final  result— 
the  relation  between  the  inertial  coordinate  system  and  the  fuselage 
coordinate  system— a  matrix  multiplication  is  performed. 


Ihis  multiplication  is 


cos  t 

-sint 

0 

sinf 

cost 

0 

C 

0 

1 

cost 

-sint 

0 

sint 

cost 

0 

0 

0 

1 

COS0 

0 

sin0 

0 

1 

0 

-sin0 

0 

COS0 

1 

0 

0 

0 

cos* 

-sin* 

0 

sin* 

cos* 

COS0 

sine  sine 

sine  cos* 

0 

COS* 

-  sin* 

-sin® 

cose  sin* 

cose  cos* 
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xl 

x2 

x3 

cost  cos0 

cost  sin0  sin* 

cost  sin0  cos* 

-sinf  cos* 

♦sint  sin0 

sint  cos0 

sint  sin©  sin* 

sint  sin©  cos* 

♦cost  cos* 

-cost  sin* 

-sin0 

cos0  sin* 

cos©  cqs4> 

Note  that  cjj®  cjj£  -  1 

where  cj|j  -  (C^)"1  -  (C^)T 

(  ) Inverse  matrix 

T 

(  )  Transpose  matrix 

IN BIT I AL  TO  FUSELAGE  AXES  RATES.  *,  0  and  t  are  measured  as  p,  and  r 

around  the  xl,  x2  and  x3  axes  respectively.  The  dot  is  the  total  deriva¬ 
tive  with  respect  to  tljne. 

t  Occurs  Around  the  N3  Axis.  To  get  the  direction  cosine  components 

refer  to  . 

*P 

p  ■  ♦(-sin9) 


q^  ■  t  cos0  sin* 
r  -  t  cos®  cos* 


&  Occurs  around  the  N'2  Axis, 
refer  to 


y*  xy )  .cn> 

x  °xp5  xp 


p  -  0  (0) 


To  get  the  direction  cosine  components 


q^-  0  cos* 
r  -  0  (-sin*) 


*  Occurs  Around  the  N,"  axis, 
refer  to  cN  *• 

P  -  *  (1) 

♦  (0) 

r  -  ♦  (0) 


To  get  the  direction  cosine  conqjonents 
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Combination  of  Thraa  Kata*  (p,  q1 ,  r). 

p  »  f  (-sind)  *  #  (0)  ♦  ♦  (1) 

•  •  • 

q^  »  t  cosQ  sin#  ♦  t  coa#  ♦  ♦  (0) 

•  *  • 

r  *  t  co s0  cos#  ♦  0  (-sin#)  ♦  ♦  (0) 

In  matrix  notation  this  becomes 


DERIVATION  OF  FORCE  AND  TORQUE  BQUATI'JXS 

The  derivation  of  force  and  torque  equations  is  dene  In  conolee 
notation.  Certain  assumptions  are  present  In  the  use  of  these  aircraft 

equations  of  motion.  These  aasiuptlona  are  ae  follows* 

(a)  The  earth  is  assumed  to  be  a  flat  plane  with  gravity  acting 
normal  to  this  plane  and  it  ie  not  rotating  or  translating 
through  space.  In  the  derivation  of  the  squat lone  of  notion, 
Newton's  lav,  F  •  ma,  la  used  in  its  aoet  staple  form— 
there  la  no  mas a  variation  with  change  in  tine. 

(b)  The  xl  and  x3  fuselage  axes  define  a  plane  of  eye  try  far 
the  aircraft. 

(c)  There  is  a  fixed  center  of  gravity. 

(d)  There  is  no  relative  notion  between  the  earth  and  its 
atmosphere . 

The  equations  will  now  be  developed  using  the  inertial  and  fuselage 
axes,  Newton's  lav,  and  direction  cosines  of  the  Euler tan  angles . 
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ACSi.C  IATMD  rf'  TH  THE  F.JSELA  >S  I  ARTICLE.  The  position  of  m  for  the 
fusel  are  i. article  with  respect  to  Ni  frame  of  reference  is  statid  act 


N  i  N  i  x  a  „N  i 

r  -  r  +  r  C 

o  P  xa 


fee  Figure  Al . 

fhe  inertial  system  velocity  is  Riven  as: 


Ni 


.Ni 

.Ni 

♦  fx'1C1‘i 
p  xa 

xa5Nii 

*  r 

»  r 

0 

r  o 
p  xa 

.Ni 

xa 

iN  i  .  X  ri 

.  .  1 
0  since  r 

’ 

“  fo 

♦  r 

P 

u  ,r 
xa  >  p 

Ni 


constant 


,  .  _Ni  Ji  xc  . ,  „ 

but  C  «  C  ,  u>  e  ,  ,  therefore 

xa  xb  abc’ 


vNi 

p 


xar_wi  xc  i  ,  ,  ,  . .  . 

r  [  C  e  ,  j  by  substitution 
p  xb  abc 


(Al) 


Note:  For  vector  treatment  of  this,  refer  to  Reference  3,  The oretical 
Phys ics  -  Me chanjcs  by  F.  W.  Constant,  pages  69  to  71. 


The  inertial  system  acceleration  is  given  by  a  similar  procedure  as 
Ni  ’Ni 

“p  -VP 

aNi  ..fNi  ,  rxa  (0NCxc  .  jNi  xc  , 
p  o  p  xb  abc  xb  abc 

Ni  ..Ni  .  xa  r„Ni.xc  ^  nNi  xg  xc  , 
p  o  p  xb  abc  xh  bhg  abc 


Ni  _..Ni  .  xa-Ni.xc  ^  xa-Ni  xg.  xc 

r  *  r  C  e  ♦  r  eui_  e 

p  xh  w 


a 


(A2) 


(A3) 


(A3) 


p  o  Jp  “xb  v'nbc  *p  “xh  ~  “bhg  abc 

For  a  particular  mass  particle  mp  it  can  be  said  that 

-Ni  Ni 
f  *  m  a 
P  P  P 

and  transforming  to  the  x-axis 

f50  •  n  aN1C^  . 
p  p  p  Ni 

Substituting  (A2)  into  (A3)  yields 

rxl  -  m  rfN1  ♦  rxaCNHxce  .  -  rxaC "‘.VV  e  ')<$ 

p  pk  o  p  xb  abc  p  xh  bhg  abc  Ni 

Interchange  b  and  h  where  they  both  appear  in  the  same  term  and 
then  by  changing  b  to  j , 

?  ■  *  -JWc.hJge^c].  (AM 


I 


I 
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Now  is  written 

=  fNi 
o 

^xj  ,  fNi  x.j 
o  Ni 

XI  _.,NUxj  .Ni  xk  xc 
r  -  fo  '•Ni  ro  CNt"’  ,ikc> 

v**  -  fNic;k . 

o  Ni 

Substituting  ( A6)  into  (A£)  gives 

V*J.  V  NicJ'j  *  Vxk..ixce  . 

o  Ni  jkc 


(A5) 

(A6) 


Rearranging  and  letting  e.^c  =  ~e-c'K» 

.^Ni^xj  =  ♦  VXka)XCe.  ,  .  (A7) 

o  N  t  ,1  ck 

Substituting  (A?)  into  ( AJ4)  offers 

gJ  .  m  ♦  V*k»xbe.hl,)  .  rxVce  .  .  rJWc.  .  eh,  ]. 

p  p1'  Jbk  p  ajc  p  ahc  h.ig 

Y  -j 

Integrating  over  r  ,  rn^  and  m^r^  finally  yields. 

F'X'i  -  ♦  VXk-XC«Jck  ♦  *  V^—abcW- 


S  i  nee 


P  e  =  fih^  -  6hje  =  . 

'nhc  hjg  ahe  hca  ca  * 


f  -  •y**J  *  **•.«._> .  yjvc*,Jt  •  «i!)  • 


(AS) 


TOR ,tLJK  !MK  TO  Hit  FUSELAGE  PARTICLE.  The  definition  of  torque  offers: 

,xm_  xi  xd  pci  xd 

t  _i  m  a  r  e.  ,  =  r*r  e. 

V  ;  p  imd  v  r  .  r'd 


n  at  is,  oy  Knowing  the  force,  the  tor  rue  of  a  particle  can  be  developed 
oy  i  ib:;t  ituting  f  l  r  fx.i  the  terns  n  a*-  from  (A8),  real  iz  ing,  of  course, 

p 

♦  hat  this  a/'  is  the  acceleration  exrress ion  level  ore  d  for  the  fuselage 

rarticle.  Then  it  follows  from  tn.e  substitution  that 

**ri  _  _  „xcif/X’  .pK  xc  xa.xc  xa  xc  xgtjg, 

t  -  -i,  rr  [CT  ♦  r  x.  e<  ck 1  *  rp  *  eaic  ♦  *>  o  o  ^Je1md 


i  A9) 
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By  substituting  (A9)  into  the  expression 


.  xiti 

yt  =  T 


,xm 


(A10) 


There  is  obtained 

Tf  -  /m  rxd  (?3  ♦  VxXoixce .  „)e,  „ 

F  p  p  jck  jmd 

+  /m  rXarxd(J>Xce  .  e.  (All) 

P  P  P  ajc  jmd 

r  xa  xd  xc  xgEjg 
+  /n  r  r  «  w  F’6UDe.  , 

J  p  p  p  ca  jmd 

Consider  the  second  and  third  terms  of  the  expression  (All).  These 
will  be  developed  using  the  definition  of  inertia.  The  definition  is: 


Jpq  "  >prpaVWeqdk'  "Uh  P-  q  fiXed 

I  -  /b  rXW®  , 

pq  p  p  p  qd  » 


( A12) 

(A13) 


or, 


W  *  -rad  *  7(IU  *  rS2  *  V6*’ 


f  xa  xd.xc 

Am  r  r  <o  e.e,  ,  * 

J  p  p  p  ajc  jmd  pcm 


T  J.XC 

I__<5>  » 


(Aia) 

(A15) 


Now  is  written 

r  xa  xd  xc  xg,  j g 
Jm  r  r  w  w  ' b^e .  . 
J  p  p  p  ca  jmd 


*  (-1 


ad 


Vhl  *  *22  * 


(A16) 


and  knowing  that 


6dejmd 

6.1g 

ca 


'jma 
-  6^6g 


6Jb 
c  a 


6°6( 
a  c 


(A17) 


Substitution  of  (A17)  into  (A16)  gives 

f  xa  xd  xc  xg,  jg 
/m  r  r  w  B6“°e,  , 

J  p  p  p  ca  jmd 


m  a) 


x0.X*I_Je,_J(6;<6f  -  6^6«) 


ad  jmdv  c  a 

!<*11  4  J22  *  *»>  (B=6a  '  6i6f>  “XC“XS 


ma 


108 


NAVTRADBVCXN  1205-2 


^XCa>X«I  ,e 


xc  XCT 

w  ci)  I  ,e  , 
gd  cmd  ad  and 


♦  *  i„„  .  i 


xc  XgT 
-a)  oo  bI  ,e 


xc  xe  1/T 

22  T  "33  7  6cmg  "  “  “  *(I 


pgd  cmd  * 

Then ,  in  conclusion, 

?  •  *  7XVC8Jck)'j.d 

♦i  ®xc  .  „xV*x  i  _ 

£ca  pgd  cmd 


?KLll  *  l22  *  X33)# 


(Al8) 


(A1 9) 


which  is  obtained  by  substitution  of  (AI5)  and  (A18)  into  (All). 

KIPAMSIOK  OF  FORCE  AND  TORQUE  EQUATIONS.  Equation  A8  is  the  force  equation 
and  equation  A19  is  the  torque  equation.  The  terms  in  these  equations  are 
defined  as  follows* 

(a)  Mp  is  the  aircraft  mass. 

(b)  v*1  are  body  axes  velocities. 

(c)  is  an  angular  velocity  around  the  appropriate  body  axis. 

(d)  r~  is  the  displacement  of  the  center  of  gravity  fro* 
origin  of  the  x-axes „  Here  r^  “  0. 

(e)  F^  are  body  axes  forces. 

(f)  t“  are  body  axes  torques. 

Dropping  the  subscript  F  these  equations  are  then  written  ast 

,«J  .  h(?j  1 

T“  -  I  <3XC  -  u**I  ..  . 

cm  gd  cmd 

Expanding  theae  equations  there  1st 

F*1  -  MfV*1  ♦  /V*  -  v*2*/3) 

F*2  -  M(  v*2  ♦  vxlux3  -  y313^1) 


r13 

-  M(V*3  ♦ 

- 

^V2) 

r11 

■  hi 

4>xl 

*  hf' 

x2  xl 

♦  (1)  0) 

R) 

x2  x3t 

*  “  “  X33  ■ 

«x\%2  in 

T12 

•  X 

l22 

*x2 

xl  xl 

—  U)  uJ 

T  Xl 

rn  -  •’ 

x  1 

a> 

I3,  ♦  ux3axl 

ll1  .  »x\,x3 

T*5 

'  *33 

ciX3 

.  -  *xl 
^  A31  J 

x2  x) 

-  U)  03 

Si 

x2  xl  r 

-  “  “  3.1 

*  wxlcx2  Ijj 
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Since  the  aircraft  is  symmetric  with  the  xl  -  x3  plane  only  the  I,.,  products 
of  inertia  exist.  1^  carries  a  minus  sign.  ^ 

Finally,  by  using  more  conventional  nomenclature  these  force  and  moment 
equations  will  be  more  readily  recognized.  Fx^  ,  Fx  ,  Fx^  ,  Tx^  ,  Tx‘^  and 
T*3  are  the  total  external  forces  and  torques,  and  as  such  denote  aerodynamic 
effects  (also  gravity  and  thrust  effects  if  applicable). 


F*1 

F*2 

pX3 

TX1 

T*2 

T*3 

V*1 

yX2 

V*3 


x3 

O) 


Xa,  total  aerodynamic  x-force. 

Y  ,  total  aerodynamic  y-force, 

Zfl,  total  aerodynamic  z-force. 

L  ,  total  aerodynamic  rolling  moment. 

2i 

M  ,  total  aerodynamic  pitching  moment. 

d 

N  ,  total  aerodynamic  turning  moment. 

U,  velocity  along  x-body  axis. 

V,  velocity  along  y-body  axis. 

W,  velocity  along  z-body  axis. 

p,  angular  velocity  around  x-body  axis, 
q^,  angular  velocity  around  y-body  axis, 
r,  angular  velocity  around  z-body  axis. 


M 


m,  the  aircraft  mass. 


The  equations  for  the  basic  aircraft  then  become  the  following  expressions. 
Xa  -  m(U  ♦  Wqx  -  Vr) 

Ya  -  m( V  ♦  Ur  -  Wp) 

Za  -  m(W  ♦  Vp  -  Uqx) 

La  "  *11^  ’  X13(f  *  ^1*  *  (I33  I22)qlr 
M.  ■  *22*1  -  V*  -  p2)  *  du  -  I33>Pr 
Na  -  I33  *  -  I13^  "  V>  *  (X22  '  1ll^pql 
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